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EXPLANATION  OF  REVISIONS  TO  ASAP 

Since  the  issuance  of  the  original  ASAP  report  in  December 
19  74,  and  distribution  of  the  source  program,  feedback  from  many 
users  and  direct  assistance  from  Dr.  Robert  Bevensee  of  Lawrence 
Livermore  Labs,  University  of  California  enabled  the  preparation 
of  this  revised  version.  During  the  past  2  years,  no  additional 
discrepencies  have  arisen  and  it  is  now  felt  that  the  code  can  be 
safely  called  "revised". 

The  nature  of  the  improvements  was  as  follows: 

1.  To  correct  the  manner  in  which  structures  were 
treated  when  elevated  over  a  ground  plane. 

2.  Improvements  and  corrections  in  calculating  finite 
ground  plane  effects. 

3.  Text  corrections  in  explanation  of  ground  effects. 

4.  Instructional  changes  in  the  DESCRIPTION,  FREQUENCY, 
and  CHANGE  cards. 

5.  Corrected  sample  problem  outputs. 

The  user  should  note  that  some  facilities  will  experience 
printing  errors  when  NEAR  FIELDS  are  called  for.   Statement  58  in 
the  MAIN  program  should  be  suitably  rewritten  if  that  occurs. 


I«   INTRODUCTION 

Although  many  thin-wire  computer  programs  have  been 
developed  for  the  purpose  of  analyzing  antennas  and 
scatterers,  few  of  these  programs  have  been  directed  toward 
the  student  of  elect ro- magnetic  theory.  The  majority  of  the 
programs  are  directed  to  the  engineer  or  advanced  student 
for  the  purpose  of  analyzing  designed  structures  or 
verifying  experimental  data. 

The  purpose  of  the  study  is  to  develop  a  computer 
program  by  modifying  an  existing  computer  code  which  can  be 
utilized  as  an  educational  method  to  develop  insight  into 
radiating  structures  by  the  beginning  student  of 
electro-magnetic  theory. 

The  modified  Ohio  State  University  Antennas-Scatterers 
Analysis  Program  (OSUMOD  or  ASAP)  is  directed  toward  the 
beginning  student  who  does  not  yet  have  the  expertise 
necessary  to  manipulate  the  input  data  for  proper  execution 
of  the  larger  more  comprehensive  analysis  program.  Even 
though  ASAP  is  small  in  core  requirements  and  is  fast  in  run 
time,  it  is  capable  of  analyzing  structures  to  assist  the 
engineer  with  design  problems. 

Since  the  resulting  program,  ASAP,  is  primarily  directed 
toward  students,  the  program  has  been  limited  to  structures 
which  contain  less  than  50  monopoles  (segments) ,  no  longer 
than  one-fourth  of  a  wavelength,  and  which  have  less  than  51 
nodes  (intersections  and  endpoints) .   If  a  ground  plane, 
either  perfect  or  finite  is  present;  the  stated  limits  above 
are  halved  due  to  the  generation  of  an  image  structure. 


II«       ORIGINAL    PROGRAM 

A.  THEORY 

Reference  1  presents  the  electro-magnetic  theory  for  the 
analysis  of  antennas  and  scatterers  in  an  isotropic,  linear, 
and  homogeneous  ambient  medium.  The  analysis  is  performed  in 
the  freguency  domain  with  an  excitation  caused  by  either  a 
generator  or  an  incident  wave. 

In  the  analysis,  a  piecewise-sinusoidal  expansion  is 
used  for  the  current  distribution.  The  matrix  eguation 
Z  I  =  V  is  generated  by  enforcing  reaction  tests  with  a  set 
of  sinusoidal  dipoles  located  in  the  interior  region  of  the 
wire.  Since  the  current  distribution  has  the  same  form  as 
the  expansion  mode,  this  formulation  is  known  as  the 
"sinusoidal  reaction  technigue" . 

B.  COMPUTER  PROGRAM 

Reference  2  presents  the  computer  program  corresponding 
to  the  theory  presented  in  Ref.  1. 

1.  Input  Format 

In  the  program,  the  input  data  must  specify  the 
freguency,  wire  radius,  wire  conductivity,  the  parameters  of 
the  exterior  medium,  coordinates  of  the  points  to  describe 
the  shape  and  size  of  the  wire  configuration,  a  list  of  the 
wire  segments,  and  the  indicators  for  the  various  outputs. 
Table  1  is  the  input  data  necessary  to  analyze  a  half-wave 
dipole. 

2.  Output  Format 

In  the  original  form,  the  only  outputs  which  could 
be  reguested  by  the  input  data  stream  are  the  following: 


a.  Antenna  Problems 

(1)  Current  Distribution  on  the  Structure. 

(2)  Input  Impedance. 

(3)  Radiation  Efficiency. 

(4)  Near-Zone  Field. 

(5)  Far-Zone  Field. 

b.  Backscatter ing  Problems 

(1)  Absorption  Cross  Section. 

(2)  Scattering  Cross  Section. 

(3)  Extinction  Cross  Section. 

(4)  Complex  Elements  of  the  Polarization 
Scattering  Matrix 

c.  Bistatic  Scattering  Problems 

Echo  Area. 

Table  2  is  an  example  of  the  output  data  available  for 
data  of  table  1. 

3.   LIMITATION 

Although  the  program  can  analyze  a  structure  with  up  to 
50  segments,  55  points  and  60  dipoles  modes;  it  can  not 
analyze  a  structure  in  the  presence  of  a  finite  ground 
plane. 
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Ill*   MODIFIED  COMPUTER  PROGRAM 

A.  Input  Format 

As  illustrated  in  table  1  the  format  for  the  input  data 
cards  is  not  self  explanatory.  This  format  can  be  determined 
by  referring  to  the  FORMAT  statements  of  the  program  of  Ref. 
2.   Since  the  modified  program  is  directed  toward  the 
student,  the  input  data  format  was  changed  to  allow  free 
format.  Reference  2  was  written  in  a  form  which  permitted 
modifications  to  allow  flexibility  in  specifying  input  data 
for  the  analysis  program.   Appendix  B,  titled  "User's 
Manual",  discusses  the  input  data  cards  necessary  for  proper 
execution  of  an  analysis  problem.   Appendix  B  is  self- 
contained  and  may  be  used  independently  of  the  remainder  of 
this  document. 

B.  Output  Format 

In  the  original  computer  program,  the  absence  of  labels 
encumbered  the  output  data  and  lessened  the  usefulness  of 
the  program.  To  improve  the  usefulness  of  the  modified 
version,  detailed  labels  were  added  to  the  output  data.  As 
with  the  input  data,  Ref.  2  was  written  in  a  form  which 
enabled  modification  to  allow  more  specific  output  data  for 
the  analyzed  problem.  With  the  addition  of  the  polar 
plotting  package,  the  far-zone  electric  field  intensity 
polar  radiation  and  reradiation  patterns  can  be  plotted.  A 
sample  problem  can  be  found  on  page  120  in  Appendix  B, 
User's  Manual. 

C.  Finite  Ground 

To  enable  the  student  or  the  engineer  to  have  an 
improved  analysis  program,  the  finite  ground  effects  were 
added  to  ASAP.  The  theory  corresponding  to  the  ground 


11 


effects,  which  utilize  Fresnel  reflection  coefficients,  is 
discussed  in  Appendix  A,  titled  "System  Manual".   Also 
discussed  in  Appendix  A  is  the  modified  computer  program  and 
the  corresponding  theory.  The  electro-magnetic  theory  was 
developed  in  Refs.  1,  2,  and  3;  and  it  is  restated  with  its 
corresponding  computer  code  to  assist  in  the  understanding 
of  the  methods  applied.   Appendix  A  is  self-contained  and 
may  be  used  independently  of  the  remainder  of  this  document. 


IV*   CONCLUSION 

The  addition  of  ground  effect  techniques  to  the  original 
program  did  not  alter  the  accuracy  or  the  computational 
capabilities  of  the  program.  The  ground  effect  techniques 
utilized  the  results  of  the  original  program  and  modified 
these  results  to  account  for  the  effects  of  the  presence  of 
the  finite  ground. 

To  verify  the  numerical  results  of  ASAP,  the  input 
impedances  of  both  a  horizontal  and  a  vertical  dipole  were 
compared  to  the  solutions  of  the  exact  form  of  the 
Sommerf ield ' s  equation.   As  can  be  seen  in  table  3  the 
finite  ground  treatment  of  ASAP  agrees  favorably  with 
Sommerf ield' s  solutions.  The  ASAP  finite  ground  results  are 
also  in  excellent  agreement  with  the  previous  computer 
solutions  of  Refs.  4  and  5. 


V«   RECOMMENDATIONS 

Although  the  program  is  a  general  analysis  tool  for 
students,  several  future  modifications  will  enhance  the 
program  as  a  design  tool  for  engineers.  These  items  include: 
varying  the  wire  radius  on  the  structure;  incorporation  of 
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non-radiating  elements  such  as  transmission  lines;  varying 
the  wire  insulation  radius,  conductivity,  and  dielectric 
constant;  and  a  geometry  generation  package  such  as  dipole 
array  or  helix.  One  major  change  that  would  both  improve  the 
speed  and  reduce  the  core  requirement  is  that  of  symmetry. 
No  attempt  was  made  to  utilize  the  symmetry  in  the 
admittance  matrix  when  the  ground  plane  is  present.  If 
symmetry  were  applied,  the  structure  size  limit  with  the 
ground  plane  present  would  be  approximately  that  of  the 
structure  without  the  ground  plane. 


VERTICAL  DIPOLE 
FREQUEMCY    3MHZ 
LENGTH    .5  WAVELENGTH 
RADIUS    .005  METERS 
DIELECTRIC  CONSTANT  (RELATIVE) 
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CONDUCTIVITY    HEIGHT/WAVELENGTH 


.1 


,001 


.00001 


.25 

.30 
.35 
.45 
.25 
.30 
.35 
.45 
.25 
.30 
.35 
.45 


ASAP 

123. 750  68.30 
98.620  38.26 
87.600  35.64 
78.690    *1.79 

107. 180  55.36 
91.800  40.18 
85.150  39.30 
80.  lb* J    43.52 

103.200  57.99 
91.830  41.99 
85.900  40.14 
86.780    43.23 


EXACT* 


126. 50  83.89 
ID0.2O  49.52 
88.500  46.52 
79.210  52.60 
119. 40  71.46 
94.300  49.02 
85.410  48.93 
80.170  54.83 
115. lO  73.89 
94.090  50.69 
86.010  49.83 
80.720    54.63 


HORIZONTAL    DIPOLE 
FREQUENCY         3    MHZ 
LENGTH         .5    WAVELENGTH 
RADIUS         .001    METERS 
DIELECTRIC    CONSTANT    (RELATIVE) 
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CONDUCTIVITY  HEIGHT/WAVELENGTH 


.1 

.001 
,00001 


.5 

.3 
.1 
.5 

.3 

.1 
.5 
.3 
.1 


ASAP 

84.200    23.69 

132. 450    46.09 

39.460    84.99 

88.090    31.80 

117.690    42.3  9 

71.920    71.07 

90.370    31.33 

115. 770   45.06 

70.070    63.07 


EXACT* 

87.740  41.77 
136.00  67.91 
40.540  88.71 
91.410  50.65 
120.30  62.69 
78.770  69.91 
93.850  50.28 
117. 90  66.07 
77.670    60.81 
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SYSTEM  MANUAL 

INTRODUCTION:       The  Ant ennas-scatt erers  Analysis  Program 
(ASAP)  for  thin  wire  structures  in  a  homogenous  conducting 
medium  performs  a  frequency  domain  analysis  of  antennas  and 
scatterers.   The  program  is  applicable  in  the  presence  of  a 
ground  either  perfect  or  finite.  This  appendix  will  describe 
the  computer  program  which  accomplishes  this.   Although  the 
program  was  written  for  the  IBM  360  computer  system  it  can 
be  executed  on  another  system  with  minor  modifications. 

A  piecewise-sinusoidal  expansion  is  used  for  the  current 
distribution.   The  matrix  equation  ZI  =  V  is  generated  by 
enforcing  reaction  tests  with  a  set  of  sinusoidal  dipoles 
located  in  the  interior  region  of  the  wire.   Since  the  test 
dipoles  have  the  same  current  distribution  as  the  expansion 
modes,  this  may  be  regarded  as  an  application  of  Galerkin*s 
method.  Rumsey's  reaction  concept  was  most  helpful  in  this 
development,  and  therefore  the  formualtion  is  known  as  the 
"sinusoidal  reaction  technique". 

The  main  routine  and  each  subroutine  is  discussed 
separately  in  this  appendix.  The  writeups  for  the 
subroutines  are  arranged  alphabetically  by  subroutine  name 
after  the  main  program.  Each  of  the  discussions  includes  the 
purpose  of  the  subroutine,  brief  description,  and  a  listing. 
After  the  subroutine  writeups  is  a  table  of  the  more  common 
symbols  used  in  this  program. 

The  input  data  and  program  limits  are  discussed  in 
detail  in  the  next  appendix  titled  "USERS  MANUAL". 
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GROUND  EFFECTS:      In  the  modified  antenna  analysis 
computer  program  finite  and  infinite  ground  effects  were 
added  by  using  the  reflection  coefficient  technigue.   The 
method  in  which  this  technigue  was  used  reguired  the 
generation  of  an  image  structure.   In  this  section  the 
reflection  technigue  will  be  discussed  in  detail. 


In  order  to  apply  ground  effects  to  the  electric  field, 
the  field  for  the  image  structure  was  first  calculated  as  if 
a  ground  were  not  present.  Then,  the  field  was  decomposed 
into  parallel  and  perpendicular  components.  (A  parallel 
component  is  the  component  which  is  parallel  to  the  plane  of 
incidence.  A  perpendicular  component  is  one  which  is 
perpendicular  to  this  plane.  The  plane  of  incidence  is  the 
plane  containing  the  normal  to  the  reflecting  surface  and 
the  incident  ray.) 

Z  1 


■7-Y 


*-V 
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Consider  an  image  monopole  with  the  electric  field  in 
the  1  direction.   The  ray,  G,  is  a  vector  which  is 
perpendicular  to  1  and  passes  thru  the  point  of  interest.  To 
apply  reflection  technique,  the  plane  of  incident  must  be 
found.  It  is  advantageous  to  define  a  new  coordinate  system 
(a#/3#z)  where  a  and  p   are  parallel  to  the  xy  plane  with  a   in 
the  plane  of  incident  and  /3  perpendicular. 


If  the  direction  cosines  (cosx,  cosy,  and  cosz)  are 
known,  it  can  be  shown  that  the  components  of  the  field  in 
the  a/3  (xy)  plane  have  the  following  relationship: 


E 


I  i 


sin  $    -cos  <t> 
where  $   =  arctan (cosy/cosx)  . 


cos  <(>   sin  4> 


E 


E 


Now  the  reflection  coefficients  for  the  interface  can  be 
applied  as: 

11(H)  || 

E       =  R    E 
II 

±(R)        jL 

E       =  R    E 

where  R    and  R   will  be  defined  later  in  this  section. 
II       1 

Applying  the  matrix  equation  above  yeilds: 
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E 

X 

= 

*(R) 

E 

y 

cos  4>    sin  $ 


sin  $   -cos  4> 


(R) 


(R) 


(the  square  matrix  is  unique,  in  that,  the  inverse  is  equal 
to  the  oriqinal  matrix)  .  Since  the  image  direction  is 
opposite  to  the  original  monopole,  that  is, 

(1  x  z)         =  -  (T   x  z)       , 

original  image 

the  2  component  of  the  field,  which  is  in  the  plane  of 
incident,  is  given  by: 

(R) 
E     =  -  R    E   . 
z        I  I   z 

Prom  electro-magnetic  theory  the  reflection  coefficients 
for  the  fields  in  medium  (1)  at  the  interface  with  another 
medium  (2)  are  defined  as: 

for  perpendicular 

cos  0  -  /  6"  -  sin2  e 
H    cos  9  +  S  &*    -    sin2  0 

and  for  parallel 

&•  cos  9  -  v/  &•  -  sin2  e 

R   =  ===== 

V      &•  cos  e  +  v/  &•  -  sin2  e 

where  0  is  the  angle  of  incident  as  measured  from  the  normal 
to  the  interface  and 


&•  =  (&   ♦  a/j6j)/(&   +  a  /jco) 

2  2  1       1 


where  the  subscripts  correspond  to  the  mediums  above. 
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To  determine  the  relationship  between  R   ,  R   and  R  ,  R 

I  |    J_       V    H 

a  perfect  ground  (&      =  0,  a  =oc  )  was  investigated. 

r 

limit  R   *  -1 
H 

limit  R   =  +1 
V 

But,  for  a  perfect  ground  the  contributions  to  the  field 
from  the  image  monopole  would  be  equal  to  the  field  of  the 
original  monopole  but  opposite  in  sign  due  to  the  chosen 
reference  direction, 

II  (R)         ]  I       II 

E       =  R    E    =  -  E 
I  I 

Ei(S>   =  R    ^   -  -  E^ 
J_ 

therefore 

R    =  -  R 
II        V 

ti    =  R 
1      H 


In  summary,  the  contribution  to  the  electric  field  of  a 
monopole  over  a  ground  plane  at  a  given  point  is  given  by: 

(R)      (R)  (R)  (R) 

E    =  E     cosx  ♦  E     cosy  +  E     cosz 
x  y  z 

where 

(H) 

E     =  R   E  cosx  +  (R    -  R  )  E  cosx  cos2  $ 

x      _L  II    _L 

+  (R    -  R  )  E  cosy  sin  4>  cos  $ 
II    -L 
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E     =  R    E  cosy  -  (R    -  R  )  E  cosy  cos2  <J> 

y     II  1 1    _L 

+  (R    -  R  )  E  cosx  sin  <p   cos  $ 
I  I    X 

(R) 

E     =  -  R    E  cosz 

z         I  I 

where  E  is  the  field  without  the  ground  plane  present  and 

-  &•  cos  e  +  v/  &•  -  sin2  6 


II 


cos 


&•  cos  6  ♦  \/  &'  -  sin2  e 
e  -  v/s»  -  sin2  e 


X    cos  0  ♦  y/  &•  -  sin2  e 


&'  =  &   -  j  (a/&  o>) 
r        o 
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MAIN 


PURPOSE:       to  control  the  input,  output,  and  the  flow  of 
calculations. 

METHOD:         The  main  program  controls  the  flow  of  the 
required  calculations  by  calling  only  a  few  subroutines. 
These  subroutines  in  turn  call  other  subroutines  which 
actually  do  the  required  calculations.  The  order  of  the 
calling  sequence  is  diagramed  after  the  listing  for  the  main 
program. 

The  DIMENSION  statements  at  the  beginning  of  the  main 
routine  provides  the  required  storage  for  a  wire  structure 
with  up  to  50  segments,  60  nodes  and  60  dipoles  without  the 
presence  of  a  ground  plane.  If  a  ground  plane  is  present 
one-half  of  the  reserved  storage  is  required  for  the  image, 
therefore  a  wire  structure  with  up  to  25  segments  and  30 
nodes  can  be  analyzed. 

NM  denotes  the  actual  number  of  monopoles  (segments) , 
INM  is  the  corresponding  dimension,  and  the  dimension  for 
CG,  VG,  and  ZLD  is  twice  INM.  The  second  subscript  for  MD 
always  has  a  dimension  of  4  to  correspond  to  the  number  of 
segments  meeting  at  a  given  node. 

N  denotes  the  number  of  simultaneous  linear  equations 
and  ICJ  is  the  corresponding  dimension.  The  dimension  for  C 
is  (ICJ  *  ICJ  +  ICJ)/2. 

In  the  statements  above  statement  4,  the  initial 
conditions  and  defaults  are  established.  After  calling 
subroutine  READ  to  determine  the  input  parameters,  the  IF 
statements  output  the  parameters  to  be  used  for  the 
calculations.   In  the  DO  LOOP   ending  at  statement  7,  the 
the  input  data  of  the  structure  geometry  is  stored  in  order 
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to  recall  if  the  structure  is  to  be  moved  for  ground  plane 
calculations. 

After  the  image  structure  is  generated  and  structure 
location  is  moved,  subroutine  SORT  is  called  to  determine 
the  dipole  modes.   Prior  to  calling  SGANT,  the  load  and 
generator  information  is  established. 

Subroutine  SGANT  is  then  called  to  calculate  the 
elements  of  the  impedance  matrix.  If  FEEDS  or   GENERATORS 
are  specified  by  the  input  data  stream,  subroutine  GANT1  is 
called  to  solve  for  the  current  distribution  due  to  these 
forcing  functions. 

In  the  DO  LOOP  ending  with  statement  29,  subroutine 
GNFLD  is  called  to  calculate  the  near-zone  field  for  the 
current  distribution  of  the  subroutine  GANT1. 

The  subroutine  GFFLD  is  called  for  the  far-zone  field  of 
the  current  distribution  of  the  subroutine  GANT1  in  the  DO 
LOOP  ending  at  statement  35.  The  subroutine  GFFLD  is  called 
again  in  DO  LOOPs  ending  at  statements  42  and  51,  if 
bistatic  and  backscattering  calculations  are  requested  by 
the  input  data  stream. 


CALLS  TO: 


GANT1 

GFFLD 

GNFLD 

POLPRT 

READ 

SGANT 

SORT 
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DIMENSION  X(bO),  Y(60).  Z(60),  XGC60),  YG(60),  ZG(60I  0001 

OIMENSION  11(60),  12(65),  13(60),  JA(63>,  JB<60>,  KFLAG(30I  0002 

OIMENSION  CPHH500),  CTHET(500I  0003 

DIMENSION  0ATYH360),  DATY2I360),  DATY3O60J.  DATY*(36Q)  000* 

DIMENSION  0(50).  IA(50I.  18(50),  ISC<50),  M0(50,*>,  NO(SO)  0C05 

DIMENSION  LZD(60I,  KGEN(60»  0006 

COMMON  IWL  0007 

DIMENSION  XNP150I,  YNP(50),  ZNP(50)  0008 

CGHPLEX  CU830)  0009 

COMPLEX  COATl(500),CDAT2(500),COAT3(500lfCOAT*(500l  0010 

COMPLEX  CJ(60) ,EP(60) .EPP(60I,ET(60) ,ETT<60I  0011 

COMPLEX  CG0(5Qj,SG0(6&).CG( 105l,VG(100) ,ZLD(100I  0012 

COMPLEX  VOLT(60l,ZLlD(607  0013 

COMPLEX  EPPS,EPT$,ETPS,ETTS,EX,EYtEZ  001* 

COMPLEX  EP2,EP3,EP*,ERR.ETA,GAM,Yll,Zll,ZS  0015 

OATA  PI, TP/J. 1*159, 6. 28318/  0016 

DATA  EQ,U0/8.85*E-l2,l.2566E-6/  0017 

NGEN  -  -1  0018 

0019 
0020 
0021 
0022 
0023 
002* 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
003* 
0035 
0036 
0037 
0038 
0039 
00*0 

88*2 

00*3 
00** 
00*5 
00*6 
00*7 
00*8 


00*9 

0050 

0051 

0052 

0053 

005* 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

006* 

CALL    READ    ( IA.I8,ieiSCtICAR0,IGAIN,IGR0,INEAR,INT, ISCAT.  IMRflFLAG,  0065 

1KFLAG,KGEN,L0AD,LJ0,MSG,NBAP.NBIP, NF  FP, NGEN ,NM,NP, AB AP. ABAT. AFFP. A  0066 

2FFT.ABIP, ABIT.AH,BM,CMM,ER2,ER3iER*,FMC,HGT.PHAF. PHAI,PH1F.PHII,PH  0067 

3SF.PHSI ,THAF,THAl.THIF,THII . THSF , THSI , S I G2 , SIG3, $ IG*,T02 , T&3 , VOLT ,  0068 

*X,XNP,Y,YNP.Z,ZIl6,ZNP,STEPI  0069 

WRITE    (6,56)  0070 

IF    (MSG.LT.l)    SO    TO    5  0071 

IF    (MSG.EQ.l)     WRITE    (6.70)    KFLAGOO)  0072 

IF    (IFI.A&.EQ.*)    GO    TO    I  0073 

5    IF    (IFLAG.EQ.5)    STOP  007* 

IF    (AM.LT.O)    W*ITE    (6,1271  0075 

IF    (AM.LT.O)    GO    TO    6  0076 

IF    ((NM.GT.O).AND.(NP.GT.O) )    GO    TO    7  0077 

WRITE    (6,116)  0078 

S.    IF    (IFLAG.EQ.l)    GO    TO    1  0079 

MSG   ■    2  0080 

GO    TO    *  0081 

7    WRITE    (6,11*)  0082 

WRITE    (6,113)  0083 

WRITE    (6,112)  008* 

WRITE    (6,83)     FMC  0085 

WRITE  (6,8*1  AH  0086 

r  mri.Uiji,cw.w  WRITE  (6,85)  CMM  0087 

IF  (KFLAG(20I.NE.l)  WRITE  (6,87)  0088 

IF  (KFLAG(*).EQ.l)  WRITE  (6,66)  0089 

IF  (KFLAG(*).Ew.l)  WRITE  (6,881  BM  0090 

IF  (KFLAG(5).E0.1)  WRITE  (6,89)  SIG2  0091 

IF  (KrLAG(6).EQ.li  WRITE  (6,90)  ER2  0092 

IF  (KFLAG(7).EQ.l)  WRITE  (6,91)  TD2  0093 

IF    Jr 

IF    UFLAG<10f.eQ.i)    WRlfE    (6,9*)     ER3  0096 


IGRD  -  -I 

LOAD  -  -1 

BM  «  -1 

ICARD  »  0 

AH  «  -1 

I  FLAG  ■  0 

VOLT(l)  -  Jl.,0.1 

HGT  «  0. 

NH  «  0 

NP  -  0 

HSG  »  0 

SIG2  »  -I. 

TD2  «  -I. 

SIG3  -  -1 

ER3  «  I 

TD3  ■  0. 

CMH  ■  50. 

ER2  -  1. 
FMC  "  300. 

INM  «  50 

ICJ  -  60 

c 

WRITE  (6,7*1 

DO  2  1-1,30 
2  KFLAGIM  •  -1 

C 

c 

DO  3  J-l, INM 
ISCU)  ■  0 

VG(J)  -  (.0..0) 

ZLO(J)  «  (  .0,  .0) 

JJ  «  J*INH 

VGIJJt  -  (.0..0) 
3  ZLDUJI  -  (.6, .6) 

r 

L 

*  NFFP  -  0 

NBIP  ■  0 

NBAP  ■  0 

AFFP  -  1000. 

AFFT  -  1000. 

AB1P  ■  1000. 

ABIT  »  1000. 

A6AP  =  1000. 

A8AT  =  1000. 

STEP  -  1. 

KNM  =  0 

MKIIC       IOiIUI 

IF  (KFLAG< ll.EQ.l) 
IF  (KFLAG(2  .EO.l) 
IF    (KFLAG1 3) .EO.l) 


!'sIrI  88?? 
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If 

IF 


(KFLAG( lll.EU.l)    WRITE 

'  .NE    " 


6,95)    TD3 
AND.(KFLAG(  25I.EQ.1) 


1)    WRITE 
ANO.(KFLAG 


6,12< 

25) .EQ.l)) 


WRITE    (6,1231 

WRITE    (6,125) 

WRITE    (6,  124)    ER4.SIG4 

WRITE    (6,126)    HGT 


(KFLAGI26) 

(< I3RD.GT.il. 

<(  I&RD.EQ.l). 

((  IGRO.GT.l  j  .  ANO.  (KFLAG»  25)  .ECU 

((  IGRD.GT.OJ.ANO. (  KFLAG(  25) . EQ. 1) ) 

(KFLAGI  21). EQ.l  )    WRITE    (6,121)     INT 
WRITE     (6,111) 
IF    KFLAGI  12).  EQ 

it),  i: 

WRITE    ( 6. Ill  ) 

IF    (KFLAGI24J.GT.0)    WRITE    (6,119)    (LZO(  I ), ZLLO( I ) , t-l.LOAD) 

IF    (KFLAG(14).GT.0t    WRITE    (6,118)    (LZD( 1 ) , ZLLD( I ) , I- l.LOAD) 


HI. EQ.l)    WRITE    (6.96)     ( I  ,  IA(  I  ) ,  X  (  I  A(  I  )  )  ,  Y(  I  A(  1 1  )  ,  Z(  I  A<  I 
B(I>),V(IB(I)),Z(IB(l)),I«l,NM) 


,NE, 
.EQ. 


II 


WRITE     (6.111) 

IF    (KFLa6(23) 

IF    (KFLAGU3) 

WRITE     (6,111) 

WRITE    (6,114) 

WRITE    (6,98) 

WRITE    (6,112) 

IF      KFLAG(22) 

IF    (KFLAGU5) 

IF    (KFLAG( 16I.EQ.1 

IF    (KFLAGI17). EQ.lt 

IF    (KFLAGI  181.63.11 
~    (KFLAG( 19J.EQ.I) 
(AFFP.LT.500 
( AFFT.IT.500 
I  A8AP.lT.500 
[ABAT.LT. 500 

IF    IABIP.LT.500 

IF  (ABIT. LT. 500 


GT.O)  WRITE  (6,120)  (KGEN( I ) , VOLTI I ) , 
GT.O)  WRITE  (6,97)  ( KGEN(  I  I , VOLTI I  I , j 


I-l.NGEN) 

•l.NGEN) 


IF 

Vf 

IF 
IF 


WRITE  (6,110) 

WRITE  (6,99) 

WRITE  (6,100)  PHftl, PHAF, THAI ,THAF, STEP 

WRITE    (6,101)    PHIlfPHIF,THII,THlF,STEP 

WRITE    (6,102       PHSI.PHSF.THSI.THSF.STEP 

WRITE    (6,103)     UNPt  I  ),YNP<I),ZNP(iT,1-l,INEAR) 

WRITE    (6,165)    AFFP 

WRITE    (6,104)    AFFT 

WRITE    (6,109)     ABAP 

WRITE    (6,108)    ABAT 

WRITE    (6,107)     ABIP 
"     16,106)     ABIT 


WRITE 
IF    (( 1BISC. GT.O) .ANO. (ISCAt.LT.O))    WRITE    (6,731 
IF    (KFLAG(4).LT.l)    GO    TO    129 
00    128    I>1,INH 

128  ISCIII-I 

129  FHZ    •    FMC*1.E6 
OMEGA    «    TP*FHZ 

IF    (SIG2.LT. 0.)    EP2»ER2*E0*CMPLX(1.,-TD2) 
JF    (TD2.LT. 0.)    EP2    -    CMPLX(ER2»E0,-SIG2/OMEGA) 

F    (SIG3.LT.0.)     EP3»ER3*E0*CMPLX(l..-TD3T 

F    (T03.LT.0.)    EP3    -   CMPLXj ER3*§0,-$IG3/OM§GA J 


IGRD.GT. 


EP4    -    CHPLX(ER4*E0.-S1G4/0HEGA) 


IF    (IGRO.GT.ll    ERR    ■    EP4/EP3 

IF    (KFLAG(21).GT.0I    WRITE    (6,121)    INT 

ETA  -    CSQRT(U0/EP3I 

GAM  -    0MEGA*CSQRT1-U0*EP3) 

IF  (KFLAG(12).NE.l)  GO  TO  9 


0097 
0098 
009V 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
0110 
0111 
0112 
0113 
0114 
0115 
0116 
0117 
0118 
0119 
0120 
0121 
0122 

mi 

0  25 
0  26 

mi 

0129 
0130 
0131 
0132 
0133 
0134 
Q  35 
0136 
0137 
0138 
0  39 
0140 
0141 
0142 
0143 
0144 


NPG 
NMG 


NP 
NH 


00  8  I 
XG(I 
VG  " 
8  ZG 


III  :  !!ll 


9  DO  10  1-1. NPG 
X(l)  ■  XGi  I) 
Yd)  -  VG(  I) 
10  Z(I)  -  ZG(  I) 


NP 
NH 


NPG 
NHG 


IWL  *    0 
IF  (IGRO. 
SET  UP  IN 


ZHIN 
K  ■  0 


LE.O)  GO  TO  15 
AGE  FOR  GROUND  PLANE 
2(1) 


IF  (Z(  I  J.LT.ZMIN)  ZHIN-Zdl 

DO  11  I'l.NP 

Z(I)  *  7.(II*HGT 

IF  (Z( I I.GT.1.E-60I  GO  TO  11 

IWL  *  IWL* I 

11  CONTINUE 

IF  (2HIN.GE.O.0IGOTO12 

WRITE  (6,117) 

IF  (IFLAG.EQ.l) 

IF  (IFLAG.EQ.2) 

HSG  •  2 

GO  TO  4 

12  DO  13  J-l.NN 
K  •  J*NM 
IA(Ki  -  IA(J) 


GO  TO  1 
STOP 
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IF  <|A( JJ.GT. IWlJ  IA(K)-IA(J)*NP-IWL 


IWLP 


IB(J)*NP*IWL 
IWL*1 


DO  14  I-IWLP.NP 
J  ■  I*NP-|WL 
X(J)  -  X( 
Y(J)  - 


[ill 
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14  Z(J)  •  -Z(  II  019} 
C  0I9-* 

KNM  «  NM*1  0195 

NM  «  2*NN  0196 

NP  *  2-NP-IWL  0197 

15  CALL  SORT  < I  A, I 8 , 1 1 , 12, 13 , JA, JB, MD.NO, NM.NP iNi MAX ,HI N, IC Jt INMI  0196 
IF  (HAX.LE.4)  GO  TO  16  0199 
WRITE  (6.711  0200 
IF  (IFLAG.EQ.ll  GO  TO  1  0201 
IF  (IFLAG.EQ.2)  STOP  0202 
NSG  ■  2  0203 
GO  TO  4  0204 

16  IF  lt1IN.GE.ll  GO  TO  17  0205 
WRITE  (6,721  0206 
IF  (IFLAG.EO.ll  GO  TO  1  0207 
IF  UFLAG.Ed.2)  STOP  0208 
MSG  «  2  0209 
GO  TO  4  0  210 

17  WRITE  (6.56)  0211 
IF  (MAX.GT.4.0R.NIN.LT.1.0R.N.GT.ICJ)  GO  TO  34  0212 
J12  -  I  0213 
IF  (LOAO.GT.O)  GO  TO  19  0214 

C  0215 

DO  18  I-l.NM  0216 

16  ZLD(I)  -  (O. ,0.1  021? 

19  IF  (NGEN.GT.OJ  GO  TO  21  81 iS 
C  0220 

00  20  I-l.NH  0221 

20  VGII1  ■  10. ,0.)  0222 
C  0223 

21  KN  >  NM  0224 
IF  (IGRD.GT.OI  KN  -  NM/2  0225 
J  -  1  0226 

C      ANTENNA  CALCULATIONS  0227 

IF  (LOAD.LE.O)  GO  TO  24  0228 

IF  (KFLAGI24I.GT.0)  GO  TO  22                                            0229 

C  0230 

DO  23  J«l,KN  0?3l 

22  DO  23  1-1 .LOAD  0233 
K  -  LZD(I)  0234 

0235 
0236 
0237 

23  CONTINUE  0239 

0240 


K   -   LZD(t) 

IF    ((  IAUt.EO.Kl  .AND.(KFLAG(14I.GT.0I)    ZLD(  J) -ZLLDU  ) 

JF    JKFLAG(24).GT.0I    ZLD< K) -ZLLDI 1 1 

IF    {<KFLAGil4>.GT.0).AN0.{lGR0.GT.0)l    ZLOI J»KNI-ZLD( J| 

IF    <<KFLAG(24).GT.0).AN0.(lGRD.GT.0lJ    ZLD( K»KnJ-ZLD<KI 


24  IF  (NGEN.LT.01  SO  TO  27  0241 
KN  -  NM  0242 
IF  (IGRD.GT.O)  KN  •  NM/2  0243 
IF    (KFLAG(23).GT.O)    GO    TO    25  0244 

0245 
DO    26    J-1,KN  0246 

0247 

25  DO  26  I-l.NGEN  0248 
K    -    KGEM  1)  0249 


IF    ((  IA(  Jl.EQ.K)  .ANO.  (  KFLAGU  3  1  .GT.O)  )    VG(  J  1-VOLT I  1 1  0250 

IF    <( I5RD.GT.6) .ANO.(KFLAG( 23) .GT^ol)    VG< K»KN)  — VG(K 1 


IF    (KFLAG<23).GT.O)    VGI K) -VOLT ( 1 1  0251 

IF    ((KFLAGI13). GT.O). ANO. (IGRD.GT.0|)       VG( J *KN) — VG( J)  0252 


26  CONTINUE  0254 

0255 

27  CALL  SGANT  ( I A , I  6. INM , INT , I SC, I  1 . I  2, 13, JA. JB.MD. N. ND.NM.NP , AM ,BH,C  0256 
1 , CGO.CMM, 0, E d, £ fi j ETA,FHZ,GAM,SGD,X,r,Z ,ZL0,ZS, ERR, |GRD»  02  57 

IF  (N.GT.oI  GO  TO  28  0258 

IF  1IFLAG.EQ.2I  STOP  0259 

MSG  «  2  0260 

IF  (IFLAG.E3.il  GO  TO  1  0261 

GO  TO  4  0262 

28  IF  (NGEN.LE.O)  GO  TO  36  0263 
WRITE  (6,751  0264 
WRITE  (6,76)  0265 
WRITE  (6,77)  0266 
WRITE  (6,82)  0267 
CALL  GANtl  (IA, IB, INM, IWR, II, 12, 13,112, JA.JB.HD.N, NO, NM, AM,C, CJ.CG  0268 

l,CM*,0,EFF,GAM,GG.CG0,SG0,VG,yIltZll,ZL0,ZS,lGR6)  0269 

WRITE  (6,57)  EFF,GG,Zll  0270 

NEAR  FIELD  0271 

IF  (INEAR.LE.OI  GO  TO  30  0272 

WRITE  16,75)  0273 

WRITE  (6,78)  0274 

WRITE  (6,77)  0275 

0276 

DO  29  l-l. INEAR  0277 

XP  -  XNP(I)  0278 

rP  -  VNP( I  I  0279 

ZP  -  ZNPlU  0280 

CALL  GNFLD  t  I  A, I B, I NM , 1 1 . I  2 , 1 3. MD, N, NO, NM, AM, CGD, SCO ,E TA ,GAH,CJ ,D ,  0281 

lX,Y,Z,XP,rP.ZP,EX,EY.EZ,IGRO,ERRT  0282 

WRITE  (6,581  XP,rP, ZP  0283 

WRITE  (6,59)  EX.EY.EZ  0284 

29  CONTINUE  0285 


FAR    FIELO 


8IH 


30    IF    (IGA1N.LE.0)    GO    TO    36  0268 
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0289 
00  31  1-1.360  0290 

DATYK  II  -  0  0291 

0ATY2( ||  -  0  0292 

DATY3UI  -  0  0293 

31  0ATY4II I  »  0  0294 

0295 
WRITE  (6,751  0296 

WRITE  16,791  0297 

WRITE  (6,771  0298 

WRITE  (6,821  0299 

INC  *    0  0300 


NPL  «  -1  0301 

AG,<  16I.E3  11  W~: 
IF  (NFFP.EQ.ll  GO  TO  32  0303 


IF  (KFLAG( l6l.E3.ll  WRITE  (6,691  0302 


NPHA  »  (PHAF-PHAI l/STEP»l  0304 

NTHA  ■  (THAF-THAI J/STEP*!  0305 

GO  TO  34  0306 


32  IF  (AFFT.GT.500.1  GO  TO  33  0307 
NPL  «  I  0308 
NPHA  •  360  0309 
NTHA  «  1  0310 
PHAI  -  0.  0311 
THAI  »  AFFT  0312 
STEP  «  1.  0313 
GO  TO  34  0314 

33  NPL  «  2  0315 
NPHA  «  1  0316 
NTHA  »  360  0317 
PHAI  ■  AFFP  0318 
THAI  -  0.  0319 
STEP  -  1.  0320 

34  PH  «  PHAI-STEP  0321 
DO  35  K-l.NPHA  0322 
PH  «  PH»STEP  0323 
TH  «  THAI-STEP  0324 
00  35  I -1, NTHA  0325 

•  0.  0326 

PH5TH  «  0.  0327 


TH  »  TH*STEP  0328 

0329 

CGD.CC.CJ.I    ;iECSP,EC5T,EP;ET7EPP',El   EPPS; gPTS.ETPS.ETTS.GG,  0331 

2GPP,GTf .PH,SG0,SCSP,SCST,SPPM,SPTM,STPN,STTM,TH,X,Y,i,2LD,ZS,ETA,G  0332 

3AM, ERA,  IGRDI  0333 

ETMAG    •    CABSJETTSI  0334 

EPMAG  -    CABSIEPPSI  0335 

IF(ETMAG.GT.1.E-32I    PHSTH-57.295779*ATAN2( AIMAG(ETTS) .REAL(ETTS)  I  0336 


IF(EPMAG.GT.1.E-32I    PHSPH-57.295779*ATAN2« A IMAG( EPPS I , RE AL(EPPS) I  0337 

IF    (NPL.EQ.ll    DATYUK1-EPMAG  0338 

IF    (NPL.EO.ll    DATY2(KI«ETMAG  0339 

IF    (NPL.EQ.2I    DATYHl  l-EPHAG  0340 

IF    (NPL.EQ.21    0ATY2(I l-ETMAG  0341 

IF    JKFLAG(16l.NE.lJ    GO    TO   35  0342 

WRITE    (6,601    TH.PH.GTT, GPP, ETTS, ETMAG, PHSTH, EPPS, EPMAG, PHSPH  0343 

35  CONTINUE  0344 

0345 

WRITE  (6,561  0346 

IF  (NPL.LE.O)  GO  TO  36  0347 

CALL  POLPRT  (1.0ATY1I  0348 

CALL  POLPRT  (2.0ATY2I  0349 

BACK  SCATTERING^  0350 

36  IF  (ISCAT.LE.OI  GO  TO  54  0351 
WRITE  (6,751  0352 
WRITE  (6,801  0353 
WRITE  i6,77>  0354 
WRITE  (6,821  0355 
L  ■  0  0356 
NPL  ■  -1  0357 
INC  -  1  0358 
IF  (NBAP.EO.ll  GO  TO  37  0359 

0360 
0361 


I! 


NPHI  ■  (PH1F-PHI I»/STEP»1 
NTHI  »  (THIF-THin/STEP*! 
IF    (IWR.LE.O)    WRITE    (6,621 


IF    (IWR.LE.OI    WRITE    (6,621  0362 

GO   TO    39  0363 


37  IF  (ABAT.GT. 500.1  GO  T 

0  38 

NPL  -  1 

NPHI  -  360 

NTHI  -  1 
PHII  -  0. 

THII  -  ABAT 

STEP  «  1. 

GO  TO  39 

38  NPL  •  2 

NPHI  -  1 

NTHI  -  360 
PHII  ■  ABAP 

THII  -  0. 

39  PH  ■  PHli-STEP 

00  42  K-l.NPHI 
PH  -  PHfSTEP 
TH  •  THlT-STEP 

00  42  I -1, NTHI 


0364 
0365 
0366 
0367 
0368 
0369 
0370 
0371 
0372 
0373 
0374 
0375 
0376 


28 


IF    (UGRO.GT.0).ANO.((TH.GT.9O)  .AND. ( TH .LT .270) ) J    GO   TO    42 

CALLLGFFLO  ( IA,IB,INC.INM,IWR,I1,I2,I3,I12,MD,N,N^,NM,AM,AC$P,ACST 
l.CCGD.CG.CJ.CMM.O.ECSP.EdST.ep.eT.tPP.ETT.EPPS.EPTS.ETPS.ETTl.GG. 
>GPP,GTTtPH,SGD,SCSp\scST,SPPN,SPTN,STPH,STTM,TH,X,V,Z,ZL0,2S,ETA,G 


40 


3AM, ERR, IGRD) 
IF  (IWR.GT.O)  GO  TO 
IF  (NPL.LT.O)  rfRiTE 
IP,ECST, SCSP.SCST 
PH 
TH 


63)  PH,TH,SPPN,SPTM,STPK,$TTM,ACSP,ACST,ECS 


41 


41 


42 


43 


45 


cphkl) 

CTHET(L) 

CDAT1ILI  » 

CDAT2(L)  - 

C0AT3UI  ■ 

CDATMLI  » 
IF  (NPL.NE 

DATYKKI  • 

DATY2(K)  » 

DATY3IK)  « 

DATY4(K)  » 
GO  TO  42 

OATYl(I)  =• 

DATY2U)  ■ 

0ATY3I  I  J  - 

DATY4U)  » 
CONTINUE 


WRITE  (6,82) 

IF  (NPL.LE.O)  GO  TO  43 

CALL  POLPRT  17,DATY1) 

CALL  POLPRT  (8.0ATY2) 

CALL  POLPRT  (9.DATY3) 

CALL  POLPRT  (16.0ATY4) 

IF  (KFLAG(17).n£.1)  GO  TO  45 

WRITE  to, 641 


EPPS 
EPTS 
ETPS 
ETTS 
.1)  GO  TO 
CABSI EPP 
CABSI EPTS) 
CABSIETPS) 
CABS(ETTS) 


CABS(EPPS) 
CABS( EPTS) 
CABSt  ETPS) 
CABS(ETTS) 


00    44    l«l,L 

WRITE    (6,65)    CP-iI(  I ) ,CTHET ( I > ,COAT H I ) , CDAT2I I ) ,C0AT3< I > , CDAT4J  I) 


8ISTATIC    SCATTERING 

IF  (IBISC.LE.O)  GO  TO  54 

WRITE  (6,75) 


WRITE  (6,81. 
WRITE  (6,771 
WRITE 
WRITE 
WRITE 


6,771 

46, 82) 


CPHKL), CTHET(L) 


0385 

0386 

0387 

0398 

0389 

0390 

0391 

0392 

0393 

0394 

0395 

0396 

0397 

0398 

0399 

0400 

0401 

0402 

0403 

0404 

0405 

0406 

0407 

0408 

0409 

0410 

0411 

0412 

0413 

0414 

0415 

0416 

0417 

0418 

0419 

0420 

0421 

0422 

0423 

0424 

0425 

0426 

5*21 

0428 

0429 

0430 

0431 

0432 


46 


L  «  0 

INC  -  2 

NPL  »  -1 

IF  (NBIP.EO.l 1  GO  TO  46 

NPHS  -  (PHSF-PHSI)/STEP*1 

NTHS  -  (THSF-IHSH/SfEP*! 

IF  (IWR.LE.O)  WRITE  (6,671 

GO  TO  48 

IF  (ABIT. GT. 500. I  GO  TO  47 


NPL 

NPHS 

NTHS 

PHS1 

THSI 

STEP 


1 


GO  TO  48 


360 

1 

0. 

ABIT 

1. 


47  NPL 
NPHS 
NTHS 
PHSI 
THSI 
>TEP 


48 


SI 
PH 


00  51 
PH 


TH 
IF 

?8 


-  i 

-  360 

■  A6IP 

-  0. 

■  1. 
PHSI-STEP 

K«l .NPHS 
PH«-STEP 


«  THSl-STEP 

(UGRO. GT.O). AND. ((TH.GT. 90). AND.UH.LT. 270)1)  GOTO  51 


51 


1-l.NTHS 
TH»STEP 


L»l 


CALL  GFFLD  ( I A , I B , I NC , INN, IWR . I  1,12,13 , 1 12 , MD.N.NO.NH. AM. ACS? , ACST 
l.C.CGD.CG.CJ.CMM.O.ECSP.EdST.EP.ET ,EPP,ETT,EPPS,  EPTS, ETPS, ETU.GG, 
2GPP,CTf .PH.SGO, SCSP.SCST, SPPH, SPTM.ST PM , ST TM.TH.X.Y, 2,ZLD,ZS, ETA, d 
3AM,  ERR, IGRD) 

IF  ( IWR.GT.O)  GO  TO  49 

IF    (NPL.LT.O)    WUTE    (6,63)    PH,TH , SPPM, SPTH, STPM, STTH 


49  CPHKL) 
CTHEHLi 
COATUL 
T2(L 


PH 
■    TH 


COAT 

CDAT3(L 

C0AT4U 

IF    (NPL.NE.l) 

DATY1CK 

DATV2 

0ATY3 


EPP| 


EP1 

ETPS 

ETTS 

GO 


(K)  -  CA8S(EP 
!(Kt  -  CABSJEP 
I'M    -    CABS1ET 


TO  50 

PPSI 

TS 

PSl 


0433 
0434 
0435 
0436 
0437 
0438 
0439 
0440 
0441 
0442 
0443 
0444 
0445 
0446 
0447 
0448 
0449 
0450 
0451 
0452 
0453 
0454 
0455 
0456 
0457 
0458 
0459 
0460 
0461 
0462 
0463 
0464 
0465 
0466 
0467 
0468 
0469 
0470 
0471 
0472 
473 
474 
475 
0476 
0477 
0478 

0480 


29 


50 


51 


52 


DATY41KJ-    CAB51EUS) 


IF  (NPL.NE.2)  GO  TO 
OATYHI)  ■  CABS(EPPS) 
CABS(EPTS) 
CABS(ETPS) 
CABS(ETTS) 


0ATY2UI 
0ATY3U  ) 
DATY4( I) 
CONTINUE 


WRITE    (6,82) 

IF    (NPL.tE.OI    GO    TO    52 

CALL    POLPRT    (3.0ATY1) 

CALL    POLPRT    (4.0ATY2 

CALL    POLPRT    (5.DATV3) 

CALL    POLPRT    (6.DATY4) 

IF    (KFLAGI18I.NE.1)    GO    TO    54 

WRITE    (6,66) 


53 

5* 


DO    53    J«l,L 
WRITE    (6,65) 


:PH  I  (  I  >  ,C  THET  I  I )  ,  C0AT1  (  I ) ,  COAT  2(  1 1  .  C0AT3  (  I  ) ,  C0AT4(  1 1 


5S 


IF    (IFLAG.EQ.l)    GO    TO    1 
IF    (IFLAG.E0.2I    STOP 

KKFLAG-0 

KJFLAG=0 

KMFLAG=0 

KNFLAG'O 

IF    (KFLAGI 13). GT.O)    KKFLAG-1 

IF    (KFLAG(23).GT.0J    KJFLAG-1 

IF    (KFLAG(  14J.GT.0)    KMFLAG-1 

IF    (KFLAG(24).GT.0I    KNFLAGM 

DO    55    1-1,30 

KFLAGI I)    -    -1 

KFLAGI8I    -    1 
KFLAG(20)    «    I 
KFLAGI26I    -    1 
IF    (KKFLAG.GT.O) 


If 


(KJFLAG.GT.O) 
(KMFLAG.GT.O) 
^KNFLAG.GT.O) 


KFLAGI 13 
KFLAGJ23 
KFLAGI14 
KFLAGI24 
WRITE 


»1 

: 

68} 
115) 


UFLAG.EQ.3l    WRITE      6 
IF    IIFLAG.E0.6)    WRITE    (6, 
GO  TO   4 

FORMAT  (1H0) 

FORMAT  (10X,* THE  RADIATION  EFFICIENCY  IS 


•,F15.7//10X,'THE    TIHE-AV 


0481 
0482 
0483 
0484 
0485 
0486 
048T 
0488 
0489 
0490 
0491 
0492 
0493 
0494 
0495 
0496 
0497 
0498 
0499 
0500 
0501 
0502 
0503 
0504 
0505 
0506 
0507 
0508 
0509 
0510 
0511 
0512 
0513 
0514 
0515 
0516 
0517 
0518 


runnai     nuAt'inc    kmuimiiupi    crniicraiT    is    '  i  r  IP.  if  /  iua,  ■  i  nc    iinc-av  Witt 

1ERAGE    POWER    INPUT    IS    ' ,F15. 11) 10X,  'THE    ANTENNA    IMPEOANCE    IS    «,F15.  0526 


27.'    ♦J',F15.7//I 
58   FORMAT    (lOX,'TH£    NEAR-FIELO    ELECTRIC _ FI ELD    INTENS I TY _AT_THE    OBSERV 


1ATI0N    POINT    '.EU.St'.'.Ell.Si'i'iEll.Si'    IX.Y.Z    RESPECTIVELY)    IS: 
2'//) 

59  FORMAT    (  20X,  •  EX»= '  ,F15  .  7,  •     ♦  J«  ,  F  1  5  .  7/20X  ,  •  EY-'  ,  Fl  5.  7,  •    ♦J',F15.7/20 
IX, 'EZ*' .F15.7,"    ♦J,.F15.7////i 

60  FORMAT  (3X,F5.1,2X,F5.1,3X,E10.4,2X,E10.4,2(3X,3(E10.4,2X),F6.1,1X 
1)) 

61  FORMAT  <T41,'F0R  BISTATIC  SCATTERING  THE  INCIOENT • /T41 , 'PLANE  WAVE 


62 


GMA' ,T66,' ABSORPTION ',T90,' EX 
"  I  INCIDENT- SCAT TEREO) ',  " 


runnni     \  i  hi  #  -  r  u^     di^iaiiu     3d  ick  in\j 
1     IS    PHI«'  ,F5.1.«     THETA«',F5.1///» 

FORMAT    (•     INCIDENT', T27, 'ECHO    AREA    SI 
1TINCTI0N',TU4,' SCATTERING'/'       PLANE* , T2 5, '  I  INC  I DENT-SCATTEREO) ' , 1 
24X, 31 5X, 'CROSS    SECTION' ,6X) /•       WAVE    • ,52X ,3  1 10X , • FOR' , 1 1  X) /•       PHI 
3    THET*'  ,3X,'  PHI- PH I ' , 3X , • PHI-THETA' ,4X, 'THETA-PHI' , 2X, • THE TA- THET A 
4'  ,31  5X, 'PHI' ,7X,' THETA' ,4X)  ) 

63  FORMAT     I1X.2JF5.1, 1X1 , l6l E 10.4, 2X  ) ) 

64  FORMAT  ( T54, • BACKSCATTER ING' /•  I NC IDENT • ,T37, • EL ECTR IC  FIELD  POLAR 
1IZAT10N  SCATTERING  MATRIX'/'  PL ANE ' . T49 , • ( INCI OENT-SCATTERED ) • /3X 
2, 'WAVE'  ,T 23, 'PHI -PHI'  , T49,' PHI -THE TA'.T 7 5, 'THETA-PHI ',T 102, 'THET A- 
3THETA'/'       PHI       THETA' ,3X, 41 3X, 'REAL' ,8X, • lMAG',8X) ) 

65  FORMAT    I1X,2(F5.  1,  IX) . 2X , 41 E 1 1 . 5,2X. E 1 1 .5,3X1) 

66  FORMAT    IT54,'BISTATIC'/T37, 'ELECTRIC    F I E LO    POLARIZATION    SCATTERING 

1  MATRIX'/'    OBSERVATION', T56,' I  I NC IDENT- SCAT TEREO) • /•  POINT'. 14X, 

2  'FHI-PHI' ,T49,' PH I- THE t A ',T76.' THETA-PHI' ,  TlOl , • THETA-THETA' /•       P 
3HI       THETA' ,4X,4(3X .'REAL' ,8X,' tMAG'.SX) ) 

67  FORMAT    I'    6BERSV AT  I ON' ,T27 , • ECHO    AREA    SIGMA'/'  POINT « ,T25. « I INCI 

1  DENT-SCAT TEREO) •/•       PHI       THETA' . T 14, • PH I -PHI • , T24, 'PHI -THETA' ,T37, 

2  'THETA-PHI' ,T48, 'THETA-THETA' i 

68  FORMAT  1 1H1 , 5X ,  •  CfiNT INUE  EXECUTION  WITH  THE  FOLLOWING  ADDITIONS  AN 
ID/OR    CHANGES'//) 

69  FORMAT  ( 54X, • ELECTRIC  FIELD  I NTE NSt TV' / 5X, • DEGREE S« ,11X, • POWER  GAI 
IN'  ,28X,'THETA',42X,'PHI'/3X,'THETA',3X,'PHI '.7X, 'THETA' ,8X,'PHI  •  , 1 
2X.2(8X,'REAL' ,8X. 'J MAG' ,8X, • NAGN' ,5X, • PHASE ' ) ) 

70  FORMAT    ( 10X, • ♦•♦♦•ERROR    IN    DATA    CARD   NJMBER    ',12,'       EXECUTION    STOP 


1PED******* • ) 

71  FORMAT  (40X,'» 
1UR     ♦•/40X,«» 
2R0      ♦•/40X,'^ 
3         ♦•) 

72  FORMAT  (40X, •♦ 
IS      ♦«/40X,'^ 
2R0      •»/40X,'^ 
3         •') 

73  FORMAT  OOX.'A  BACKSC ATTERING  CALL  MUST  BE  INCLUDED  FOR  A  BISTATIC 
1  CALLV/50X, 'REQUEST  IGNOREO'///// ) 

74  FORMAT  I • 1' , T50. 371 • •• 1/T50, •♦• , T86, ••• / 

1  T50,«^       OHIO  STATE  UNIVERSITY         ♦•/ 

2  T50,'^     ANTENNA  ANALYSIS  PROGRAM        »•/ 


A  WIRE  SEGMENT  MAVNOT  BE  SHARED  BY  MORE  THAN  FO 

DIPOLE  MODES CHECK  DESCRIPTION  OATA  CA 

EXECUTION  STOPPED 

AN  ISOLATED  WIRE  MUST  HAVE  AT  LEAST  TWO  SEGMENT 

AND  THREE  POINTS CHECK  DESCRIPTION  OATA  CA 

EXECUTION  STOPPEO 


0529 
0530 
0531 
0532 
0533 
0534 
0535 
0536 
0537 
0538 
0539 
0540 
0541 
0542 
0543 
0544 
0545 
0546 
0547 
0548 
0549 
0550 
0551 
0552 
0553 
0554 
0555 
0556 
0557 
0558 
0559 
0560 
0561 
0562 
0563 
0564 
0565 
0566 
0567 
0568 
0569 
0570 
0571 
0572 
0573 
0574 
0575 
0576 


30 


3  T50,'« 

4  T50,'* 

5  T50,'* 

6  T50,'- 

75  FORMAT 

76  FORMAT 

77  FORMAT 
1*'  )) 

73  FORMAT 

79  FORMAT 

80  FORMAT 

81  FORMAT 

82  FORMAT 

83  FORMAT 
8*  FORMAT 

85  FORMAT 

86  FORMAT 

87  FORMAT 

88  FORMAT 

89  FORMAT 

90  FORMAT 

91  FORMAT 

92  FORMAT 

93  FORMAT 

94  FORMAT 

1  Ell. 5) 

95  FORMAT 

96  FORMAT 
It 18XI/T 
22(12, 5X 

97  FORMAT 

1  'REAL' 

98  FORMAT 

99  FORMAT 

100  FORMAT 

1  'AND  T 
2T50.UN 

101  FORMAT 
lit'  AND 
2T50.MN 

102  FORMAT 
1' ,F5.1, 
2T50,'IN 

103  FORMAT 
11.5,5X1 

104  FORMAT 

105  FORMAT 

106  FORMAT 


MOOIFIEO  FOR  USE  AT  ••/ 

NAVAL  POSTGRAOUATE  SCHOOL  *•/ 

3  SEPTEMBER  1975  ••/ 

T86,'*'/T50,37('*')l 


(«1ST50,29('^'1/T50,'»ST78,'*'» 

(T50,'»«  , 1 IX, 'ANTENNA', T 78,'* •  ! 

(T50,'*«  ,8X,' CALCULATIONS ',T78t'«'/T50,'^«  ,T78, • *• /T50 ,29( • 


(T50, 
(T50, 
(T50. 
(T50, 
(//// 
iT30, 
(T30, 
(T30, 
(T30, 
(T30, 
(T30. 
<T30, 
(T30, 
U30, 
<T30, 
(T30, 
(T30, 


'♦' ,9X 
•♦' ,9X 
•♦' ,7X 
••• ,4X 

) 

•FREQU 

•WHE 

•WIRE 

•WIRE 

'UIU 

• INSUL 

' INSUL 

•  INSUL 

■  INSUL 

•EXTER 

'EXTER 

•EXTER 


,'NEAR  FIELD' ,T78,'*'» 

,'FAR  FIELD' ,f78,'»' ) 

,' BACK  SCATTERING'  ,T7B,'*'I 

,'BISTATIC  SCATTERING'  ,T78,'»'  I 

ENCY  (MHZ! '.TBI. Ell. 5) 

RADIUS  (METERS! '.TBI, Ell. 5! 

CONDUCTIVITY  (MEGAMH0S/METERJST81.E11.5) 

INSULATED    «N0/VESIST85,'YES'! 

INSULATED    < NO/YES! ',T85, ' NO'I 

ATION    RAOIUS    <  METERS!  ST8  1.E1 1.  5) 

ATION   CONDUCTIVITY    (MHOS/METER  I S T81 , E 1 1 . 5) 

ATION   DIELECTRIC    CONSTANT     (REL AT  I VE! ', T8 1 , E 11 .5 1 

ATION   LOSS    TANGENT'. TBltEU. 51 

IOR    MEDIUM' .TBI, 'FREE     SPACE'! 

IOR    MEDIUM   CONDUCTIVITY    ( MHOS/METER! '.TBI . Ell . 5  I 

IOR    MEDIUM   DIELECTRIC    CONSTANT    ( RELAT I VE 1 • , TBI , 


IT30, 'EXTERIOR  MEOIUM  LOSS  TANGE 
(T50.'W1RE  STRUCTURE'//T20,'SEG' 
21,'NO.' ,3X,2(  •  NO. • ,9X.'X' ,13X, 
,E11.5,4X,Ell.5,4X,Ell.St  1X1) I 
<T50, 'ANTENNA  FEEDS' /T40 , 'NODE S 
,  7X,'  IMAGINARY' /<T4ltI2.6Xt 21 4X, 
(T50,'*S  6X, 'OUTPUT  REQUESTED', 
(T30, 'STRUCTURE  CURRENTS'! 
(T30,'FAR  FIELDS  FOR  PHI  VARYING 
HETA    VARYING   FROM    SFS.lf(    TO    S 

STEPS    OF    '.F5.lt'    OEGREES.'l 
<T30,'BACKSCATTERING   FOR    PHI    VAR 

THETA    VARYING    FROM    '.F5.lt'    TO 

STEPS    OF    SF5.1.'    DEGREES. '! 
(T30. 'BISTATIC    SCATTERING    FOR    P4 
•     AND    THETA    VARYING    FROM    SFS.lt 

STEPS    OF    '.F5.li'    DEGREES. •» 
<T30.'NEAR    FIELDS    FOR    FOLLOWING 

T30,'PL0T  FOR  FAR  FIELD  THETA-' 
T30,'PL0T  FOR  FAR  FIELD  PHl-'.F 
T30,'PL0T    FOR    BISTATIC    SCATTERI 


NT'  , TBI, Ell. 51 

,4X,2< 'NODE ',19X, 'LOCATION' 

'  V'  ,  13X,'ZS  7X1/1  T2 1,1  2,5X, 

I6X, 'V0LTS'/T4l,'N0.«  tl2X, 
Ell. 5)11 

T78,'*') 

FR0MSUtF5.1,'  T0»,F5.1t 

YING  FROM  SF5.1.'  TO  SF5. 
SF5.1/ 

I  VARYING  FROM  SF5.lt*  TO 
•  TO  SF5.U 

POINTS  (X,Y,Z!'/50(T40,3<El 


.F5.11 

5.11 

NG-FOR  THETA-SF5.1) 


0577 

0578 

05  79 

0580 

0581 

058i 

0583 

0584 

0585 

0586 

0587 

0588 

05b9 

0590 

0591 

059? 

0593 

0591 

0595 

0596 

0597 

0598 

0599 

0600 

0601 

0602 

0603 

0604 

0605 

0606 

0607 

0608 

0 1 1  ■ 

0610 

0611 

0612 

0613 

0614 

0615 

0616 

0617 

0618 

0620 
0621 
0622 
0623 
0624 


i    FOR    PHI-SF5.il 

rA-',F5.1) 

*• .F5.ll 

:d'I 


107  FORMAT  (T30,'PL0T  FOR  BISTATIC  SCATTERIN 

108  FORMAT  (T30,'PL0T  FOR  BACK SCATTE R ING  THE 

109  FORMAT  (T30,'PLOT  FO*  BACK  SCAT TER I NG  PHI 

110  FORMAT  (T30,'N0  OUTPUT  OR  PLOTS  REOUESTE 

111  FORMAT  (//) 

112  FORMAT  1T50,'*'  , T 78 , • ♦ • /T50  .29 ( • *•  >) 

113  FORMAT  (T50,'»S  8X, 'INPUT  DATA  ST78,'* 

114  FORMAT  (T50,29('*« J/T50,'*' ,T7tt,'*'l 

115  FORMAT  (10X, 'SINCE  THIS  DATA  BLOCK  DOES  NOT  HAVE  A  TERMINATION  CAR 
ID  A  CHANGE  CARJ  IS  ASSUMED' ! 

116  FORMAT  (//10X,40( ••• I/IOX.'THE  DESCRIPTI 
1  STRUCTURE'/IOX.'MUST  8E  STATEO  IN  THE  F 
2*  EXECUTION  STOPPED  •♦*•> 

117  FORMAT  (//lOX.'NO  PART  OF  THE  WIRE  STRUC 
I  UNO  PLANE.  VIOX,  •♦♦♦♦EXECUTION  STOPPED^ 

118  FORMAT  (T50,« STRUCTURE  LOADS' /T40  ,' NODE ' 


I 


ON  ANO  THE  GEOMETRY  OF  THE 
IRST  DATA  BLOCK. «/lOX,'**» 


1  t'REAL't 7X,'  IMAGINARY'/(  T41  ,  I  2  •  6X  ,  2(  4X  , 

119  FORMAT  (T50, • STRUCTURE  LOADS' / T 39 , «SEGME 
lX,,fiEAL,,7X,,IMAGINARY'/(T4l,I2,6X,2(4X, 

120  FORMAT  (T50,' ANTENNA  FEEDS' /T39 ,' SEGMENT 
IX.  'REAL  S7X,  •  IMAGINARY'  /(T4  1,12,  6X  ,  2 I4X, 

121  FORMAT  (//TSOt'THE  NUMBER  OF  INTERVALS  F 
1S'/T30,'IN  THE  IMPEDANCE  MATRIX  WITH  SIM 
2i/T30,I3,S   IF  CLOSED  FORM  INTEGRATION 

122  FORMAT  (T30,'GR0UN0  PLANE  (  NO/YE  S!  ST85  , 

123  FORMAT  JT30,'GRQUN0  PLANE  (NO/YES) '.T85 , 

124  FORMAT  (T30, 'GROUND  OIELECTRIC  CONSTANT 

1  T30, 'GROUND  CONDUCTIVITY  (MHOS/METER)', 

125  FORMAT  (T30, 'GROUND  PLANE • ,T83. • PERFECT • 

126  FORMAT  (T30. 'ANTENNA  HEIGHT  (METERS! ',T8 

127  FORMAT  < // lOX ,40? ' •• J / 10X ,  •  THE  WIRE  RADI 
l'*'ll 

ENO 


TURE    CAN    LIE    BELOW    THE    GRO 

•  •••  ) 

,16X,'0HMS'/T41,'N0.S12X 

Ell. 5)1) 

NTS14X,  '0HMS'/T41,'N0S12 

Ell.  51!) 

S14X.'V0LTS'/T41,'N0.',12 

Ell  .5)1) 

OR    CALCULATING    THE    ELEMENT 

PSONS-RULE    INTEGRATION    IS' 

IS    REQUIRED    SET    INT-O'///! 

•NO*  I 

•YES') 

(RELATIVE)ST81tE11.5/ 

TBI, Ell. 5) 

) 

1.EU.5! 

US  MUST  BE  STATED'/10X,40( 


0625 
0626 
0627 
0628 
0629 
0630 
0631 
0632 
0633 
0634 
0635 
0636 
0637 
0638 
0639 
0640 
0641 
0642 
0643 
0644 
0645 
0646 
0647 
0648 
0649 
0650 
0651 
0652 
0653 
06S4 
0655 
0656 
0657 
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BLNK 


PURPOSE:       to  compress  data  to  the  left  by  removal  of 
the  blank  spaces  on  the  input  data  cards. 

METHOD:         A (I)  character  is  compared  to  the  blank;  and 
if  it  is  true,  the  A  (1+1)  character  is  shifted  to  the  A  (I) 
position. 


CALLED  BY: 


READ 


CALLS  TO: 


NONE 


SUBROUT INE  BLNK  (A) 
DIMENSION  4(80) 
DATA  BLANK/'  */ 
K  »  0 

00  1  I  -  I  .80 
J  «  l-K 
A(J|  «  All) 
1  IF  <A<  I  I. EQ. BLANK)  K-K»l 

IF    (K.EQ.O)     RETURN 
A(81-K)     -    BLANK 
RETURN 
END 


1 

2 

3 

* 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1* 
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CBES 


PURPOSE:        to  calculate  the  quantity  B01  where 

B01  =  J  (z)  /J  (z)  . 
o         1 


METHOD:         If  the  absolute  value  of  the  argument  for  the 
Bessel  function  is  less  than  12,  B01  is  calculated  via  the 
power  series  expansion  for  the  Bessel  function  in  the  DO 
LOOP  ending  at  statement  3.  If  greater  than  12,  the 
asymptotic  expression  is  utilized  at  statement  4.  If  the 
magnitude  of  the  complex  part  of  the  argument  for  the  Bessel 
function  is  greater  than  20,  B01  is  set  to  (0.,-1) .   If  the 
complex  part  of  the  argument  is  negative,  the  sign  of  E01  is 
changed  prior  to  returning  to  the  calling  program. 

CALLED  BY:      SGA NT 
CALLS  TO:       NONE 


SUBROUT INE  C8ES  ( Z.801  )  1 

COMPLEX  ARG.CC.CS.EX  2 

COMPLEX  B01 ,1, HRMJ,TERMN,MZ2«f JN(2J  3 

OATA  PI /3. K159/  * 

IF  (CABS1Z I .GE.12.0)  GO  TO  4  5 

FACTOR  *    0.0  6 

TfckMN  =  <0.,0.)  7 

MZ24  *  -0.25*Z*Z  8 

TERMJ  *  (  1.0,0.01  9 

10 

00  3  NP=1,2  11 

N  =  NP-  1  12 

JN(NP)  '    TERMJ  13 

M  «  0  1* 

M  =  MM  15 

TERMJ  «=  TERMJ»HZ2«/FL0AT(M*(N»M) )  16 

JN(NP)  =  JN(NP) »TERMJ  17 

IF  (NP.NE.l)  &0  TO  2  18 

FACTOR  =  FACT0«» 1 .0/FlOATIM)  19 

TERMN  =  TERMN»TERMJ*FACT0R  20 

ERROR  =  CA6S(TFRMJI  21 

IF  (ERROR. GT. 1.0E-10I  GO  TO  I  22 

TERMJ  =  0.5*1  23 

801  *  JN(1)/JN(2)  25 

RETURN  26 

Y  «  AIMAG(Z)  27 

IF  (ABS( Y)  .GT.20.1  GO  TO  5  28 

ARG  «  ( .0, 1.1*1  29 

EX  «  CEXP(ARG)  30 

CC  *    EX*1  ./EX  31 

CS  -  (.0,-l.)*(EX-l./EX)  32 

B01  *  ICS»CC)/(CS-CC)  33 

RETURN  3J 

B01  •  (.0,-1.)  35 

IF  (Y.LT.6.)  B01  -  (.0,1.)  3* 

RETURN  37 

ENO  38 
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DSHELL 

PURPOSE:        to  calculate  the  mutual  impedance  term 
contributed  by  the  dielectric  insulation  on  the  surface  of  a 
thin  wire. 


METHOD:         The  contribution  to  the  impedance  matrix  is 
calculated  utilizing  the  equation  below 


mn 


(S   -  &)  ln(b/a) 

2 

2n  jw&& 


/ 

m,  n 


F«  (1)  F»  (1)  dl  , 
m      n 


where  z    is  defined  in  subroutine  SGANT,  &      is  the 
mn  2 

dielectric  constant  of  the  insulation,  b  is  the  outer  radius 

of  the  insulation,  a  is  the  inner  radius,  &  is  dielectric 

constant  of  the  external  medium,  and  F  is  the  sinusoidal 

expansion  function. 


CALLED  BY, 


SGANT 


CALLS  TO 


NONE 


SUBROUTINE    DSHELL    ( AM, BM, OK .CGDS , SGOS , EP2 . EP. ET A ,GAM,P I 1.P12 ) 

COMPLEX    CGCS,SGDS,EP2  .EP.ETA,  GAM.PU.PIZ  tGO,CST 

OAIA    PI /3. 1*159/ 

GO    -    GAM'OK 

CST    •    I EP2-EP)»ETA»AL0G(BM/AHI/I4.»Pl»EP2»SGDS*SGDS) 

Pll    -    -CSr»IGD»SGOS*CGOSI 

P12   ■    CST»(GO*CGOS»SGDS) 

RETURN 

END 
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EQUAL 


PURPOSE:       to  determine  position  (location)  of  the  equal 
symbol  on  input  data  card. 

METHOD:        The  character  search  begins  in  the  column 

passed  to  the  subroutine.   On  returning  to  the  calling 

program,  the  argument  passed  is  the  column  following  the 
equal  symbol. 


CALLED  BY: 


READ 


CALLS  TO 


NONE 


SUBROUTINE  EQUAL  (Nl 
INTEGER  A,  EQULS 
COMMON  /A/  A(80l 
DATA  EQULS/'-'/ 
K  •  N 

00  1  I-K,B0 

N  ■>  1*1 

IF  (Atll.EQ. EQULS)  GO  TO  2 

1  CONTINUE 

H   -  1 

2  RETURN 
END 


0001 
0002 
0003 
000* 
0005 
0006 
0007 
0008 
0009 
0010 

28H 

0013 
0014 
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EXPJ 


PURPOSE:        to  calculate  the  exponential  integral  with 
complex  limits. 


METHOD, 


The  exponential  integral  is  defined  as: 


V2 
W12  =   /  e 

V*1 


dv  =  E  (V1)  -  E  (V2)  +  j2mi  , 
1         1 


where  the  integration  path  is  the  straight  line  from  V1  to 
V2  on  the  complex  v  plane  and 


Et(z) 


-00 


dt 


The  integration  path  is  a  horizontal  line  in  the  w  plane 
or  an  inclined  straight  line  from  V1  to  V2  the  v  plane.   The 
integer  n  is  zero  unless  this  path  intersects  the  negative 
real  v  axis  at  a  point  between  V1  and  V2.   When  there  is 
such  an  intersection, 

a)  n  =  1  if  Im(V1)  >  Im  (V2) 

b)  n  =  -  1  if  Im(V1)  <  Im(V2). 

The  term  j2nii  is  caluclated  below  statement  12. 


CALLED  BY 


GGMM 


CALLS  TO 


NONE 
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SUBROUTINE    EXPJ     (Vl.V2.W12)  1 

COMPLEX    EC.fcl5,S, T ,UC,VC,V1 ,V2,W12,Z  2 

OIMENSION    VI21I,     W(21),     D(16),     E<l6)  3 

DATA    V/0.22?8  46  67EOO,0.  11  8  8  9  3  2  1  E  0  1  ,  0.  299  2  73  63  EO  1  ,  0 . 5 77  5 1436E 0  1  ,  0. 9  4 

18374674E01.0. 159826 74 E02,0.93307812E-01,0.49269l74EOO,0. 1215595 4E0  5 

21. 0. 226 99  495E01 .0. 366 76 22 7E01, 0.542 53 366E01, 0.756 591 62 E01, 0.101202  6 

328E02,0.1313  02  32E02,0.166  5440  8E02,0.207764T9E02,0.25  62  38  94E02,0.31  7 

44 075 1 9E 02, 0.3853 06 83E02. 0.48026 08 6E 02/  8 

OAIA    W/0.45896460E00.0.41 70008 3 E 00 ,0 . 1 1 3  3  73 38EOO ,0. 1 03  9  9 1 9  7E-01 ,0.  9 

126101 72 OE-O 3, 0.89854791E-06.0. 2182 3487E00.0. 34221017 EOO, 0.26302758  10 

2E0O,0.12642582E0O,O.4O206  86  5E-Ol,0.8  56  38778E-O2,O.12l24361E-02,O.l  11 

3  1  16  7440E-03,0.64  5  9926  7E-05,0.22  26  316  9E-0  6,0.422743  04E-08,0.392l89  7  12 

43E-10,0.14565152E-12,0.14830270E-15,0.16065949E-19/  13 

DATA    D/0.22495842E02,0.744U568E02,-0.41  43  1  5  76E  0  3  ,  -0.  7  8  754  3  39  E  02  ,  0  14 

1. 11 2 54 744 E 02,0.  16  02 1 76 1 E 03  ,  -0. 23  86  2195E 03, -0. 500946  8 7E0  3 ,-0. 6 848 78  15 

2  54E02,0.12  254  778E02,-0. 10 16 1 976E 02 .-0. 472 1 9 59 1E01 , 0. 79 72 9681 E01 ,-0  16 
3. 2 10695 74 £02, 0.2 2 0464 90E0 I, 0.89 72 8244E01/  1  7 

OATA    E/O.21l63l07E02,-O.3  79  59  767E03,-0.9  74  892  2OE02,O.12  90O6  72E03.O  18 

1. I  79492 26E 02,-0.  I  2  9 1093  IE  03 .-0. 55  705574E 03, 0. 13  524801E02 .0. 1 46967  2  19 

21 E0  3.0.  17949528E 02,-0. 3298 l6l4E 00 ,0.3 10288 36E02, 0.8 165 765 7E01  ,0.22  20 

3  2  36961E02,0.39I24  89  2E02,0.81636  799E01/  21 
I    *    VI  22 

23 

00    12    JIM=l,2  24 

X  «  REAlU)  25 
Y  =  AIMAGI Z)  26 
E15  «  (  .0..0)  27 
AB  «  CABSJ  Z)  28 
IF  (AB.EO.O.I  GO  TO  11  29 
IF  (X.GE.0..AN0.A8.GT.10.I  GO  TO  10  30 
YA  «  ABS(Y)  31 
IF  (X.LE.0..AN0.YA.GT.10. )  GO  TO  10  32 
IF  IYA-X.GE.1 1.5.0ft. YA.&E. 6.5. OR. X*YA.GE. 5.5. OR.X.GE. 3.1  GO  TO  2  33 
IF  IX.LE.-9.)  GO  TO  6  34 
IF  IVA-X.GE.2.5)  GO  TO  7  35 
IF  IX*YA.GE.1.5>  GO  TO  3  36 
N  *  6.*3.»AB  37 
E15  »  1 ./<N-1.I-Z/N**2  38 

1  N  «  N-l  39 

E15  «  1./(N-1.1-Z*E15/N  40 

IF  IN.GE.3)  GO  TO  1  41 

E15  -  Z»E15-CMPLX(.5772l6*AlOG<A8),ATAN2(Y,Xn  42 

GO  TO  11  43 

.  2  J  1  «  1  44 

J2  ■  6  45 

GO  TO  4  46 

3  Jl  •  7  47 
J2  -  21  48 

4  S  '  (  .0..0I  49 
YS  »  Y»Y  50 

51 

00  5  I=J1,J2  52 

XI  =  V(  1 t*X  53 
CF  =  Wl I I/(XI»X1*YS)  54 

5  S  =  S*CMPLX(Xl*CF,-YA*CF)  55 

56 

GO  TO  9  57 

6  T3  =  X*X-Y*Y  58 
T4  »  2.»X*YA  59 
T5  =  X»T3-YA*T4  60 
T6  =  X*T4*YA*T3  61 
UC  «  CMPLXIOt  1  1  )  *0(  12)*X*D<  13)*T3*T5-E(  121*YA-E<  1  3  )  ♦  T4  ,  E  (  It  >*EU2I  62 

1*X*EU3)»T3«T6*D(12)*YA*0(13)*T4)  63 

VC  =  CMPLXt D(  14) »0«  15)*X*D<  16)*T3*T5-E<  15I*YA-E< 1  6  I  •  T4  ,  E  ( 14»*E(15»  64 

1*X*E<16)*T3*T6*OU5)*YA*0(16)*T4)  65 

GO  TO  8  66 

7  T3  =  X»X-Y*Y  67 
T4  =  2.»X*YA  68 
T5  -  X* T3-YA*T4  69 
T6  «  X*  T4*YA»T3  70 
T7  '  X»T5-YA*T6  71 
T8  =  X* T6*YA»T5  72 
T9  -  X»T7-YA*T8  73 
T10  -  X*T6*YA*T  7  74 
UC    =    CNPIXIDI  1) *0<  2 )*X*0(  3)*T3*D<4)*T5*D< 5» *T7* T 9- ( E < 2 ) * YA*E ( 3)*T4  75 

l*EI4)*r6*El5»»T8),EUi*E<2)*X*E<3»*T3*EI41*T5*E<5)*T7*T10*I0<2l*YA  76 

2*0I3)*T4*D(4)»T6»D( 5)*T8I )  77 

VC    =    CMPLX(0(6> »0<7>*X*D(8I»T3*D(9I*T5*D( 10  I • T7* T9- i E < 71*YA* E(8I*T  78 

14*E(9l*T6*E<l0i*T8),E(6)*EI7)*X*E<8)*T3*E(9»«T5*EU0l*T7*T10*(D<7)  79 

2*YA*D<8)*T4*D(9)*T6*D<10)*T8))  80 

8  EC  »  UC/VC  81 
S    =    EC/CMPi_X(X,YA)  82 

9  EX  «  EXP(-X)  83 
T  -  EX*CMPLX(C0S(YA),-S1N« VA)  )  84 
E15  «  S»I  85 
IF  (Y.LT.O.j  E15  •  C0NJGIE151  86 
GO    TO    11  87 

10  E15  »  .409319/1 Z*. 1930441*. 421831/IZM. 026661*. 147126/<Z*2.  567881  *  88 
1.206335E-l/(Z*4.90035»*.107401E-2/(Z*8.  1 82 1 51 ♦. 1 58654E-4/I Z* 1 2 . 734  89 
221*. 317031E-7/<Z*19. 3957)  90 

E15    -    E15»CEXP(-Z)  91 

11  IF    (JIM. EG. I)    H12    -    E15  92 

12  Z    •    V2  93 

94 

Z    »    V2/V1  95 

TH    -    ATAN21A1MAG(Z),REAL(Z) ) -AT AN2 ( AIMAG I V2 ) , RE AL I V2 I) ♦ AT AN2 ( Al MAG  96 

KVll.REAKVm 

AB    ■    ABS(TH)  21 

IF    UB.LT.l.l    TH    -    .0  22 

If  IJHH:eJ:ii!,T?HV:i?238!?8353  iB 

END  |03 
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GANT1 


PUfiPOSE:       to  consider  the  wire  structure  as  a 
transmitting  antenna  and  calculate  the  input  impedance  and 
current  distribution. 


METHOD: 


If  a  wire  antenna  is  driven  by  a  voltage 


generator  v   located  at  one  of  the  current  sampling  points 
i 

1   and  if  displacement  currents  are  neglected,  Ampere's  law 
i 

yields 


V   =  V   F  (1  ) 

m    i  m   i 


where  F  is  the  sinusoidal  expansion  function.  Thus,  the 

excitation  voltages  V   will  vanish  everywhere  except  where 

m 

v   is  not  zero, 
i 


The  DO  LOOP  ending  with  statement  2  uses  the  delta-gap 
model  defined  above  to  determine  the  excitation  voltage 
CJ (I)  for  all  the  dipole  modes.   These  are  stored 
temporarily  in  CG(I).  Then  subroutine  SQROT  is  called  to 
obtain  a  solution  of  the  simultaneous  linear  eguations. 
SQROT  stores  the  solution  (the  loop  currents)  in  CJ  (I)  . 

In  the  DO  LOOP  ending  at  statement  6  ,  the  complex  power 
input  and  input  impedance (s)  are  calculated.  The 
time-average  power  input  (PIN)  is  the  real  part  of  the 
complex  power  input. 

Subroutine  RITE  is  called  to  make  the  transformation 
from  the  loop  currents  to  the  branch  currents.  If  IWR  is  a 
positive  integer,  RITE  will  write  out  the  list  of  branch 
currents. 
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Finally,  GANT1  calculates  the  radiation  efficiency  by 
calling  subrouinte  GDISS  to  obtain  the  time-average  power 
dissipated  in  the  lumped  loads  and  the  imperfectly 
conducting  wire. 


CALLED  BY: 


MAIN 


CALLS  TO: 


GDISS 


RITE 


SQROT 


subroutine  gantl  ( i  a. i b. jnm, imr, 1 1. 1  2, i  3. 1 12, ja. jb . md.n.no.nh. am.c 

..yII.zIi.zlo.zs.Igr6) 

2ld(1).y11.z11.zs.gam,cg(1) 

N0I1I 


i.CJ.CG.CMM.D.EFF .GAM.GG.CGD.SGD.VG 
COMPLEX  VY  ,  CGEN 

COMPLEX  CUttCJUl  iCGDIllfSOOd)  .VGdl.ii. 
OIMENSION  D(I7,  IAIlI.  IBllTi  JftlU;  Ji(l 
DIMENSION  Hill.  12(11.  13(1),  M0(INM,4>, 
COMMON  IWL 


00  3  I-l.N 
CJ(I)  -  1.0, .01 
K  -  JA(  I) 


IA(K| 


00  2  KK> 

KA 

KB  •  IB(K1 

it '  I 

IF  (KB.EQ.I2( 

IF  IKB.EQ. 

CJ(I)  -  CJI 

GO  TO  2 

IF  (KA.EQ. 13(1)1 

JJ  ■  K*NM 

CJ(I)  «  CJ(I)*FI*VGUJ) 

K  -  JB( I) 


).I2f III 
).  Ill  1 1  I 

.  J  (  I  I *F  I  * 


GO  TO  I 
FI  — I. 
VG(JJ) 

FI  — 1. 


3  CONTINUE 


DO  4  I-l.N 
CG(I)  -  CJ(I) 


CALL  SQROT  ( C ,CJ ,0  ,1 12 ,N) 

112  -  2 

VU  -  (  .0..0) 

NNN  ■  N 

IF  (IGRD.GT.O)  NNN  -  (N*IWL)/2 


00  6  I-l.NNN 

NN  -  IA(JB(I)I 

CGEN  -  CG(I) 

IF  (I.LE.tML)  CGEN-CGEN/2. 


0001 
0002 
0003 
000* 
0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0016 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 

mi 


YY  «  CJ(I)*CONJG(CGEN) 


IF  (CABS(VY).LT .l.E-20)  GO  TO  5 
ABS(CGEN7*»2l 
o.d)  r- 
YU  -  Yll*YY 


Zll  -  ( l./YY)»(Ci 
MRITE  (6.8)  NN.Zll 


CONTINUE 

IF  (IMR.GT.OI  WRITE  (6.7) 

CALL  RITE  (IA.IB.INM.IMR.ll.I2.I3.MD.ND.NM.CJ.CG.IGRD! 

GG  -  REAKYll) 

IF  (IGRD.GT.O)  GG-2.»REAL(Y1U 

PIN  -  GG 

CALL  GOISS  (AM.CG.CMN.D.OISS.GAM.NM.SGO.ZLO.ZS) 

PRAO  ■  PIN-OISS 


EFF  «  lOO.^PRAO/PIN 
RETURN 


FORMAT  (SOX. 'ANTENNA  BRANCH  CURRENTS* I 
FORMAT  (lOX.'THE  INPUT  IMPEDANCE  AT  NODE 


1F15.T//1 
END 


•,13.'  ISSF15.7,'  ♦  J', 


0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
00  58 
0059 
0060 
0061 
0062 

88$* 

0065 

>68 
0069 
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GDISS 


PURPOSE:       to  calculate  the  time-average  power 
dissipated  in  the  imperfectly  conducting  wire  and  in  the 
lumped  loads. 

METHOD:  The  time-average  power  dissipated  by  the  wire 
is  calculated  in  the  DO  LOOP  ending  at  statement  1  utilizing 
the  equation  below: 


1 

b  r  * 

P   =  s_       I  I  dl 

d    2na  J 

0 

where  R    is  the  surface  resistance  of  the  wire  and  a  is  the 
s 

radius  of  the  wire. 


The  power  dissipated  by  the  lumped  loads  is  calculated 
by  the  DO  LOOP  ending  at  statement  3.  If  the  wire  is 
perfectly  conducting,  CMM  <0,  the  first  calculation  is 
by-passed. 

CALLED  BY:      GANT1 
CALLS  TO:       NONE 
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SURRGUT 1NE  GOISS  ( AM, CG ,CMM,D, 0  I SS.GAM, NM. SGDtZLO t IS  1  1 

COMPLEX  CG< 1  ) ,SGD(  1  I  ,ZLD( 1 )  , C J  A  ,  C JB  ,  GAM  , it  2 

DIMENSION  Oil)  3 

DATA  PI  /3.  1M59/  * 

DISS  *  .0  5 

IF  (CMM.LE.0.1  GO  TO  2  6 

7 
8 
9 
10 
11 
12 
13 
14 
15 


ALPH 
BETA 
RH    « 

*     REAL(GAM) 

■     A1MAGIGAM) 

REALUSI/<4.*PI»AM» 

00     I     K'l.NM 

DK    ■    Dt  Kl 

OEN    •    CABS(SGD(K  »  )»*2 

EAD    «     EXP(  AlPH*l>K  1 

CAO    «     • EAD*l./EADI/2. 

CBO    *    C0S(BETA*DK) 

SAU    -    OK 

IE     IALPH.NE.O.)     SAO-IEAD 

I* 


•1 ./EAD)/(2.»ALPHl  18 

SBD    =    OK  19 

IF    (BETA. NE. O.J     S8D=SIN<BETA*DK»/BETA  20 

FA    «    RH*<SAO*CAD-SBD»CBD)/DEN  21 

FB    '    2.»RH*(CAD»SBD-SAD*CBD)/DEN  22 

CJA    «    CGIKI  23 

L    »    KtNM  2* 

CJB    «=    CG(L)  25 

1  DISS    «    DISS»FA»<CABS(CJA1**2*CA8S(CJB»»*2)*FB*«REAL(CJA)*REALICJB)  26 
H-AIMAG<CJA»»AIMAG(CJ3)  I  27 

C  28 

C  29 

2  00    3    JM.NM  30 
K    »    J*HM  31 

3  OISS    -     DISS*ftEAL  (ZLO( J)  I*(CABS(CG(J1  ) ♦♦  2 ) *RE AL ( Z ID < K ) ) • ( C ABS C CG( K )  32 
1)*»2)  33 

C  34 

RETURN  35 

ENO  36 
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GFP 


PURPOSE:       to  calculate  the  far-zone  field  of  a 
sinusoidal  electric  monopole. 

METHOD:         If  an  electric  line  source  has  length  d  and 

endpoints  at  (x  ,y  ,z   )    and  (x  ,y  ,z  ),  then  the  coordinates 

111         2   2   2 

of  any  point  on  the  source  are 

x  =  x   +1  cosx 
1 

y  =  y   +1  cosy 
1 

z  =  y   +  1  cosz 
1 

where  cosx,  cosy,  cosz  are  the  direction  cosines  of  the  1 

axis,  and  1  is  the  distance  along  the  source  measured  from 

the  endpoint  (x  ,y  ,z  ).   Let  the  current  distribution  on 

111 

the  monopole  be 


I   sinh  /(d  -  1)  +  I   sinh  yi 
KD  =  —1 2 


sinh  yd 

where  I   and  I   are  the  endpoint  currents.  The  far-zone 
1       2 

field  of  this  source  is 

E   =  (cosx  cosG  cos<b  -  cosy  cos6  sin*  -  cosz  sinQ)  E 

E   =  (-cosx  sin$  +  cosy  cos$)  E 
fl  1 

where 

-yr 

r/e  ygd  /f  < » > 

E  =   —  [  (e    -  g  sinh  /d  -  cosh  yd)     I  e 

1   4nr  (1-g2)  sinh  yd  » 

yd  Kf<2) 

♦  (e    +  g  sinh  yd   -    cosh  yd)    I   e      ] 
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f  <  *  >  =  x   sinG  cos<|>  +  y   sinO  sin<(>  +  z   cosO 
1  1  1 


fC2>  =  x   sinO  cos*  +  y   sinO  sin<i>  +  z   cosO 
2  2  2 

g  =  cosx  sinG  cos<t>  +  cosy  sin©  sin*  +  cosz  cosd 


and  (r#Q,<|>)  are  the  spherical  coordinates  of  the  observation 
point. 

In  this  subroutine  the  range  dependence  has  been 
suppressed.  The  far  field  vanishes  in  the  endfire  direction 
where  GK  =  0.  If  a  ground  plane  is  present  (IGRD  >  0)  the  E 

equation  above  is  decomposed  into  the  x,    yr  and  z  components 

and  the  reflection  coefficients  are  applied  before  E   and  E 

e      <)> 

field  components  are  returned  to  the  calling  program. 
CALLED  BY:      GFFLD 


CALLS  TO:       NONE 
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SUBROUTINE    GFF    < X A.YA .ZA.XB ,YB,ZB.OtCGD, SGO tCTHt STHtCPH, SPH.GAN.ET 
lA.ETl,ET2,EPl,EP2,tGRO,ERRl 
COMPLEX    ERR,RV,RH,RR,EX,EY,EZ,EE 
CCMPLEX    ETl,ET2,EPl,EP2,GAM,ETA 
COMPLEX    GD.CGD.SGO.EGD 
COMPLEX    EGFA.EGFB.EGGD.ESA.ESB 
COMPLEX    CST 
FP    «    12.56637 
XAB    «=    XB-XA 
VAB    -    Y8-YA 
ZAB    <    ZB-ZA 
CA    ■    XAB/O 
CB    «    YAB/D 
CG    -    ZAB/D 

G    -    (CA*CPH»CB*SPH)*STH*CG*CTH 
GK    -    l.-G*G 
ET1    «    ( .0,.0) 
ET2    ■    < .0,.0) 
EP1    -    I  ,0,.0» 
EP2   >    I  .O..OI 
IF    (GK.LT. .001 )    GO    TO    3 
FA    «     <XA»CPH»YA*SPHI*STH*ZA»CTH 
FB   «    (XB»CPH»Y8*SPH)*STH»ZB*CTH 
EGFA    -    CEXP(GAM*FAI 
EGF8    «    CEXPIGAM»F8» 
EGGO    »    CEXPIGAM*G»DI 
CST    «     ETA/(GK*SGD»FPI 
ESA   -    CST»EGF»*(EGGO-G*SGO-CGO» 
ESB   <    CST»EGFB*I  1 ./EGGD»G*SGO-CGD) 
IF    < IGRO.LE.Ol    GO    TO    2 


RV 

RH 

IF 

RR 

RV 

RH 

EX 

EY 

EZ 

EE 

EX 

EY 

EZ 

ESA 

EX 

II 
EE 


*CTH*RRI 


C-l.fOI 

I-l.tOJ 
IGRO.ffl.ll    GO    TO    1 

CSGRT(ERR-STH*STHj 

-(ERR*CTH-RR»/(ERR 

<CTH-RR)/(CTH*RR) 

CA*ESA 

CB*ESA 

CG»ESA 

( EX*SPH-EY*CPH)» (RH-RVJ 

EX*RV»EE»SPH 

£Y»RV-EE*CPH 

-EZ»RV 
— EX»CA-EV*CB»EZ»CG 

CA*ESB 

Mh 

iEX*SPM-EV*CPHI*«RH-RVI 


0001 
0002 
0003 
000* 
0005 
0006 
0007 
0008 
0009 

881? 

0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 

mi 

0035 
0036 
0037 

8811 

0040 
0041 
0042 
0043 
0044 
0045 

88*7 

0048 


EX  ■  EX*RV*EE»SPH 

EY  ■  £Y*RV-EE*CPH 

EZ  -  -EZ*RV 

ESB  — EX»CA-EY*CB*EZ»CG 

T  -  <CA*CPH»CB*SPH)»CTH-CG*STH 

P  -  -CA*SPH*CB*CPH 

ET1  «  T*ESA 

ET2  -  T*ESB 

EP1  -  P»ESA 

EP2  «  P*ESB 

CONTINUE 

RETURN 

END 


885-3 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 


uu 


GFFLD 


PU3P0SE:       to  calculate  the  far-field  for  the  thin  wire 
structure. 

METHOD:        The  far-field  for  the  structure  is  calculated 
from  the  loop  currents.  The  loop  currents  are  either  the 
currents  produced  by  the  transmiiting  antenna  calculations 
of  subroutine  GANT1  or  the  currents  produced  by  an  incident 
plane  wave. 

If  the  incident  field  is  generated  by  a  distance  source 

with  spherical  coordinates  (r  ,0  ,<t>  )  ,  the  excitation 

0   0   0 

voltages  induced  by  a  incident  plane  wave  are 


=  /  F  E 

J     m  i 


V   =  /  F  E   dl 
m 


where 


E   =  E  exp (/  r-r   ) 
i    0  0 


where  E   is  a  vector  constant,  f   is  a  vector  from  the 
0  0 

coordinate  origin  to  the  distance  source,  and  r  is  the 

radial  vector  from  the  origin  to  the  observation  point 


The  field  E   is  generated  by  test  dipole  m  when 
M 

radiating  in  the  homogeneous  medium.  Using  the  vector 

potential,  the  field  at  the  distance  point  (r  ,0  ,4>  )  is 

0   0   0 

iwu  e    o 
E   = /F   exp(/  f.f  )  dl 


/' 


m       Unr        /  ra  ° 
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where  the  radial  component  is  to  be  suppressed.  From  the 
above  equations, 


4nr    /r 

V   = Q-    e   o  e   e 

m      jwu         °   m 


If  an  antenna  gain  calculation  is  desired,  INC  is  set  to 
zero.   PH  and  TH  denote  the  spherical  coordinate  direction 
of  the  distance  observation  point.  The  phi-polarized  (EPPS) 
and  the  theta-poiarized  (ETTS)  components  of  the  electric 
field  intensity  are  returned  to  the  calling  program. 

If  INC  =  1,  a  backscattering  calculation  is  desired.  In 
this  case  PH  and  TH  denotes  the  incident  angles  for  the 
incident  plane  wave.  These  are  also  the  spherical 
coordinates  of  the  distance  source.  The  outputs  returned  to 
the  calling  program  include  absorption,  extinction,  and 
scattering  cross  section  for  each  polarization;  scattered 
electric  field;  and  echo  areas. 

If  INC  =  2,    a  bistatic  calculation  is  desired.  In  this 
case  PH  and  TH  denote  the  spherical  coordinate  of  a  distance 
observer.  Since  this  calculation  uses  the  induced  loop 
currents  (EP  and  ET)  ,  a  backscattering  call  must  preceed 
this  calculation.  The  outputs  returned  to  the  calling 
program  consist  of  the  scattered  electric  field  components 
and  echo  areas. 

EPP(I)  and  ETT  (I)  denote  the  phi-polarized  and 
theta-polarized  far-zone  fields  of  dipole  mode  I  with  unit 
terminal  current.  In  a  backscattering  situation,  the 
excitation  voltages  EP (I)  and  ET (I)  are  obtained  by 
multiplying  EPP  and  ETT  by  the  constant  CJI.  Then  calls  are 
made  to  SQROT  which  stores  the  solution  (the  induced  loop 
currents)  in  EP(I)  and  ET(I).  RITE  is  called  for  the  branch 
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currents  CG  ( J) ,  and  GDISS  is  called  for  the  time-average 
power  dissipated  in  the  imperfectly  conducting  wire  and  the 
lumped  loads.  This  power  is  denoted  PDISS  and  TDISS  for 
phi-polarized  and  theta-polarized  incident  waves, 
respectively. 

In  scattering  problems,  the  incident  plane  wave  has  unit 
electric  field  intensity  at  the  origin.  GGG  denotes  the 
time-average  power  density  of  the  incident  wave  at  the 
origin.  ACSP  and  ACST  denote  the  absorption  cross  sections 
for  the  phi  and  theta  polarizations. 

PIN  and  TIN  denote  the  time^average  power  input  to  the 
wire  structure,  delivered  by  the  equivalent  voltage 
generators  VP  and  VT  at  the  terminals.  PIN  and  TIN  apply  for 
the  phi  and  theta  polarizations,  respectivity.  The 
time-average  power  input  is  regarded  as  the  sum  of  the 
time-average  power  dissipated  and  the  time-average  power 
radiated  or  scattered  by  the  wire.  ECSP  and  ECST  denote  the 
extinction  cross  sections  and  SCSP  and  SCST  denote  the 
scattering  cross  sections. 

The  distance  field  is  calculated  in  the  DO  LOOP  ending 
with  statement  7  for  scattering  situations,  and  in  the  DO 
LOOP  ending  with  statement  9  for  the  antenna  situation. 

The  radar  cross  sections  (echo  areas)  SPPM,  SPTM,  STPM, 
and  STTM,  are  defined  as 

2   2ar 
a  -    limit   Unr   e     S   /  S 
r  -ao  s     i 

where  S   and  S   denote  the  time-average  power  densities  in 
s      i 

the  scattered  and  incident  fields  evaluated  at  the  origin. 

For  an  antenna,  the  following  definition  is  employed  for 
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the  power  gains: 

2   2ar 
G  («,♦)  =  limit   Unr   e    S(r,e,4>)  /  P 

p        r  —  •  i 

where  P  ,    GG#  denote  the  time-average  power  input  and 
i 

S(r,6,$)  is  the  time-average  power  density  in  the  radiated 

field.   GPP  and  GTT  denote  the  power  gains  associated  with 

the  phi-polarized  and  the  theta-polarized  components  of  the 

field,  respectively. 

The  use  of  the  variables  JFLAG  and  KFLAG  are  described 
in  subroutine  SGANT. 


CALLED  BY:      MAIN 

CALLS  TO:       GDISS 

GFF 
RITE 

SQROT 
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lunncA     tiiiitjiincriiMcililicrriili 
COMPLEX    ETA,&AM,CGD( 1  1  , SGO (  l),CG( I) 
DIMENSION    14(1),     16(1),     I  l<  Ml     |2«l»' 
01MENSI0N    311),     X(l),     YI1I,    2(1) 


SUBROUTINE    GFFLO     (  I  A ,  I  B , I NC , INM. IWR,  1 1 , I  2, I  3, 1 1 2, MO.N, ND.NH, AH. ACS  1 

lP,ACST.C.CGO,CG,CJ,CMM.6.ECSP,EdST,EP.ET ,EPP, ETT , EPPS, EP TS.ETP$,Ef  2 

jTS.GG.GPP.GTt.PH.Sdo.SCSP.SCST.S^PMtSPTM.STPM.STTM.TH.X.V.Z.iLOtZS  3 

3,ETA,GAM,ERR,lGRD)  4 

CCMPLEX    err  s 

COMPLEX    CJI.ETl,ET2,EPl.£P2,EPPS,ETTS,EPTS,ET»>S.ZS,Vr>,VT  6 

COMPLEX    c7l),CJ(  1),EP(iJ,Et(1).EPP(1),ETTUT,ZL&(  II  7 

8 

13(1),    NOll),  ND(INM,4)                               9 

DATA    Pi ,TP/3. 14159.6. 28318/    "  11 

CJI    «    -4.*PI/(ETA*GAM)  12 

GGG   ■    REAL! 1 ./ETA)  13 

THR    «    .0174533»TH  14 

CTH    ■    COS(THR)  15 

STH    -     SIN(THR)  16 

PHR    «    .0174533*PH  17 

CPH   «    COS(PHR)  18 

SPH   »    SIN(PHR)  19 

C  20 

00    I     1*1, N  21 

ETT( I)     •    (  .0..0I  22 

1  EPP(I)  »  (.0,.0)  23 
C  24 
C  25 

00  3  K=1,NM  26 

KA  «  IA(KI  27 

KB  -  IB(K)  28 

NGRO  >  IGRO  29 

IF  (K.LE.NM/2)  IGRO— I  30 

CALL  GFF  (XlK.A),Y(KA).Z(KA).X(*Bi,Y(KB).Z(KB),D(K)iCGD(K>,SGO(K),C  31 

lTH,STH,CPH,SPH,GAM,ETA,ETl,ET2,EPl,EP2,IGRO,EKR)  32 

IGRO  »  NGRO  33 

NOK  ■  ND(K)  34 

C  35 

00  3  11*1. NOK  36 

1  -  MO(K,II)  37 
FI  >  1.  38 
IF  (KB. EQ. 12(1))  GO  TO  2  99 
IF  (KB. EO. U(  I))  FI  — 1.  40 

m:. 

GO   TO    3  43 

2  IF    (KA.EO.m  I))    FI— 1. 
EPP(J)    -    EPP( I)»FJ*EP2  45 
ETT(I)    «    ETTlIi»FI*ET2  46 

3  CONTINUE 


HO.CU.U'HII      IjU      IU     C  It 

(KB.EU.1K  I»)    FI  — 1.  40 

m\ :  i?Tp{ll:FF{:iTpt  l! 


APP    «    CABS(EPPS) 
ATT    «    CABSIETTSI 


tl 


EPPS    •    (.0..0)  49 

ETTS    -     (.0..0)  50 

IF    ( INC.EQ.O)    GO    TO    8  51 

IF    (  INC.EQ.2)     30    TO    6  52 

53 

00    *    1=1. N  5* 

ETUI     "     ETT(  I  )*C  J  I  55 

EPU)    ■     EPP(I)»CJI  56 

57 

CALL  SCROT  (C.EP.O, 112, N)  58 

112  *  2  59 

CALL  SOROT  <  C  ,E  T  ,  0,  I  1  2  ,  N)  60 

IF  (IWR.GT.O)  WHITE  (6,10)  PH.TH  61 

IF  (IWR.GT.O)  WRITE  (6,11)  62 

CALL  RITE  ( IA, IB, INM,  IWR.I1 .12,  I  3 ,MD,NO, NM. EP.CG, IGRO)  63 

CALL  G01SS  lAN,CG,CMM,O.POlS,GAM,NN,SGO,ZL&,ZSI  6* 

IF  (IWR.GT.O)  WRITE  (6,12)  65 

CALL  RITE  (I  A, IB, INM,  IWR, f 1,12, I  3 , MO,ND , NM , ET ,CG,  IGRD)  66 

CALL  GOISS  ( AH,CG,CMM,D,TDI S ,GAN, NM, SGO , ZLD , ZS)  67 

ACSP  ■  POIS/GGG  68 

ACST  «  TOIS/GGG  69 

PIN  -  .0  70 

TIN  «  .0  71 

72 

00  5  1=1, N  73 

VP  ■  CJI*EPP( I )  74 

VT  «  CJI*ETT( I )  75 

PIN  ■  PINtREAL (VP*CONJG(EP(  I  )  )  )  76 

TIN  -  TIN»REAL(VT*CONJG(ET(  lit j  77 

78 

ECSP  -  PIN/GGG  79 

ECST  »  TIN/GG&  80 

SCSP  ■  ECSP-ACSP  81 

SCST  -  ECST-ACST  82 

EPTS  «  (.0,.0)  83 

ETPS  -  (.0,.0)  84 

85 

00  7  1-1, N  86 

EPPS  •  EPPS»EP(  l)*EPP(  I  »  87 

EPTS  ■  EPTS»EP(  ll-FTTI  I  88 

ETTS  ■  ETTS*ET(  I  |»ETT( I  89 

ETPS  -  ETPS'ETl I  >*EPP  I  )  90 

l\ 

93 
9* 
95 
96 

97 

Q  Q 

8  00  9  1-l.N  oS 
ETTS  -  ETTS'CJI  I  )»ETT(  I 

9  EPPS  «  EPPS»CJ( I )*fcPP( I) 


100 
101 

102 
103 


GPP  «  *.*PI*APP*APP*GGG/GG  J  XI 

GTT  -  *.»Pl*ATT*ATT*GGG/GG  }2? 

RETURN  ,0° 

10  FORMAT  (10X. 'BRANCH  CURRENTS  ASSOCIATED  WITH  PLANE-WAVE  SCATTERING  108 
1  FOR  THE  INCIDENT  ANGLES,  PHI-',F5.1«'  AND  THET A- ' , F5. 1 // I  10* 

11  FORMAT  I44X. 'CURRENTS  INDUCED  BY  THE  PHI  POLARIZED  WAVE')  }}0 

12  FORMAT  (44X, -CURRENTS  INDUCED  BY  THE  THETA  POLARIZEO  WAVE')  1 
END  lu 
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GGMM 


PURPOSE:       to  calucate  the  mutual  impedance  between  two 
filamentary  monopoles  with  sinusoidal  current  distribution. 

METHOD:         As  stated  in  subroutine  SGANT,  the  mutual 
impedance  of  coupled  dipoles  may  be  expressed  as  sum  of  four 
monopole-raonopole  impedances.  This  subroutine  calculates  the 
mutual  impedance  with  closed-form  expressions  in  terms  of 
exponential  integrals. 

For  skew  monopoles  it  can  be  shown  that  the 
monopole-monopole  mutual  impedance  is  given  by: 


t  -z         z 

i+j      n  mm 

z    =  (-1)     B[e  (F    -e     G    +  e    G) 
ij  J1  12        22 

-t  -z         z 

n  mm 

-  e     (F    -e     G    +e    G   )  ] 
j2  11         21 

where  m=  2/i,  n  =  2/j  and 


B  = 


16  n  sinh  d   sinh  d 
1        2 


The  functions  F   are  defined  by: 

ik 

qz   cos  V 
i 
F    =2  sinh  d   e  E (R   +  qz   cos  9  -  qt) 

ik  *  i     i 


k 
where  q  =  (-1)  .  d   and  d   are  the  lengths  of  the  monopoles 

1       2 

beinq  considered.  The  functions  G    are  defined  as  follows: 

ik 
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G    =  E(R   +  gz   +  g't  -  jg")  +  E  (R   +  gz   *  g't  +  jg") 
ik      2     2  2      2 

-  E  (r   ♦  gz   +  g't  -  jg")  -  E  (R   +  gz   +  g»t  +  jg") 
11  11 

i  k 

where  g  =  (-1)  ,  g*  =  (-1)  ,  and  g"  =  gb  +  g»c  with 

b  =   c   cos  9  and  c  =  d/sin  #.  The  angle  9  is  the  angle  formed 

by  the  apparent  intersection  of  the  two  monopoles.   This 

will  be  discussed  later  in  detail. 


ik 
In  the  above  equation  for  G   ,  t  denotes  the  position  of 

an  observation  point  somewhere  on  monopole  2.  R   and  R   are 

1       2 

the  distances  from  the  endpoints  of  monopole  1  to  this 

observation  point.  Finally,  the  E  functions  are  defined  as 

follows: 


E(a  +  jg 


">  =  e^9"  / 


a  +jg» 
2    -yv 

e    dw 


a  +jg" 


where  a  and  g"  are  real  quantities  with  dimensions  of 

length,  a  is  a  function  of  t,  a   =  a (t  ) ,  a   =  a(t  )  and 

1       12      2 

/  =  jw  yu&.   The  integral  above  is  evaluated  by  subroutine 

EXPJ. 


♦ 

««-- 

-  D 

---► 

i 

I 

I 

S2 

i 

T2 

I 
I 

i    si! 

T              A 

T1 

j 

PARALLEL 

S2 


S1 


T2 


SKEW 


To  explain  the  input  data  for  GGMM,  refer  to  the  above 
figure.   If  the  monopoles  are  parallel,  then  the  new 
coordinate  system  is  defined  such  that  the  new  z  axis  is 
parallel  to  the  monopoles.   The  coordinate  origin  may  be 
selected  arbitrarily.   S1  and  S2  denote  the  z  coordinates  of 
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the  endpoints  of  the  test  monopole,  T1  and  T2  are  the 
coordinates  of  the  endpoints  of  the  expansion  monopole,  and 
D  is  the  perpendicular  distance  (displacement)  between  the 
monopoles.   The  mutual  impedance  of  parallel  monopoles  is 
calculates  in  the  last  part  of  GGMH  below  statement  5. 

For  skew  monopoles,  let  the  test  monopole  s  lie  in  the 
xy  plane  and  the  expansion  monopole  t  in  the  plane  z  =  D.  (D 
is  the  perpendicular  distance  between  the  parallel  planes.) 
If  the  monopoles  are  viewed  along  a  line  of  sight  parallel 
with  the  z  axis,  the  extended  axes  of  the  two  monopoles  will 
appear  to  intersect  at  a  point  on  the  xy  plane.  Let  s 
measure  the  distance  along  the  axis  of  the  test  monopole 
with  the  origin  at  the  apparent  intersection.   S1  and  S2 
denote  the  s  coordinates  of  the  endpoints  of  the  test 
monopole.   Similarily,  let  t  measure  the  distance  along  the 
axis  of  the  expansion  monopole  with  the  origin  at  the 
apparent  intersection.   T1  and  T2  denote  the  t  coordinates 

of  the  endpoints  of  the  expansion  monopole.  Let  s  snd  T  be 

unit  vectors  parallel  with  the  positive  s  and  t  axes, 

respectively.   Then  CPSI  =  s*^  =  cos  |J.  The  monopole  lengths 

are  d   and  d  . 
s      t 

The  output  data  from  GGMM  are  the  impedances  P11,  P12, 

P21,  and  P22.  In  defining  these  impedances,  the  reference 

direction  is  from  S1  to  S2  for  the  current  on  monopole  s, 

and  from  T1  to  T2  for  the  current  on  monopole  t.  In  the 

impedance  P   ,  the  first  subscript  is  1  or  2  if  the  test 
ij 

dipole  has  terminals  at  S1  or  S2  on  monopole  s.   The  second 

subscript  is  1  or  2  if  the  expansion  dipole  has  terminals  at 

T1  or  T2  on  monopole  t.   The  monopole  lengths  d   and  d   are 

s      t 

assumed  positive  in  defining  the  input  data  CGDS,  SGD1  and 
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SGD2, 


For  parallel  monopoles,  CPSI  =  1  or  -1.  S1,  S2 ,  T1,  and 
T2  are  cartesian  coordinates  for  parallel  monopoles  and 
spherical  coordinates  for  skew  monopoles.  For  skew 
monopoles,  the  radial  coordinates  S1,  S2,  T1,  and  T2  tend  to 
infinity  as  the  angle  <J  tends  to  zero  or  n.   Therefore,  if 
the  monopoles  are  within  4.5°  of  being  parallel,  they  are 
approximated  by  parallel  dipoles. 

CALLED  BY:      GGS 

SGANT 

CALLS  TO:       EXPJ 


SUBROUTINE    GGMM    ( S I , S2 , Tl , T2, 0, CGDS , SGD1 , SG02 ,CPS I , ETA, GAM ,P UtP 1 2  1 

1,P21,P22)  2 

DOUBLE    PRECISION    R 1 ,R2 ,DPQ, SI S , TS  1 , TS2 , ST  I , ST2, CD , BD.C PSS, SK , TL 1 , T  3 

IL2,  FDl ,T02tS0I ,DPS1 ,CD,Z0  * 

COMPLEX    C&OS,SGOS,s6oT,SGDl , SGD2 , ET A,GAM , PI  1 , PI  2 , P2 1 ,P22  5 

COMPLEX    CSr,EB,EC,EK.EL,EKL,EG2I,ESl,ES2,ETl,ET2,EXPA,EXPB  6 

COMPLEX    E(2,2),F(2,2)  7 

COMPLEX    EGZ<2,2I ,GM(2),GP(2)  8 

DATA    PI/3.14159/  9 

DSO   «    D»D  10 

SGDS   «    SG01  11 

IF    (S2.LT.Sl)    S3DS    -    -SGOl  12 

SGOT    »    SGD2  13 

IF    (T2.LT.Tl)     SGOT    «    -SG02  14 

IF    (ABSICPSI J.GT..9971    GO    TO   5  15 

ESI    =    CEXP1GAM*S1 J  16 

ES2    •    CEXP(GAM*S2)  17 

ET1    =    CEXP(GAM*T1|  18 

ET2    «    CEXP(GAM*T2)  19 

DD    =    0  20 

OPSI    «    CPSI  21 

101    =    Tl  22 

TD2    =    T2  23 

CPSS    =    0PSI*0PSI  24 

CO    =    DO/DSQRTU.DO-CPSS)  25 

C    *    CO  26 

80    '    C0*DPSI  27 

D        K        Q  n  >  Q 

EB    »    CEXP(GAM»CMPLX<  .0,8)1  29 

EC    '    CEXP<GAM*CMPLX(.0,C) )  30 

DO  1  K«l,2  32 

33 

DO  1  L«l,2  34 

1  E(K,L)  *  <.0,.0I  35 

36 

TS1  «  TDl*T01  37 

TS2  ■  TD2*T02  38 

0P0  -  DD«DD  39 

SI  -  SI  40 

41 

DO  4  1-1,2  42 

F I  -  4-1 1**5  43 

SOI  -  SI  44 

SIS  ■  SDI*S0I  45 

STl  «  2.*SDI*T01*DPSI  46 

ST2  -  2.*SDI*T02*OPSI  47 

Rl  >  0S0RT(0P0*SI S*I Sl-ST 1 )  48 
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DSQRT(0PQ»SIS*TS2-ST21 
EB 

K*l.2 
l-li«*K 
FK*SOI 
EC 


R2 

EK 

00 
FK 
SK 
EL 

DO 

FL 

EKL    =    EK*EL 

XX    «    FK*BD»FL*CO 

TH    «    FL*T01 

IL2    =    FL*TD2 

RRl    =    RHSKUll 

RR2    =    R2*SK»TL2 

CALL    EXPJ    (GAM*CMPLX( RR1,-XX),GAM*CMPLX( RR2.-XX)  ,  EXPA) 

CALL    EXPJ    IGAM*CMPLX(RRl,XXI ,GAM»CMPLX(RR2,XXI,EXPB) 

E(K,L)     =    EIK,L1»FI»IEXPA»EKL*EXPB/EKL) 


L'l.2 
(-  1»**L 


2  EL 

3  EK 


l./EC 
l./EB 


10   •=    SOI*DPSI 
ZC    «    10 

EGZl    -    CEXP(GAM»ZC) 
RRl    «    RUZD-TDl 
RR2    »    R2*ZD-TD2 

CALL    EXPJ     (GAM*RR1,GAM*RR2,EXPB> 
RRl    «    R1-ZD*T01 
PR2    «=    R2-ZO»TD2 

CALL    EXPJ    »GAM»RRl,GAM*RR2,EXPAJ 
F(I,l)    «    2.*SGDS»EXPA/EGZl 
FU,2)     «    2.»SGOS*EXPB*EG2I 
*    SI    -    S2 

CST    «    ETA/(16.*PI»SGDS*SG0T) 

Pll  -  CST*((F<  1,1 )»E<2,2I*ES2-E<1,2)/ES2J«ET2*I-FU,2)-E<2,1 (♦ES2» 
1E(  1,1)/ES2)/ET2) 

PI 2  •  CST»<(-F( 1,1>-E<2,2>*ES2+EII,2>/ES2>»ETU(F<1,2)»E(2,1)*ES2- 
lE(ltl)/ES2>/ETll 

P21  »  CST*U-F<2,ll-E<2,2)*ES1»EU,2)/ES1)»ET2»<F<2,2)»EI2.1)»ES1- 
1E(1,11/ES1)/ET2) 

P22  «  CST«(<F(2| II *E< 2,21 ♦  ESl-EU, 21 /ESI) ♦  ETUI -F( 2,21- E12,1)*ES1* 

RETURN 
5  IF  (CPSI.LT.O.)  GO  TO  6 
TA  ■  Tl 


49 
50 
51 
52 
53 
5* 
55 
56 
57 
58 
59 
60 
61 
62 
63 
6* 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 

11 
89 
90 
91 
92 
93 
94 
95 
96 


TB  «  T2 

GO  TO  7 

6  TA  «  -Tl 

TB  ■  -T2 

SGDT  «  -SGOT 

7  SI  «  SI 

c 

DO  9  1=1,2 

TJ  *  TA 

c 

00    8    J= I , 2 

ZU   »    TJ-SI 

R    »    SQRT(DSQ*Zl J*Z I Jl 

M    '    R»Z  I  J 

IF    UIJ.LT.O.)     W    -    OSQ/CR-ZIJ) 

V    *    R-Z  I J 

IF    (ZU.GT.O.)    V    "    OSQ/tR*ZIJ) 

IF    ( J.EQ. 1)    VI    «    V 

IF     IJ.EO.l)    Wl     =    H 

EGZ(I.J)    *    CEXPIGAM»ZI  J) 

8  TJ  «  TB 

CALL  EXPJ  IGAM*V1 ,GAM*V,GP( I)) 
CALL  EXPJ  (GAM*W1,GAM*M,GN( II) 

9  SI  »  S2 

CST  «  -ETA/I 8.*PI*SG0S*SG0T ) 

Pll  -  CST*(GM(21»EGZl2,2)»GP<2)/EGZt2,2l-CGOS*<GM<l)*EGZU,2l»GP< 1 

1P12G-,ci?*I-GM(2)*EGZ«2,l)-GP(2)/EGZ(2,l)»CG0S*IGM(l)*EGZll,l)»GPI 

-  lp2lE«ZCsfilGMU)*EGZU.2)»GP(ll/EGZ«l,2>-CG0S*<GM(2l*EGZ<2,2)»GP(2 

lP22G-<c{Tii-GM(l)»EGZll,l)-GPm/EGZIl,l)»CGOS*(GN(2)»EGZ«2,l)»GP( 
12)/EGZ(2.1I)) 

RETURN 
END 


97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 

il2 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 

132 
133 
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GGS 


PURPOSE:       to  calculate  the  mutual  impedances  between 
two  filamentary  monopoles  with  sinusoidal  current 
distributions. 

METHOD:        The  monopole-monopole  mutual  impedance  as 
defined  by  SCANT  is  calculated  using  the  equations  defined 
in  subroutine  GNF.  The  endpoints  of  the  axial  test  monopole 
s  are  (XA,YA,ZA)  and  (XB,YB,ZB),  and  the  endpoints  of  the 
expansion  monopole  t  are  (X1,Y1,Z1)  and  (X2#Y2,Z2) .   DS  and 
DT  denote  the  lengths  of  monopoles  s  and  t,  respectively, 
CAS,  CBS  and  CGS  are  the  direction  cosines  of  monopole  s, 
and  CA,  CB  and  CG  are  the  direction  cosines  of  monopole  t. 

The  effects  of  ground  for  vertical  co-linear  monopoles 
are  applied  in  a  slightly  different  manner  than  mentioned 
previously.  As  with  self  impedance  calculations,  the  test 
monopole  and  the  expansion  monopole  are  laterally  displaced 
by  the  wire  radius.  This  lateral  displacement  is  used  to 
determine  the  angle  of  incident.  This  technique  is  applied 
at  statement  8. 

If  INT  =  0,  GGS  calls  GGMM  for  the  closed  form  impedance 
calculations.  Otherwise  GGS  calculates  the  mutual  impedance 
via  Simpson1 s-rule  integration  with  the  following  number  of 
sample  points:    IP  =  INT  +  1.  If  the  monopoles  are  parallel 
with  small  displacement,  GGS  calls  GGMM  to  avoid  the 
difficulties  of  numerical  integration. 

Since  the  point  (X,Y,Z)  of  subroutine  GNF  lies  on  the 
expansion  monopole  t,  T  is  the  integration  variable  and  is 
measured  from  (X1,Y1,Z1).   C1  is  the  current  at  T  for  the 
mode  with  terminals  at  (X1,Y1,Z1),  and  C2  is  the  current  at 
T  for  the  mode  with  terminals  at  (X2,Y2,Z2) .  C  denotes  the 
Simpson' s-rule  weighting  coefficient. 
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Below  statement  7,  GGS  performs  some  analytic  geometry 
in  preparation  for  calling  GGMM.  The  remainder  of  this 
section  is  concerned  with  this  preparation. 

Let  s  denote  a  unit  vector  in  the  direction  from 

(XA#YA,ZA)  toward  (XB,YB,ZB) .  Also  let  T  denote  a  unit 

vector  from  (X1,Y1,Z1)  toward  (X2,Y2,Z2).  Then  s«£  = 

cos  6   =  CC  where  Q    is  the  angle  formed  by  the  axes  of  the 

two  monopoles.  Let  monopole  s  lie  in  one  plane  P   and 

s 

monopole  t  lie  in  another  parallel  plane  P  .  CAD,  CBD  and 

CGD  are  the  direction  cosines  of  the  unit  vector  d  =  t"  x  s  / 

sin  d   which  is  perpendicular  to  both  planes.  To  obtain  the 

distance  DK  between  the  two  planes,  a  vector  E   is 

11 

constructed  from  (XA,YA,ZA)  to  (X1,Y1,Z1)  and  take  DK  = 


I  .d\ 
1 1 


A  line  is  constructed  from  (X1,Y1,Z1)  to  the  test 
monopole,  such  that  the  line  is  perpendicular  to  the  test 
monopole.  SZ  denotes  the  s  coordinate  of  the  intersection  of 
this  line  with  the  test  monopole,  and  the  cartesian 
coordinates  of  this  intersection  are  XZ,  YZ,  and  ZZ.  The 

direction  cosines  of  s  x  3  are  CAP,  CBP,  and  CGP. 


From  the  point  (X1,Y1,Z1)  in  plane  P  ,  a  line  is 
constructed  perpendicular  to  the  point  (XP1 , YP1,ZP 1 )  in  the 
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plane  P  .   This  line  is  parallel  with  Z   and  has  length  DK. 
s 

Let  B  represent  a  vector  from  (XZ,YZ,ZZ)  to  (XP1 , YP 1 , ZP1) . 

P1  denotes  E  (s  x  3) .  S1  and  T1  are  defined  in  subroutine 
GGMM. 


CALLED  BY: 


SGANT 


CALLS  TO: 


GGMM 


SUBROUTINE    GGS    (  XA .YA, ZA.XB.YB.ZBiXl .Yl . II . X2.Y2 1 Z2t AM.DSfCGDSt SGD 
XS.OT.S&OT, INI,ETA,GAH,PlI.Pl2jPai,P2i,EftR, 1GR6) 
COMPLEX    EXl,EYl,EX2,EY2,EZltEZ2 

COMPLEX    Pll.P12tP21,P22,EJA,EJB,EJltEJ2,ETA,GAM,Cl,C2,CST 
COMPLEX    EGO,CG0S,SGDStSGOT,ERl,ER2,ETl,ET2 
CCMPLEX    ERR 
COMPLEX    EEtEXX.EYY 
CCMPLEX    PP.PX.PY.PZ 

CCMPLEX    RR1,RR2,RR3,RR*,RH1,RV1.RH2,RV2,RH3,RV3,RH*,RV* 
DATA    FP/12. 56637/ 
CA    •    IX2-X1I/0T 
CB    «    (Y2-Y1I/DT 
CG    «    (Z2-Z11/0T 
CAS    «    ( XB-XAI/OS 
CBS   ■    1 YB-YA1/0S 
CGS   -     1 ZB-ZAJ/DS 
CC   ■    CA*CAS*CB 
IF    UCG.LE..00 

IF    IABSICCI.GT..997)    GO    TO    6 
L    SZ    »    (Xl-XA)*CAS»(Vl-VAI*CBS»(2l-ZA)*CGS 
IF    ( INT.LE.OI    GO    TO    7 
INS    «    2»(INT/2) 
IF    (  INS.LT.2)    INS    -    2 
IP  «    INS*1 
OELT    *    OT/INS 
T    -    .0 
DSZ    «    CC'OELT 

I.OfOI 

(  .0..0I 

<  .0..0J 

( .0..0) 

AN*  AM 

-I. 


CA»CAS*CB»CBS*CG*CGS 

3».AN0.(CG5.LE..003).AN0.{ IGRD.GT.OSI    GO    TO    1 


Pll 
P12 
P21 
P22 

AMS 
SGN 


00    5 

U\i 

XXZ    - 

YVZ 

III 


IN«l,IP 
SZ 

SZ-DS 

X1»T*CA-XA-SZ*CAS 
Y1»T*CB-VA-SZ*CBS 
Z1»T*CG-ZA-SZ*CGS 
RS   ■    XXZ**2»YVZ**2*ZZZ*»2 
Rl    ■    SQRT«RS*ZZ1*»2» 
EJA    -    CEXPI-GAH*R1I 
EJ1    -    EJA/R1 
R2   -    S(1RT(RS*ZZ2**2I 
EJB  «    CEXP(-GAM*R2) 
EJ2  -    EJB/R2 


0001 
0002 
0003 
000* 
0005 

888t 

0008 
0009 

88K 

0012 
0013 
0014 
0015 

OOlo 
0017 

oola 

0019 
0020 
0021 
0022 
0023 
002* 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
003* 
0035 
0036 
0037 
0038 
0039 
00*0 
00*1 
00*2 
00*3 
00** 
00*5 
00*6 
00*7 
00*8 
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ER1  »  EJA*SG0S<-ZZ1*EJI*CGDS-ZZ2»EJ2 

ER2  «  -EJB*SGDS»ZZ2»EJ2*CGDS-ZZ1*EJ1 

FAC  ■  .0 

IF  (RS.GT.AMSI  FAC  -  < CA*XXZ»CB*YYZ»CG*ZZZ) /RS 

ETl  -  CC*<EJ2-EJ1»CGDSI»FAC»ER1 

ET2  «  CC*<EJ1-EJ2*CGD$J»FAC*ER2 

IF  (IGRO.LT.O)  GO  TO  4 

RV1  «  t-1.,01 

RHl  -  1-1 . ,01 

RV2  -  1-1. ,0) 

RH2  ■  I-1..01 

IF  1IGRO.E0.1)  GO  TO  2 

XG1  «  X1*T*CA-XA 

YG1  ■  Y1»T*CB-YA 

ZG1  «  Z1»T*CG-ZA 

XG2  »  X1*T*CA-XB 

YG2  -  Y1»T*CB-YB 

ZG2  «  Z1»T»CG-ZB 

RG1  •  SQRT(XG1*XG1»YG1*YGI) 

RG2  -  SQRTJXG2»XG2»YG2*YG2) 

TTl  -  ATANIRG1/ZG1 ) 

TT2  «  ATAN(RG2/ZG2) 

CTHl  «  C0SITT1I 

SSTHl  »  SIN( TTl )*SIN(TT1) 


CTH2 

SSTH2 

RR1 

RHl 

RV1 

RR2 

RH2 

RV2 

RG  ■ 

CPH 

SPH 

IF  ( 

CPH 

SPH 

EXX 

EYY 

EE 


EX1 
EYl 
EZ1 
ETI 
EXX 
EYY 
EE  > 


'  C0S1TT2) 

-  SIN<TT2»*SIN(TT2» 

CSQRT1EKR-SSTH1I 

(CTH1-RR1)/ICTHI«-RR1I 

-IERR»CTHl-RRll/<ERR*CTHl*RRU 

CSQRT(ERR-SSTH2> 

(CTH2-R«2»/(CTH2*RR2» 

-( ERR*CTH2-RR2)7(ERR»CTH2»RR21 
SORTUXB-XA)*(XB-XA»*(VB-YA)*(YB-V*n 

0 

0 


RG.LT.l.E-32)    GO   TO   3 

(XB-XA1/RG 

I YB-YA1/RG 

ET1*CAS 

ET1*CBS 
CEXX*SPH-EYY»CPH»*<RH1-RV11 

£XX*RV1*EE»SPH 

EYV*RV1-EE*CPH 

-ET1*RVI*CGS 
-EX1*CAS-EV1*CBS»EZI*CGS 

ET2*CAS 

ET2*CBS 
(EXX*SPH-EYY«CPHI*IRH2-RV2» 


0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
0070 

°0°072 
0073 

887$ 
0076 
0077 
0078 
0079 
0080 
0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 


C-l. 


EX2    «    EXX*RV2»EE*SPH 

EY2   •    EYV*RV2-EE*CPH 

EZ2    »    -ET2*CGS*RV2 

ET2    «-EX2*CAS-EY2*CBS»EZ2»CGS 

4  C   -    3.»SGN 
IF    (IN.EQ.1.QR.IN.EQ.IP1 
EGO   «    CEXP(GAM*(OT-Tn 
CI   ■    C*lEGD-l./EGD»/2. 
EGO   -    CEXP(GAM»T1 
C2    ■    C»(EG0-l./EGDI/2. 
Pll    ■    PU»ET1*C1 
P12   «    P12»ETl*C2 
P21    «    P21»ET2»C1 
P22   -    P22*ET2*C2 
T   -    T»DELT 
SZ   -    SZ*OSZ 

5  SGN   «    -SGN 


CST 

Pll 

P12 

P21 

P22 

RETURN 

SZl    -    <X1-XA)*CAS*<Y1-YA»*C8S*IZ1-ZA»»CGS 


-ETA*DELT/(3.»FP*SGDS*SGDT) 

CST*Pll 

CST*P12 

CST»P21 

CST*P22 


0R1 

liiz 

0R2 
12) 

000 


SQRTUX1-XA-SZ1«CAS»««2M 


VA-SM*CBS)**2+(Zl-ZA-SZl»CGSJ»* 


SZUOT»CC 
SQRTUX2-XA-SZ2*CAS1*«2»«Y2-VA-SZ2«CBS)**2»«Z2-ZA-SZ2*CGS1** 


<ORl*OR2)/2. 
F    <ODD.GT.20.*AM.AND.INT.GT.OI    GO   TO   1 
IF    (ODD. LT. AMI    DDD    -    AM 

CALL    GGMM    < .0, OS , SZl ,SZ2 .OOO.CGDS, SGOS, SGDT , 1. ,ETA ,GAM,P 11 ,P 12, P21 
1.P22) 
IF    (IGRD.LE.l)    RETURN 
IF    (IGRO.GT.l)    GO    TO   8 

ss  ■  sort  (1. -cocci 
cao  -   (cgs*cb-cbs*cg)/ss 
cbo  «  icas*cg-cgs*caj/ss 
cgd  •  <cbs*ca-cas*cb»/s$ 

OK  -  (X1-XAI*CA0»(Y1-VaI*CBD*(Z1-ZA)*CGD 

OK  •  ABS(OK) 

IF  <OK.LT.AM»  OK  -  AM 

XI   -    XA»SZ*CAS 

YZ  -  YA*SZ*CBS 


0097 
0098 
0099 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
0110 
0111 
0112 
0113 

811* 

0116 
0117 
0118 

0121 
0122 
0123 
0124 
0125 
0126 
0127 
0128 
0129 

8139 

0132 
0133 
0134 
0135 
0136 
0137 

siii 

0140 
0141 

8Ui 

0144 
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Zl   «  ZA*SZ*CGS  01*5 

XP1  -  X1-0K*CAD  0146 

VP1  ■  ri-DK*C90  0147 

ZP1  «  Z1-DK*CGD  0148 

CAP  =  CBS*CGD-CGS*CBD  0149 

CBP  «  CGS*CA0-CAS*CGD  0150 

CGP  *    CAS*CB0-CBS*CAD  0151 

PI  »  CAP*(XP1-XZJ»CBP»(VPI-YZ)»CGP*(ZP1-ZZ1  0152 

Tl  -  Pl/SS  0153 

SI  «  T1*CC-SZ  0154 
CALL  GGMM  ( S 1 , 3  H-DS , r I , TH-OT ,0K ,CGDS . SGOS, SGOT ,CC ♦ ETA  ,6AM, PI  I ,P12 ,     0155 

1P21.P22)  0156 

RETURN  0157 

0158 

S  AHS  «  AM*AM  0159 

RG  *  IX1-XA»*IX1-XA)»(Y1-YAJ*IY1-YA)  0160 

IF  (RG.LT.AMSJ  RG  «  AMS  0161 

OG  «  SQRT((Z1-ZA»*(Z1-ZA»«-RG»  0162 

CPH  •  ABSIZ1-ZAJ/DG  0163 

SSPH  *  RG/IOG*DGI  0164 

RR1  ■  CSQRTIERR-SSPH)  0165 

RVl  »  -IERR*CPH-RR1)/IERR*CPH«-RR1)  0166 

Pll  «-Pll*RVl  0167 

RG  «  (X1-XB)*(X1-XB)«-(Y1-YBI*(Y1-YB)  0168 

IF  IRG.LT.AHS)  RG  >■  AMS  0169 

DG  «  S0RTI(Z1-ZbT*(Z1-ZB1»RG)  0170 

CPH  «  ABS<  Zl-ZBl/DG  0171 

SSPH  «  RG/(0G*0GI  0172 

RR1  -  CSQRTIERR-SSPH)  0173 

RVl  «  -CERR»CPH-RRI)/(ERR*CPH»RR1)  0174 

P12  »-P12*RVl  0175 

RG  -  IX2-XA)*(X2-XA)*(Y2-VAI*IY2-YAI  0176 

IF  IRG.LT.ANS)  RG  ■  AMS  0177 

OG  -  SQRT((Z2-ZAI*(Z2-ZA)+RG)  0178 

CPH  ■  ABSIZ2-ZA1/0G  0179 

SSPH  «  RG/(0G*0G)  0180 

RR1  «  CSQRT(ERR-SSPHI  0181 

RVl  ■  -<ERR*CPH-RR1)/(ERR*CPH»RRI)  0182 

P21  --P21*RV1  0183 

RG  -  IX2-XB)*(X2-XB)»(Y2-VB)*(V2-VB1  0184 

IF  (RG.LT.AMS)  RG  -  AMS  0185 

OG  -  SQRT(<Z2-ZB>*(Z2-ZB)»RG)  0186 

CPH  ■  ABS< Z2-Z81/DG  0187 

SSPH  -  RG/?DG*DGI  0188 

RR1  -  CSQRTIERR-SSPH!  0189 

RVl  -  -<ERR*CPH-RR1)/(ERR*CPH*RR1|  0190 

P22  — P22*RV1  0191 

RETURN  0192 


END 


0193 
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GNF 


PURPOSE:       to  calculate  the  near-zone  electric  field  of 
a  sinusoidal  electric  monopole. 


z 
i 

i 
i 
i 


E   i 
z 


I" "":;^  (e'2) 

2  --'- 


+   E 


METHOD 


An  electric  line  source  is  located  on  the  z 


axis  with  endpoints  at  z  and  z   as  shown  in  the  above 

12 

figure.  Let  the  electric  monopole  have  the  following  current 
distribution: 

I   sinh  /(d  -  1)  +  I   sinh  yi 

1(1)  =  _^_J ? 

sinh  yd 


where  I   and  I   are  the  endpoint  currents,  y  is  the  complex 

12 


propagation  constant  of  the  medium,  d  =  z   -  z   is  the 

source  length.  The  cylindrical  components  of  the  field  are 

E(0)  =0  and 

„  -KR        -/R 

E(e)  =  - [  (I   e    i  -  I   e    2)  Sinh  yd 

Une  sinh  /d     *  2 

-yR 
♦  (I   cosh  yd  -  I  )  e   l  cos  6 

12  1 

-yR 

+  (I   cosh  yd  -  I  )  e   2  Cos  9  ] 

2  1  2 
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77  -/R 

E(z)  = [  (I   -  I  cosh  /d)   e 

an  sinh  yd  *     2  R 

2 

-/R 
+  (I   -  I    cosh  yd)    e    *  ] 
2     1  R    J 

1 

where    is  the  intrinsic  impedance  of  the  medium  and  where 
(€,$,z)  denote  the  cylindrical  coordinates  in  a  coordinate 

system  centered  at  the  endpoint  of  z  . 


These  expressions  exclude  the  field  contributions  from 
the  point  changes  at  the  endpoints  of  the  line  source,  since 
these  charges  disappear  when  two  monopoles  are  connected  to 
form  a  dipole. 

Let  the  coordinate  s  measure  distance  along  the  test 
monopole  with  the  origin  at  (XA,YA,ZA) .  From  any  point 
X,Y,Z,  a  line  is  constructed  perpendicular  to  the  monopole. 
SZ  denotes  the  s  coordinate  of  the  intersection  of  this  line 
with  the  monopole.  The  length  of  the  line  is  the  radial 

2 

coordinate  6,  and  RS  denote  6  .  R1  and  R2  are  the  distances 
from  (XA,YA,ZA)  and  (XB,YB,ZB)  to  the  point  (X,Y,Z) . 

In  the  statements  above  statement  1 ,  the  above  equations 
are  solved;  and  after  statement  1,  the  cartesian  components 

(E  #E  ,E  )  of  the  field  are  determined.  If  a  ground  plane  is 
x   y   z 

present  (IGRD>0)  the  reflection  coefficients  are  applied  to 

the  cartesian  components  before  returning  to  the  calling 

program. 

CALLED  BY:      GNFLD 


CALLS  TO:       NONE 
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ex 


CAS  =  I XB-XAJ/DS 

CBS  «  ( YB-YAI/OS 

CGS  =  »  ZB-ZAI/DS 

SZ  =  JX-XAI*CASMY-YA)*CBS*  (Z-ZA)*CGS 

III    '    SZ 

IU   «  SZ-OS 

XXZ  *  X-XA-SZ*CAS 

YYZ  =  Y-YA-SZ*CBS 

III    '    Z-ZA-SZ*CGS 

PS  =  XXZ»*2»YYZ«*2*ZZZ**2 

Rl  «  SQRT(RS»ZZl*»2) 

EJA  •  CEXPI-GAM*R1) 

EJ1  -  EJA/R1 

R2  »  SORTIRS»ZZ2*»2) 

EJB  *  CEXP<-GAM*R2) 

EJB/R2 

EJ2-EJl*CGOS 

EJl-EJ2*CGOS 

I .0. .0) 


EJ2 
ESI 
ES2 
ER1 
ER2 
AMS 


I .0, .0) 
-     AM»AM 
IF    (RS.LT.AMSJ    GO    TO    1 
CTH1    «    ZZ1/R1 
CTH2    «     ZZ2/R2 
ER1 


ER2 
1  CST 
EX1 
EY1 
EZ1 
EX2 
EY2 
EZ2 


( EJA*SGDS*EJA*CGDS*CTH1-EJB*CTH2J/RS 
<-EJB*SG0S*EJB»CG0S*CTH2-EJA*CTHl)/RS 
ETA/14. »PI*SGOSI 


CST»<ESI*CAS»ER1*XXZI 
CSr*IESl*CBS»ERl»YYZ» 
CSr*(ESl*CGS»ER l*ZZZ» 
CST»(ES2*CAS»ER2*XXZ' 
CST*(ES2*CBS»ER2*YYZ 
-   CST*(ES2*CGS*ER2»ZZZ 

IF    (1GRD.IE.0)    RETURN 

RV1    -    t-1.,01 

RH1    -    I-l.tOI 

RV2   «    I-l.iOJ 

RH2    «    t -1 . tO» 

IF    IIGRO.EO.l)    GO    TO    2 

Rl    " 

R2 

TH 


i«.     "     l  —  itivi 

=    IIGRO.EO.l)    GO    TO    2 

L  «  SQRTC!XA-X)*(XA-X)HYA-Y)*(YA-YI) 
2  »  SQRTI  (XB-XI*IXB-X)*(YB-Y)MYB-Y>  t 
11    -    ATANIR1/C ZA-Z) ) 


1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
1* 
15 
16 
IT 
18 
19 
20 
21 
22 
23 
2* 
25 
26 
27 
28 
29 
30, 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 


>  ATAN(R2/(ZB-Z1» 
'    CSQRT(E**-SIN(TH1J*SIN(THIII 
■    CS0RT(ERR-SIN(TH2»*SIN(TH2») 

>  -IERR»COS( TH1 l-RRl I / ( ERR»COSITHl » *RR1 ) 
'     (C0SITH1 )-*Rll/(COS( THl)»RRl| 

>  -( ERR*COS( TH2)-RR2l/JERR*COStTH21»RR2l 

>  <COS(TH2l-RR2)/(COS(TH2)»RR2) 
SORT<IXA-XB>*IXA-XB)»IYA-VB)*<YA-YBII 


TH2 
RR1 
RR2 
RV1 
RH1 
RV2 
RH2 

2  RG  « 
CPH  «  0 
SPH  «  0 

IF  (RG.LT.1.E-32I  GO  TO  3 
CPH  ■  (XB-XAI/RG 
SPH  «  <YB-YA)/RS 

3  EE  ■  IEX1*SPH-EY1»CPHI*<RH1-RV1I 
EX1  ■  -EXI»RV1*EE*SPH 

EYl    »    -EYl»*vi-EE*CPH 

EZ1    -    EZ1M-RV1I 

EE  •  <EX2*SPH-EY2«CPH)*<RH2-RV2I 

EX2  •  -EX2*Ry2*EE*SPH 

E22  -  £22 

RETURN 
END 


rEY2*RV2-EE«CPH 


49 

i. 

52 

\\ 

55 

56 
57 

58 
59 
60 

62 
63 
64 
65 

» 
68 
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GNFLD 


PURPOSE:       to  calculate  the  near-zone  electric  field 
intensity  at  a  given  point. 

METHOD:        This  subroutine  calls  GNF  for  the  near-zone 
field  of  each  wire  segment,  and  sums  over  all  segments  to 
obtain  the  near-zone  field  of  the  wire  antenna.  FI  is  used 
in  a  manner  similiar  to  FI  of  subroutine  SGANT.  CJ(I)  is  the 
loop  currents  calculated  by  subroutine  GANT1. 

The  use  of  the  variables  JFLAG  and  KFLAG  are  described 
in  subroutine  SGANT. 


CALLED  BY 


MAIN 


CALLS  TO: 


GFF 


SUBROUTINE    GNFLOilA, IB, INM, 11.12. I3,M0,N,N0,NH,*H,CGO,SGD,ETA, GAM 
l.cj.d.x.y , z.xp.vp, zp.ex. ey.ez, igro.err) 
Complex  ex,ev,ez,exi,eyI,ezi,ex2,ey2,eZ2,eta,gam 


LI 


loPUllf! 


COMPLEX    _ 

COMPLEX   CJ(1 


12(11.    13(11.    0(11,    XCli,    Yd),    211 


DIMENSION  MDdNM,*),  N0(  1  I 
DATA  PI. TP/3. 1*159, 6. 28318/ 
EX  ■  < .0, .0) 
EY  ■  I.0..0) 
E2  ■  (.0..0) 


00 
KA 
KB 
NGRD 


K-l.NM 

IA(K) 

IB(KI 

IGRO 


IF  (K.LE.NM/2)  IGRD--1 
CALL  GNF  (XIKA1,YCKA>.Z(KAI,X(KB).Y(KB),Z(KBI,XP.VP,ZP,AN,DIK),CGD 
l(KI,SG0(K>fETA,6AM,ExI,EVl,£Zl,EXi,EY2,£22,IGRO,ERRj 
IGRD  ■  NGRO 
NOK  ■  NOIKI 

00  2  1 1 "1. NOK 

1  -_M0(K,II) 

IF  (KB^EQ. 12(1)1  GO  TO  1 
IF  (KB. EO. 11(1)1  FI  —  1, 
EX  -  EX*FI*EX1*CJ( II 
EY  ■  EY*F|*E) 
E2  -  EZ»FI*E< 
GO  TO  2 
L  IF  (KA.EQ.I3d))  FI  —  1, 
EX  •  6X*FI*EX2»CJ( I) 
EY  -  EY*FI*EY2*CJd) 

.  l^:."?F'*E£2»cj(i 

!  CONTINUE 

RETURN 
ENO 


1 1  1 1  i  r  !■" 

;xi*cj(ij 


0001 
0002 
0003 
000* 
0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 
001* 
0015 
0016 
0017 

88it 

0020 
0021 
0022 
0023 
002* 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
003* 
0035 
0036 
0037 
0038 
0039 
00*0 
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LEFT 


PURPOSE:       to  determine  position  (location)  of  the  left 
paren  symbol  on  the  input  data  card. 

METHOD:        The  character  search  begins  in  the  column 
passed  to  the  subroutine.  On  returning  to  the  calling 
program  the  argument  passed  is  the  column  following  the  left 
paren  symbol. 

CALLED  BY:      READ 
CALLS  TO:       NONE 


SUBROUTINE  LEFT  (Nl  I 

COMMON  /A/  A(80)  2 

OATA  PLEFT/M'/  3 

K  -  N  * 

5 

oo  i  i-K,eo  6 

N  -  I ♦!  7 

IF  (A(I).EQ.PLEFT)  GO  TO  2                                                     8 

1  CONTINUE  9 

10 

N  >  1  11 

2  RETURN  12 
END  13 
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LINECK 


PURPOSE 


to  insert  grid  characters  on  the  polar  plot. 


METHOD:        The  peroid  character  (ISM  (2))  is  inserted  in 
the  proper  position  in  the  statements  above  statement  4.  In 
the  statements  after  statement  4,  the  grid  numbers  labels 
are  inserted  on  the  horizontal  axis. 


CALLED  BY; 


POLPLOT 


CALLS  TO 


NONE 


SUBROUTINE  LINECK  IX, Y) 

THIS  SUBROUTINE  INSURES  ALL  GRID  CHARACTORS  LIE  ON  THE  POLAR  GRID 

CCMMON  ISYM.LlNE 

INTEGER  Y 

DIMENSION  ISYMI14),  LINEU30) 

IF    (Y.EQ.O)    GO    TO    3 

K    =    0 

IF    IX. LT. 10. 01    GO    TO    5 

SET    UP    AREAS    OF    "PERIOD"    POLAR    GRID    POINT    CHARACTERS 

I    »    INTIXI 

I    ■    IABSI I ) 

I    '    ABS(X) 

IF     IIZ-U.GT.0.51     l-l»l 

1  IF    (Z.LT.10.0.3R.Z.GT.U1.0)    GO    TO    2 
LINEUI    ■     ISYM12) 

LINEI60)     -     ISYM13) 

LINE162)    »    I  SYMI  3  J 

K    «    K*\ 

IF    IK.E0.2)    GO    TO    2 

I    '    12?-I 

GO    TO    1 

2  LINEI61 )  ■  ISYMI 2) 
IF  IY.NE.0I  GO  TO  5 

3  00  U=  11,111 
LINE(K)  «  ISYM12) 

4  CONTINUE 

FILL  IN  GRID  NJM8ER  LABELS  ON  HORIZONTAL  AXIS 


LINEUI  ) 
LINEI20) 
L INEl 21  J 
LINEI22) 
LINEI30) 
LINE  I  31  ) 
LINE! 321 
L  INEl 40) 
LINEUI  > 
L1NEI42) 
LINEI50) 
LINE! 51  I 
LINEI52) 


ISYMI 7) 
ISYMI 10J 
ISYM15J 
ISYMI 111 
ISYMI 9) 
ISYMI 5) 
ISYMI 11) 
ISYMI 8) 
1SYMI5) 
ISYMI 11) 
ISYMI 7) 
ISYMI5) 
ISYNIll) 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 


LINE (61  ) 

LINEI70) 

L  INEl  71  ) 

LINEI72) 

L INEl 80) 

LINE181  ) 

LINEI82) 

LINE  I  90) 

LINEI91) 

LINEI92) 

LINEI100) 

LINEU01I 

LINEI102) 

LINEI111I 

CONTINUE 

RETURN 

END 


ISYMI 1 ) 
ISYMI 7) 
ISYM15) 
ISYMI  U) 
ISYMI 8) 
ISYMI 5) 
ISYMI  11) 
ISVMI9) 
ISYMI5) 
ISYMI  11) 
'  ISYMI10) 
'  ISYMI5) 
>  ISYM(ll) 
■  ISYMI7) 


49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
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NUMB 


PURPOSE 


to  place  degree  numbers  on  the  polar  plot. 


METHOD:  The  current  line  which  is  being  printed  is 
passed  to  the  subroutine  in  the  calling  argument.  If  this 
line  contains  degree  numbers,  these  numbers  are  placed  in 
the  correct  position  by  the  IF  statements. 


CALLED  BY: 


PTPLOT 


CALLS  TO: 


NONE 


SUBROUT INE    NUMB     (V) 

THIS    SUBROUTINE    PUTS    DEGREE    NUMBERS    ON    POLAR    GRID 

COMMON     ISYM.LINE 

INTEGER    Y 

OIMENSION     ISYMI14),    LINE(130I 


IF    (Y.NE 
LINE (331 
LINEI34) 
LINEI35) 
LINEI87) 
LINEI88 ) 
LINE! 89) 
IF    (Y.NE 
LINE! I?)    • 
L I NE (  I  3  »    - 
L!NE( 14)    • 
LINEI108)     = 
LINE! 109) 


37)  GO  TO 
'  ISYMI 7) 
=  1SYMI8) 
«     ISYMI6) 

■  ISYMI6) 
»  ISYMI12) 

■  ISYMI6) 
21)  GO  TO  2 
•  1SYMI7) 

■  ISYM(ll) 
ISYMI6) 

ISYMI6) 
ISYM(9) 


LINEI110I  =  1SYMI6) 
IF  IY.NE.0)  GO  TO  3 
LINEI7)  =  ISYMI7) 


LINEI8) 
LINE  I  9)  = 
L INEI 11 31 
LINE! 114) 
LINE! 11 51 


ISYMI |3) 
ISYMI6) 
=  ISYMI6) 
«  ISYMI6) 
=  ISYMI6) 


3  I F  IY.NE.-21  )  GO  TO  4 


LIME! 12  I 
LINE! 1 3) 
LINE! 141 
LINEU08) 
LINEI 109) 
LINE! 1101 


IF  IY.NE.-37) 
LINEI33 
LINEJ34 
LINFI35 

L INE(d7 

LINEI 68 

LINEI89 

CONTINUE 

RETURN 

ENO 


ISYMI8) 
ISYMI 7) 
ISYMI6) 
'  ISYMI9) 
'  ISYMI9) 
ISYMI6I 


GO  TO  5 
ISYMI 8 i 
IsymJ 10) 

ISYMI6 
JSYMI9 
I SYMI6 
ISYMI6 


3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 

11 
39 
40 
41 
42 
43 
44 
45 
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NUMBER 

PURPOSE:       to  convert  alpha-numeric  numbers  to  floating 
or  fixed  point  numbers. 

METHOD:        After  initially  determining  the  sign  of  the 
number,  the  DO  LOOP   ending  at  statement  6  scans  each 
character  beginning  at  N1.  The   DO  LOOP  ending  at  statement 
3  terminates  the  outer  DO  LOOP  if  the  character  being 
compared  is  not  an  alpha-numeric  number.  The  DO  LOOP  ending 
at  statement  5  converts  the  alpha-numeric  number  to  an 
actual  number.   Below  statement  7,  the  multiplier  correction 
is  applied  to  the  floating  point  number  before  returning  to 
the  calling  program. 

CALLED  BY:      READ 
CALLS  TO:       NONE 


SUBROUTINE  NUMBER  (Nl,N2iX,IX>  1 

COMMON  /A/  A180)  2 

DIMENSION  B(  101  3 

DATA  8/  '0'  ,'  I  ■  ,  «2«  ,'3'  ,'4«  t'5't  'fi'T,  •8,t,9»/  4 

DATA  AMNUS^LUS^POlNU'-'  ,'**  ,'  .'/  5 

OATA  AK  .AM.AU/'K'  t'M* t'U'/  6 

N  =  Nl  7 

NSIGN  -  0  8 

I  1  ■  -1  9 

IX  *    0  10 

ISET  «  0  II 

IF  (A(N).EQ.PIUS)  N=N*1  12 

IF  (A(N).NE.AMNUS)  GO  TO  1                        .  13 

NSIGN  =  1  14 

N  «  N*l  15 

16 

1  DO  6  l*N,80  17 
IF  IA( I  I. NE. POINT)  GO  TO  2  18 
ISET  «  1  19 
GO  TC  6  20 

2  IF  IISET.E0.1)  II  -  11»1  21 

DO  3  K= 1, 10  23 

IF  (Ml  J.EQ.BIKI)  GO  TO  4  24 

3  CONTINUE  25 

26 

GO  TO  7  2  7 

28 

*  DO  5  K=l,10  29 

KK  =  K-l  30 

IF  I  A(  I  )  .EQ.B(K)  I  NUMB'KK  31 

5  CONTINUE  32 

33 

IX  =  NUMB»10*IX  3* 
N2  =  1*1  35 

6  CONTINUE  36 

37 

7  IF  (NSIGN. EQ.l )  IX  -  -IX  38 
V  »  IX  39 
IF  ( I  I. IT  .01  II  -  0  40 

X  «  V/(  10** I  I  I  41 
IF  U(N2).EQ. POINT)  N2-N2M  42 
IF  A (N2I.E0.AK)  X  -  X*1000.  43 
IF  ( A(N2) .EO.AM)  X  *  X*0.001  44 
IF  (A(N2) .EO.AU)  X  -  X*0. 000001  45 
IF( (A(N2).EQ.AK) .OR.t AIN2) .EO.AM) .OR.IAIN2) -EO.AU) )  N2-N2U  46 
Nl  ■  N2  47 
RETURN  48 

ENO  49 
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POLPRT 

PURPOSE:       to  control  the  plotting  of  the  polar  plot. 

METHOD:        This  subroutine  is  the  main  subroutine  in  the 
polar  plot  package  and  is  responsible  for  calling  the 
various  subroutines  of  the  package. 

The  scale  factor,  S,  must  be  changed  according  to  the 
printer  characteristics.  The  scale  factor  in  this  subroutine 
is  set  for  ten,  10,  characters  per  inch  for  the  abscissa  and 
eight,  8,  characters  per  inch  for  the  ordinate  axis. 
Therefore  S  =  10. /8. 

After  initializing  DATAX,  DATAY,  and  X,  the  input  data, 
Y,  is  scanned  to  determine  the  normalizing  factor.   If  this 
normalizing  factor  is  less  than  1.E-32,  an  error  statement 
is  printed  and  the  plotting  is  abortted. 

In  the   DO  LOOP   ending  with  statement  8,  each  line  of 
the  polar  plot  is  printed  after  a  call  is  made  to  PTPLOT  to 
establish  the  ploar  grid  information.   The  variable,  DIM,  is 
used  to  as  a  scaling  factor  for  the  polar  plot.  The  value  of 
1.0  will  cause  all  of  input  data  to  be  plotted,  however,  if 
only  the  values  less  than  one-half  of  the  normalizing  factor 
are  of  interest,  then  DIM  can  be  set  to  .5.  This  will 
enlarge  of  the  center  of  the  polar  plot. 

CALLED  BY:      MAIN 


CALLS  TO:       PTPLOT 

SART 
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c 


c 


c 


c 


00  5  1AM.N 


BIN  -  DIM/80.0 
ULL 
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SUBROUTINE  POLPRT  (NAME.YI  I 

COMMON  ISYM.LINE  2 

OIMENSION  Xi360),  Y(360),  DATAXO60),  DATAYI360),  LINEU30),  ISYMI  3 

114)  4 

OIMENSION    TITLAI2),    TITL2(2)  5 

DATA    TITLA/'PHI     »,'THET,/k  6 

N    ■    360  7 

DIM    •    1 .0  8 

NST    -    I  9 

KST    -    1  10 

C  11 

C               S    IS    SCALE    FACTOR    OF    PRINTER!  12 

C               ABSCISSA    CHAR.    PER    INCH    /    ORDINATE    CHAR.    PER    INCH  13 

C  I* 

S    ■    10.0/8.0  15 

C      ZERO  OATAX  AND  OATAY  17 

C  18 

C  19 

00  1  IA*1,N  20 

0  =  IA-1  21 

DATA  X(  IAI  =    0.0  22 

DATA  Y(  1A>  =  0.0  23 

1  X(IA)  =  0*3.1415927/180.0  2* 
C  25 
C  26 
C  FACTOR  IS  THE  NORMALIZING  DIVISOR  27 
C  28 

FACTOR  =  Y( 11  29 

C  30 

DO  2  IA=2,N  31 

2  IF  IFACTOR.LT. YI  I  A )  I  FACTOR-YIIA)  32 
C  33 
C  3* 

IF  INAME.EO.ll  TITL1«TITLA( 1)  35 

IF  INAME.EQ.2)  T I TL1  =  T I TL A ( 2  I  36 

IF  ((NAME. EQ. 3). OR. (NAME. EQ. 4). OR. (NAME. EQ. 71. OR. (NAME. EQ. 81)  TITL  37 

12(1)'TI TLA( 1 )  38 

IF  ((NAME. EQ. 5). OR. (NAME. EQ. 6). OR. (NAME. EQ. 9). OR. (NAME. EQ. 10) )  TIT  39 

1L2(1)=TITLA(2)  40 

IF  t(NANE.EQ.3).0R.(NAME.EQ.51.0R.(NAME.E0.7).0R.(NAME.E0.9 ))  TITL  41 

12I2I-TITLAUI  42 

IF  ((NAME. EQ. 4). OR. (NAME. EQ. 6) .OR. (NAME. EQ. 8). OR. (NAME. EQ. 10)1  TIT  43 

1L2I2)-TITLA(2)  44 

IF  (FACTOR. GT.l.E-32)  GO  TO  3  45 

IF  (NAME. LE. 21  WRITE  (6,91  TITL1  46 

IF  (NAME.GE.3)  WRITE  (6,10)  TITL2  47 

RETURN  48 


49 


C      NORMALIZE  DATA  TO  ONE  50 

r  3 1 

C  5Z 

3  DO  4  IA-l.N  H 

4  YUA)  «  Y(lA)/FACTOR  54 
f  3 ' 
C  56 

IF  (NAME.LE.2)  WRITE  (6,111  TITL1, FACTOR  C1M..                57 

IF  l(NAME.GE.3).AN0.(NAME.LE.6»)  WRITE  (6,13)  TITL2, FACTOR  58 

IF  (NAME.GE.7)  WRITE  (6,12)  TITL2, FACTOR  59 

C      FILL  OATAX  AND  DATAY  ARRAY  FROM  X  AND  Y  ARRAY 


60 
61 
62 
63 


OATA  X(IA)  =  Y(IA)*COS(X( IA) )  ** 

5  OATA  Y4IA)«  Y( I  A ) *  SI N ( X( I  A  1  I  £5 

c  br 

C      SORT  DATA  BY  ORDINATE  MAGNITUDE  £8 


69 
70 
71 


CALL  SART  (OATAX , OAT AY ,N ) 

C      OATAX  AND  DATAY  ARE  SORTED  BY  DESENDING  MAGNITUOE  ON  THE  DATAY  VAL       72 
C      SET  UP  FOR  PLOTTING  POLAR  GRIO  WITH  DATA  " 


C  75 

DO  8  IYY=l,8l  J6, 

CALL  PTPLOT  IIYY,S)  7,8 

C      LINE  IS  RETURNED  WITH  POLAR  GRIO  INFORMATION                               80 


81 


C      SET  UP  »Y»  BIN  SIZE  UPPER  AND  LOWER  LIMITS  82 

C      ULL  IS  THE  LOWER  BIN  LIMIT  83 

C      UL  IS  THE  UPPER  BIN  LIMIT  84 


85 
86 


DIM-(2*1YY-1)*BIN  87 


UL  ■  ULL»2*BIN 

C  *° 

C      CYCLE  THROUGH  OATA  TO  FIND  WHICH  ONES  FALL  IN  «Y«  BINS  91 


92 
93 
IF  (NST.GT.N)  GO  TO  7  ,  '* 

DO  6  JJ-NST.N  96 


IF  (OATAY(JJ).LT.ULl)  GO  TO  7  97 

KST  «  JJ  98 

AMAG  «  SQRTlDATAXl JJ)*DATAX( JJ)*OATAY( JJ J*OATAY( JJ) »  99 

C  100 

C      CHECK  THAT  MAGNITUDE  IS  NOT  OVER  DIM  101 

C  102 

IF    ( AHAG.GT.OIM)    GO   TO    6  103 

C  10* 

C              OK    IS    THE    FINAL    LINE    POSITION    FOR    THE    «*•  105 

C  106 

OK    •    OATAXI JJ)*S**0.0/DIM»61.0  107 

IF    (OK. LT. 10.01    GO    TO    6  108 

K    ■    INT  I  OK  I  109 

K    •    IABS(K)  110 

OK    ■    ABS(OK)  1 11 

IF     (COK-KJ.GT.0.5)    K«K*1  112 

IF    (OK. LT. 10.0. OR. OK. GT. 111. 01     GO    TO   6  113 

LINE(K)     •     ISYMI*)  11* 

6  CONTINUE  115 
C  116 

7  CONTINUE  117 
NST    '    KSTM  118 

C  119 

C      PRINT  OUT  ONE  LINE  OF  PLOT  120 

C  121 

WRITE    (6,14)    LINE  122 

8  CONTINUE  123 
C  12* 

RETURN  125 

C  126 

•9    FORMAT     IIOX.IA*,'     COMPONENT    OF    THE    ELECTRIC    FIELD    IS    LESS'/lOX,  127 

1     'THAN     l.E-64,     THEREFORE    THIS    FIELD    WAS    NOT     • /10X , 'PLOT  TED.          EXEC  128 

2UTI0N    rflLL    CONTINUE    AS    NORMAL.'//)  129 

10  FORMAT     UOX.'THE     MAXIMUM    VALUE    OF    THE    BISTATIC    PATTERN    FOR    •/  130 

1  10X, 1A*.«-* , IA*, •     ( INCIDENT-SCATTERED!     IS    LESS    THAN    •/  131 

2  10X,     •     1.E-30.I          POLAR   PLOT    NOT    CALLED.'///)  132 

11  FORMAT    ('l'.lA*,'     ELECTRIC    FIELO    ANTENNA    PATTERN    FOR    SPECIFIED   PLA  133 
lNE.'^X,  'NORMALIZING    FACTOR-    •.E10.5)  13* 

12  FORMAT     J'IBISTATIC    SCATTERING    PATTERN   F OR  •  .  I  A*. • - • , 1  A*, • ( INC IDENT-  135 
1SCATTERED)    POLAR  II  AT  I  ON. • ,9X , • NORMAL  1 Z 1 NG    F AC  TOR- ' , E10. 5)  136 

13  FORMAT    i'lBACKSCATTERING    PATTERN    FOR' .1 A*. ■ -• . 1 A*, ' ( INC IOENT- SCATT  137 
1EREDI     POLARIZATION. «,9X, 'NORMALIZING    FACTOR-' ,E10. 5)  138 

1*    FORMAT     (IX.  130*1  )  139 

ENO  1*0 
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PTPLOT 

PURPOSE:       to  establish  the  grid  information  for  the 
polar  plot. 

METHOD:        In  the  DO  LOOP  ending  at  statement  1  the 
alpha-numeric  characters  are  transferred  to  ISYN  in  order 
to  pass  via  COMMON  to  other  subroutines.  In  the  statements 
following  statement  2,    the  equations  for  the  plotted 
concentric  circles  are  established.  Below  statement  7  the 
grid  marks  on  the  090-270  axis  are  inserted. 

CALLED  BY:      POLPRT 

CALLS  TO:       LINECK 

NUMB 
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SUBROUTINE  PTPLDT  (IVY, SI 

THIS  SUBROUTINE  SETS  UP  POLAR  GRID  INFORMATION 

COMMON  ISYM.LINE 

DIMENSION  LINE!  1301  ,  ISYMI14I,  ISYNI14) 

OATA  ISYN/  IHt.lH. ,1H  ,  IH* , 1 H/, 1  HO. 1H 1 , 1 H2 , 1 H3 , 1H4 , 1 H5, 1 H6 , 1H8, 1H9/ 

INTEGER  Y,YY,W 

SET  UP  ISYM  FROM  I SYN  FOR  COMMON 

00  1  K=1,14 
ISYM(K)  «  ISYNIKl 

1  CONTINUE 

CLEAR  LINE  ANO  SET  TO  BLANK 

00  2  1=1,130 

2  LINE! I)  *  ISYMI3) 

Y  «  41-lYY 

IF  IY.EO.OI  GO  TO  7 

SET  UP  EQUATIONS  FOR  CONCENTRIC  CIRCLES 

YY    »    Y*Y 

Z    «    IYY»2.5/2)*S 

X    =    6l.0»S0«T(2500.0-Z) 

CALL    LINECK    IX, Yl 

IF    IY.GT.32.0R.Y.LT.-32)    GO    TO    3 

X    «    61.0»S0RTI 1600.0-21 

CALL    LINECK    IX, Yl 

3  IF    IY.GT.24.0R.Y.LT.-24)    GO    TO    * 
X    •    61.0*SQRTI900.0-Z) 

CALL  LINECK  IX, Y) 
*  IF  IY.GT.16.0R.Y.LT.-16)  GO  TO  5 
X  •  61 .0»S0KT(400.0-Z) 
CALL  LINECK  (X.YI 

5  IF  IY.GT.B.0R.Y.LT.-8)  GO  TO  6 
X  •  61.0»S0RTI100-Z) 

CALL  LINECK  IX, Y) 

SET  UP  EQUATIONS  FOR  MULTIPLES  OF  30  OEGREES 

6  X  -  61.0»1.732051*Y»S 
CALL  LINECK  (X,Y) 

X  •  61. 0»Y*S/1. 732051 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

1* 

\l 
17 
IB 
19 
20 
21 
22 
23 
2* 
25 
26 
27 
28 
29 
30 
31 

II 

34 
35 
36 
37 
38 
39 
40 

a 

43 

44 
45 
46 
47 

48 


CALL  LINECK  (X,Y) 

PUT  IN  POLAR  PLOT  NUMBER  LABELS 

CALL  NUMB  IYI 
W  »  IABSIY) 

FILL  IN  POLAR  PLOT  AT  000,  090,  180,  ANO  270 


IF  IW.NE, 

LINEI55I 

L1NEI57J 

L INEI59) 

LINEI63) 

LINEI65I 

LINEI67) 

IF  (W.NE, 

LINEI56 J 

L1NEI5B) 

LINEI60) 

LINEI62) 

LINEI64) 

LINE  I  66) 

IF  (W.NE, 

L1NEI57I 

LINEI59) 

LINEI60) 

LINEI62) 

LINEI63) 

LINEI65 ) 

IF  IW.NE, 

LINEI58) 

LINEI60) 

LINEI62) 

LINEI64) 

IF  IW.NE, 

LINE  I  59  J 

LINEI63) 

CONTINUE 

RETURN 

ENO 


40)  GO 
»  1SYM1 

*  ISYMI 

■  ISYMI 

■  ISYMI 
«  ISYMI 
«  ISYMI 
321  GO 
«  ISY1I 
»  ISYMI 
»  ISYMI 

*  ISYMI 

*  ISYMI 
'  ISYMI 
24)  GO 

ISYMI 
ISYMI 
ISYMI 
ISYMI 
ISYMI 
ISYMI 
6)  GO 
ISYMI 
ISYMI 
ISYMI 
ISYMI 
08)  GO 
ISYMI 
ISYMI 


TO  8 

21 

2) 

2) 

2) 

2) 

2) 

TO  9 

2) 

2) 

2) 

2) 

2) 

2) 

TO 

2) 

2) 

2) 

2) 

2) 

21 

TO  11 

II 

II 

TO  12 

II 


10 


49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
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READ 


PURPOSE:       to  interpret  and  translate  the  input  data 
cards. 

METHOD:        The  program  utilizes  free  format  for  the  data 
cards,  that  is,  the  program  uses  character  recognition  to 
determine  which  parameters  are  being  read.    In  the  IF 
statements  containing  A(1),  A  (2),  A  (3),  and  A  (4),  the  first 
four  characters  on  the  data  card  are  compared  to  the  first 
four  letters  of  the  key  words.  This  will  determine  the  type 
of  parameters  that  card  contains.   The  other  IF  statements 
determine  which  parameters  are  being  read. 

Subroutine  BLNK  is  called  to  remove  the  blank  spaces  on 
the  parameter  cards.   Subroutines  EQUAL  and  LEFT  are  called 
to  determine  the  position  of  the  equal  character  and  the 
left  paren,  respectively.  Subroutine  NUMBER  is  called  to 
convert  the  alpha-numeric  characters  to  numbers,  either 
fixed  or  floating  point.   This  numerical  value  is  assigned 
to  the  parameter  just  determined. 

A  detailed  explanation  of  the  data  cards  is  found  in 
appendix  II  titled  "USERS  MANUAL". 

CALLED  BY:      BLNK 

EQUAL 
LEFT 
NUMBER 
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SUBROUTINE  READ  (  I  A , I  6, I BI SC, IC ARO, IGAI N, IGRO, INEAR, INT , I  SCAT . I WR,  1 

11FLAG,KFLAG,KGEN,L0AD,LZD,MSG,NBAP,NBIP,NFFP,NGEN,NM,NP,ABAP,ABAT,  2 

2AFFPlAFFT,ABlP,ABIT,AM,BM,CMM,ER2,ER31£R41FMC!HGT,PHAF.PHAI,PHIF,P  3 

3HU ,  PHSF.PHSI , THAF.THAI.THIF, T H  I  I  , THSF , THS I  ,  S lG2, S I G3, S IG4, TD2, TO  3  4 

4,V0LT,X,XNP,Y,Y*P,Z,ZLLD,ZNP,STEPJ  5 

CCMMON  /A/  A180I  6 

COMPLEX  VOLU  1)  .ZLLDIl)  7 

DIMENSION  1AI1I.  IBIlIt  XII),  Yll).  III!,  KGENI1),  KFLAG(l)  8 

DIMENSION  XNPIl),  YNPli),  ZNPI 1  I  ,  LZDIl)  9 

DATA  AA,A8,AC,AD,AE,AF.AG.AH,AI,AK,AL,AMA,AN,A0,AP,A0,AR,AS,AT,AU,  10 

1AW,AX/'A•,•8•,,C,,•D■,•E•,•F•,•G•,•H•,•I•,,K•,•L•,•M•,,N•,'0•,•P•,  11 

2»0»  .'R«  ,•  S«  t  '  T"  ,  'U'  ,'W'  ,  'XV  12 

DATA  eLANK, COMMA, MI NUS.PLEF T ,PQI  NT , RIGHT, SLANT/'  •  ,  •  ,•  ,  «-• t»  I  • t« . «  13 

1,  •  J  •,'/•/  14 

RAD  '    57. 295779  15 

INT  *  *  16 

I8ISC  =  -1  17 

1GAIN  *  -1                                        ,  18 

INEAR  =  -1                                       '  19 

ISCAT  =  -1  20 

IWR  =  -I  21 

IF  IIFLAG.EQ.6I  GO  TO  2  22 

IF     IMSG.NE.O)    GO    TO    4  23 

1  READ    (5,76,END=72)    A                                                             •  2* 

2  IF    H  A(  1)  .NE.ACI  .0R.IAI2)  .NE.BL  ANK  ).0R.  (  AI 3  ).NE. BL ANK 1 .OR. I  AI4I.NE  25 
l.BLANKI )     GO    TO    3  26 

WRITE    16,741    A  27 

GO    TO    I  28 

3  WRITE  I  6,75)  29 
GO    TO    5                                                            »•  •  30 

4  READ    <5,76.EN0=72)    A  31 

5  ICARD  =  ICAROM  -  32 
WRITE  (6,771  ICARD, A  33 
IE    nMSG.NE.0).AND.(IAIl).EQ.AE).AN0.IAI2).EQ.AN>.AND.I  A(3)  .EQ.AO)  34 

1))    GO    TO    70  35 

IF     HMSG.NE.01.AND.il AI 1  )  .EQ.  AS  I  . AND.I A I 2  )  .EQ.AT ) .AND. I  A I  3) .EQ.AO)  36 

t.AND. I AI4) .EO.API I )    GO    TO    69  37 

IF    II  A I  1) .EQ.ACI .AND. I AI 2  I . EQ. BL ANK ) . AWD. I  A  I  3  I . EQ.8L ANK ) .AND. I A I  4  I  38 

l.EQ. BLANK) I    GO    TO    73  39 

IF    IMSG.GT.O)    GO    TO    4  40 

CALL    BLNK     (AI  41 

N    -    4  42 

43 

INSULATION  44 

45 

IF    HAI  1).NE.AI).0R.IAI2).NE.ANI.0R.IA(3).NE.AS).0R.UI4).NE.AU))  46 

1G0  TO  10  47 

KFLAGI20)  -  1  48 


CALL  LEFT  INI  49 

6.  1£.!i*iN.,;NE;*R' •OR««*<N»l).NE.AA).OR.(AIN*2).NE.AD).OR.IAIN»3).NE  51 

I . AI I  I  GO  10  7  rj 

KFLAGI4)  -  |  53 

CALL  EQUAL  IN)  54 

CALL  NUMBER  IN,N2,X1,IX)  55 

IF  IAIN2).EQ. RIGHT)  GO  TO  4  57 

IF  I AIN2I.NE. SLANT)  GO  TO  71  58 

N  ■  N2U  59 

GO  TO  6  60 

7  IF  (IAIN) .NE.AD)  .  OR.  I  A INH ) .NE. A  I ) .OR. I  A  I N* 2) .NE . AE  I  .OR.  I  A  IN*  31  .NE  62 
1  .AL) I  GO  TO  8  63 

KFLAGI6)  »  1  64 

CALL  EQUAL  IN)  65 

CALL  NUMBER  IN, N2, XI, IX)  66 

ER2  -XI  67 

IF  IAIN2).EQ. RIGHT)  GO  TO  4  68 

IF  IAIN2I.NE. SLANT)  GO  TO  71  69 

N  •  N2fl  70 

GO  TO  6  71 

8  IE    H A<N).NE.AC).OR.IAIN*l).NE.AO).OR.<AIN»2).NE.AN).OR.IAIN»3).NE  73 
1  .AD)  )    GO    TO   9  74 

KFLAGI5)     »    I  75 

CALL     EwUAL     IN)  76 

CALL    NUMBER    IN, N2, XI, IX)  77 

SIG2    'XI  78 

IF    IA(N2).EQ. RIGHT)    GO    TO    4  79 

IF    IAIN2J.NE. SLANT)    GO    TO    71  80 

N    =    N2U  81 

GO   TO    6  82 

9  IF    !|A«N).NE.AL).OR.tAtN*l).NE.AO).OR.«A<N*21.NE.AS».OR.IAIN»3).NE  84 
l.AS))    GO    TO    71  85 

KFLAGI7)     -    1                   .  86 

CALL    EQUAL     IN)  87 

CALL    NUMBER    IN, N2, XI, IX)  88 

T02   «    X 1                                                                                        '•  89 

IF    (AIN2I.EQ. RIGHT)    GO    TO    4  90 

IF    (AIN2I.NE. SLANT)    GO    TO    71  91 

N   •    N2»l  92 

GO    TO   6  93 

C  94 

C              WIRE                                                                   »  95 

C  96 


n 


10  IF  ((  A(  1) .NE.AW) .0R.IA12).NE.AI).0R.IA<3).NE.AR).0R.<A<*).NE.AE|)  97 
1G0  TO  13  98 

CALL  LEFT  INI  99 

C  100 

11  IF  1 (A(N).NE.AA) .OR. ( A (NM ) .NE . AA > .OR.  I  A <N»2 I .NE . AO) .OR. I AIN» 3) .NE  101 
l.AI ) )  GO  TO  12  102 

•   KFLAGI2)  »  I  103 

CALL  EQUAL  (Nl  10* 

CALL  NUMBER  (N,N2,X1,IX)  105 

AM  =  XI  106 

IF  ( AIN2I .E0.RI3HT1  GO  TO  *  107 

IF  ( AIN2I.NE. SLANT  I  GO  TO  71  108 

N  =  N2*l  109 

GO  TO  11  110 

C  111 

12  IF  (I  A(N) ,NE.AC).0R.<AIN»1) .NE.A01.0R.I A  I N»2 I .NE. ANI .OR. < A<N» 3) .NE  112 
l.ADI)  GO  TO  71  113 

KFLAGI3)  =  1                   .  11* 

CALL  ECUAL  <N)  115 

CALL  NUMBER  (N,N2,X1,IX)  116 

CMM  =  XI  117 

IF  (AJN2I.EQ. RIGHT)  GO  TO  *  1  Id 

IF  (AIN2I.NE. SLANT)  GO  TO  71  119 

N  =  N2»l  120 

GO  TO  11  121 

C  122 

C      EXTERNAL  MEDIUM  123 

C  12* 

13  IF  (  ( A(  1  l.NE.AE)  .0R.IA12)  .NE.AX1  .OR.(A(  31.NE.AT1  .OR.(A(*).NE.AE»)  125 
1G0  TO  1 7  126 

KFLAGI8)  =  1  127 

CALL  LEFT  IN)  128 

C  129 

1*  IF  (U(N).NE  .AC)  .0R.(A(N»1>  .NE.A0).0R.(  A  (N*2)  .NE.  AN)  .OR  .  I  A (N»3) .NE  130 

l.AD) )  GO  TO  15  131 

KFLAGI9I  »  I  132 

CALL  ECUAL  (N)  133 

CALL  NUMBER  (N, NT,  XI, IX)  13* 

SIG3  *  XI  135 

IF  U(N2).E0. RIGHT)  GO  TO  *  136 

IF  (AIN2).NE. SLANT)  GO  TO  71  137 

N  »  N2»l  138 

GO  TO  1*  139 

C  1*0 

15  IF  ( (A(N).NE.AD) ,0R.IAIN»1) .N6.AI J.OR.I AIN*2) .NE.AE) .0R.IA(N»3).NE  1*1 

l.ADI  GO  TO  16  1*2 

KFLAGI10)  ■  1  1*3 

CALL  EQUAL  (Nl  1** 


CALL  NUMBER  IN,  N2,  XI,  IX)  1(.s 

FR3eXl 

IF  (A(N2).EQ. RIGHT)  GO  TO  *                                                147 

IF  (A(N2).NE. SLANT)  GO  TO  71                                               }Zn 

N  «  N2*l  \ll 

GO  TO  I*  }*J» 

l6l!AS))AGo'foE71U*OR*U<N*n*NE**0,*OR*<*<N*2,,NE**SKOR,U1N*3,*NE  >"« 

KFLAGI  11)  •  I  \H 

CALL  EOUAL  IN)  iiZ 

CALL  NUMBER  IN, N2, XI, IX)  }  ij? 
T  03  =  X  1 

IF  (A(N2).E0. RIGHT)  GO  TO  *                                                III 

IF  IAIN2I.NE. SLANT)  GO  TO  71                                               ,  e| 

N  «=  N2fl  }?2 

GO  TO  1*  \l\ 

162 


RMAG  ■=  XI 
N  '    N2» 1 


C 

C      LOAD  \tl 

C  I" 

171G0  |0A{^,*Ne'*L,*OR*<*<2,*NE«*°>«OR«l*<3»«NE.AA|.OR.<A<*).NE.AOH  166 

KFLAGU*)  «  1  \fl 

.   G0  T0  I*  ifcg 

1811^0'{*(^^E«AI».0R.(AI2).NE.AMA).0R.(A(3).NE.AP).0R.IA(*).NE.AEI»  170 

,.  KFLAGI2*)  -  1  !U 

19  I  -  I  I? 
CALL  LEFT  (N)  !if 

20  CALL  NUMBER  IN.N2tXl.IXt  ! j% 
IF  IIX.LE.O)  GO  TO  21  7? 
LZDCI)  =  IX  },? 
N  •  N2»l  iU 
CALL  NUMBER  (N,N2,X1,IXI  170 


179 
180 


CALL  NUMBER  (N,N2,X1,IX)  !§i 

ROEG  -  XI  JSS 

RREAL  «  RMAG*COS(RDEG/RAD)  [at 

RIMAG  «  RMAG»SINIROEG/RAD)  fas 

ZLLOII)  »  CMPLXIRREAL, RIMAG)  1B6 

LOAO  "I  ,23 

IF  IAIN2I.E0. RIGHT)  GO  TO  *  taa 

IF  (A(N2).NE. SLANT)  GO  TO  71  «o 

I  -  1*1  Jqo 

N  -  N2»l  if,? 

GO  TO  20  [4J 
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21  KFLAG(24I  -  -1  193 
LOAD  »  -I  19* 
GO  TO  4  195 

C  196 

C       FREQUENCY  197 

C  198 

22  IF  << A(  11 .NE.AF) .OR.( A(2).NE.ARI  .OR. I A(  3  )  .NE. AE ) . OR. ( A ( 4) .NE . AQI )  199 
1GO  TO  23  200 

KFLAG(l)  »  1  ,                                                   201 

CALL  LEFT  <N)  202 

CALL  NUMBER  (N,N2,X1,IX)  203 

FMC  «  XI  20* 

GO  TO  4  205 

C  206 

C      PLOT  207 

C  208 

23  IF  I < A< 1 ) .NE.AP) .OR. ( A<2I .NE.AL I .OR.(A( 3I.NE.A0) .0R.(A(4I .NE.AT) )  209 
1G0  TO  31  210 

KFLAG(22)  '    I  211 

CALL  LEFT  <N>  212 

C  213 

2*  IF  (  (A(NJ  .NE.AF)  .OR. ( A ( N* 1 ) .NE . AA ) .OR. (  A  I N* 2) .NE . AR) .OR. ( A (N* 3) . NE  21* 

l.AFII    GO    TO    25  215 

IGAIN    «     1  216 

NFFP    «     1  217 

GO    TO   27  218 

25  IF     ( (A(N) .NE.AB)  .OR. ( A ( N» I ) .NE . AI ) .OR. (  A ( N» 2) .NE . AS ) .OR. ( A( N» 3) .NE  219 
l.AT))     GO    TO    26  220 

IBISC    =     1  221 

NBIP    -     1  222 

GO   TO    27  223 

26  IF     ((  AINI  .NE.AB) .OR.  I  A  INM  I  .NE.  AA  I  .OR.  (  A(N»2).NE. ACI .OR. I A(N»3I.NE               27* 
l.AKII    GO   TO    71  225 

ISCAT    »     1  226 

NBAP    -    1  227 

C  228 

C  229 

C  230 

27  DO  28  I «N,80  231 
K  •  !♦}  232 
IF    (A(I).EQ. SLANT)    GO    TO    29  233 

28  CONTINUE  234 
C  235 
C  236 
C  237 

GO   TO    71  238 

29  N  •  K  239 
IF    UA(N) .NE.AT) .OR. ( A ( N» 1 ) .NE . AH I .OR. ( A ( N*2) .NE . AE > .OR. I  A IN* 3) .NE  240 


I .ATI)    GO    TO    30  2*1 

CALL    ECUAL    (N)  242 

CALL    NUMBER    ( N ,*J2 , XI , I Xi  243 

IF     INFFP.EO.l)     AFFT'Xl  •.                                                                                       244 

IF    <NBIP.EO.lt    A3  I T  =  X 1  245 

IF    (N8AP.E0.il    ABAT'Xl  246 

IF     ( A(N2)  .EQ. RIGHT!    GO    TO  4                                                                                                                           247 

IF     (A1N21.NE. SLANT)    GO    TO  71                                                                                                                         248 

N    ■    N2»l  2*9 

GO    TO    24  j                                                                                                                   250 

30  IF  I (AINi .NE.AP) ,0R.< A( N* I ) .NE. AH) .OR. (  A ( N» 2)  .NE . Al  )  i  GO  TO  71  251 
CALL  EQUAL  (Nl  252 
CALL  NUMBER  (N,N2.X1,IX)  253 
IF  (NFFP. EO. II  AFFP=Xl  254 
IF  (NBIP.EO.l)  ABIP'Xl  255 
IF  (NBAP. EQ. II  ABAP-Xl  256 
IF  (A(N2I  .EQ. RIGHT)  GO  TO  4  257 
IF  IA(N2).NE. SLANT)  GO  TO  71  258 
N  «  N2»l  259 
GO    TO    24  260 

C  261 

C              OUTPUT  262 

C  263 

31  IF     (( A(  1) .NE.AO)  .0R.( A(2) .NE . AU 1 .OR. » A( 3  I .NE. AT ) . OR. ( A ( 4 ) .NE . API )  264 
1G0    TO    44  265 

KFLAG(22>     »    1  266 

CALL    LEFT    IN)  267 

C  268 

32  IF    I  I Al Nl .NE.ABI .OR.I A(N*1).NE.AI I.OR.I  A(N*2).NE.AS).OR.(AIN»3).NE  269 
l.AT) )    GO    TO    33  270 

KFLAG118)     ■     1  271 

IBISC    =    1  272 

CALL    ECUAL     (N)  273 

CALL    NUMBER    (N,N2,X1,IXI  274 

PHSI     ■    XI  275 

N    •    N2»  1  276 

CALL    NUMBER    (N,N2,Xl,tX)  277 

PHSF    •    XI  278 

N    -    N2»l  279 

CALL    NUMBER    IN, N2, XI, (XI  280 

THSI    ■    XI  281 

N   -    N2»l  282 

CALL    NUMBER    (N.N2.X1.IX)  283 

THSF    •    XI  284 

IF    (A(N2I.EQ. RIGHT)    GO    TO  4                                                                                                                    285 

IF    (AIN2I.NE. SLANT)    GO    TO  71                                                                                                                 286 

N    »    N2» 1  287 

GO    TO    32  288 
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33l!AF)»AGO,TSE32F',OR,U<N*U-NE-A*KOft*,*,N*2,-NE-*R'-OR-«*«NO».NE  290 

KFLAGU6)  ■  I  2"" 

IGA1N  ■  1  292 

CALL  ECUAL  (N)  293 

PHA1  2U£?ER  ,N»*»2»«»»*>  \l\ 

N  =  N2»l  296 

"LL  NUMBER  IN, N2, XI,  IX)  "7 

^H4h  «  XI  ;»? 

N  =  N2M  299 

CALL  NUMBER  IN, N2, XI, IX)  f0,? 

THAI  =  XI  IHi 

N  =  N2»l  302 

CALL  NUMBER  IN, N2, XI, IX)  \°l 

THAF  =  XI  5s5 

IF  (AIN2J.E0. RIGHT)  GO  TO  4  ?2* 

IF  IAIN2I.NE.  SLANT)  GO  TO  71  I?,6. 

N  =  N2*l  307 

c       G0  T0  32  181 

3*UARt'AGO,foE;oN,*OR'<MN*n-NE-AE)-OR-,*<N*2,-NE-**»'OR-<*«N»3»-NE  311 

KFLAGI19)  •  1  312 

I  NEAR  =  2  3J3 

CALL  EQUAL  IN)  3  4 

{NEAR(='iEQ'PtEm  G°  T°  "  316 

I  =  1  317 

c       G0  T0  36  iiS 

C  320 

C  321 

35  00  37  L-1,50  322 
I  «  L  323 
N    =    NM  324 

36  CALL  NUMBER  IN, N2, XI, IX)  If* 
XNPI I  I  •  XI  326 
N  -  N2M  327 
tt^NUMBER  IN,  N2,  XI,  IX)  f§» 
N  -  N2»l  330 
CALL  NUMBER  IN, N2, XI, IX)  ||l 
ZNPII )  «  XI  |32 
IF  IINEAR.EO.I)  GO  TO  39  HI 
INEAR  «  L*l  33* 
IF    IAIN2I.EQ. RIGHT)    GO    TO    38  \\\ 


N  «  N2  ,'                          337 

37  CONTINUE  338 
C  339 
C  •  .  340 
C  341 

GO    TO    7  1  ,                                                                          342 

38  N2  •  N2»l  343 
INEAR    =     INEAR-1  344 

39  IF  IAIN2I.EQ. RIGHT)  GO  TO  4  345 
IF  IAIN21.NE. SLANT)  GO  TO  71  346 
N  '  N2» 1  347 
GO    TC    32  ...                                                                                                          348 

C  '                                            349 

40  IF  HA<N>.NE.AB).OR.IAINMI.NE.AA).OR.tAtN*2).NE.AC).OR.tAIN»3).NE  350 
l.AKH  GO  TO  41  351 

KFLAGI17)  «  1  352 

1SCAT  =  I  353 

CALL  EQUAL  IN)  354 

CALL  NUMBER  IN,N2,X1,IX)  355 

PHII  «  XI  356 

N  *  N2»l  357 

CALL  NUMBER  IN, N2, XI, IX)  358 

PHIF  *    XI  359 

N  =  N2*l  360 

CALL  NUMBER  IN, N2, XI, IX)  361 

THII  «  XI  362 

N  -  N2U  363 

CALL  NUMBER  IN, N2,  XI, IX)  364 

THIF  '    XI  365 

IF  IAIN2)  .EQ. RIGHT!  GO  TO  4                                                366 

IF  (AIN2).NE. SLANT)  GO  TO  71                                               367 

N  -  N2« 1  368 

GO  TO  32  369 

C  3  70 

41  IF  I<AIN).NE.AC>.0R.<AIN»1I.NE.AU).0R.IAIN»2I.NE.AR).0R.1AIN*3>.NE  371 
l.ARJ)  GO  TO  43  372 

KFLAGI 15»  •  1  373 

IWR  -  1  3  74 

C  375 

C  3  76 

C  377 

00  42  K.-N.80  378 

IF  IAIKI.EQ. RIGHT)  GO  TO  4  379 

N  -  KM  380 

IF  IAIKI.EQ. SLANT)  GO  TO  32  381 

42  CONTINUE  382 
C  383 

GO  TO  71  384 
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C  385 

43  IF  ((A(N)  .NE.AS)  .OR. (AINU)  .NE.AT).OR.(  A < N* 2) .NE . AE I  .OR. ( AIN* 31 .NE  386 
LAP))  GO  TO  71  387 

CALL  EQUAL  (N)  388 

CALL  NUMBER  (N,N2,X1,IX)  389 

STEP  '     XI  390 

IF  (AIN2I.EQ. RIGHT)  GO  TO  4  391 

IF  (At N2) .NE. SLANT)  GO  TO  71  392 

N  =  N2M  393 

GO  TO  32  394 

C  395 

C      FEEO  POINT  396 

C  397 

44  IF    (< A(  1) .NE.AF) .0R.( A(2I .NE.AEI .OR. ( A( 3 ) .NE . AE ) .OR. ( A ( 4  I .NE.AOI )  398 
1G0    TO    45  399 

KFLAGU3I     »    1  400 

GO    TO    46  401 

45  IF    ( I Al 1 ) .NE.AGI .OR.I A(2I .NE.AEI .OR.I Al 31. NE. AN) .OR. IAI4I.NE .AE) )  402 
1G0    TO    49  403 

KFLAGI23)    =    1  404 

46  NGEN  »  0  405 
CALL    LEFT     (N)  406 

47  CALL  NUMBER  (N,N2,X1,IXJ  407 
NGEN  =  NGEN»1  408 
KGEN(NGEN)  =  IX  409 
IF  (A(N2).EQ.R13HT)  GO  TO  4  410 
N  =  N2»l  411 
CALL  NUMBER  (N,N2,X1,IX>  412 
VMAG  '  XI  413 
N  =  N2»l  414 
CALL  NUMBER  IN, N2, XI, IX)  415 
VDEG  =  XI  416 
VREAL  »  VMAG*COS(VDEG/RAO)  417 
VIMAG  =  VMAG*SIN( VOEG/RAO)  418 
VOLT(NGEN)  =  CMPLX ( VRE AL , V  I  MAG)  419 
IF  IAIN2) .EO. RIGHT)  GO  TO  4  420 
IF  (A(N2).NE. SLANT)  GO  TO  71  421 
IF  <<A!N2I. EO. SLANT)  .AND. ( A IN2»  I  I  .EQ. BLANK)  )  GO  TO  46  422 
N  '  N2*l  423 
GO    TO    47  424 

48  READ  15,76)  A  425 
ICARO  «  ICARDH  426 
WRITE  (6,77)  ICARO, A  427 
N  «  1  428 
CALL  BLNK  (A)  429 
GO    TO    47  430 

C  431 

C  432 


C      DESCRIPTION  433 

C  434 

49  IF  (  (  Al  1)  .NE  .AD) .0R.IAI2) .NE.AEI  .OR. ( A{ 3) .NE.AS)  .0R.IAC4)  -NE.ACI )  435 
1G0  TO  52  436 

KFLAGI12I  «  1  437 

J  =  0  438 

CALL  LEFT  IN)  439 

50  CALL  NUMBER  (N,N2,X1,1X)  440 
J  =  J»l  441 
NM  =  J  442 
1A(J)  «  IX  443 
N  =  N2*l  444 
CALL  NUMBER  (N,N2,X1,IX)  445 
1BU)  =  IX  446 

.   IF  (AIN2) .EO. RIGHT)  GO  TO  4                                                  447 

IF  IAIN2) .NE. SLANT)  GO  TO  71                                                 448 

IF  (IA(N2) .EQ. SLANT). ANO.  (  A  I N* 1  J . EQ. BLANK) )  GOTO  51  449 

N  -  N2»l  450 

GO  TO  50  451 

51  READ  (5,76)  A  452 
ICARD  ■  ICAROM  453 
CALL  BLNK  (A)  454 
WRITE  (6,77)  ICARO, A  455 
N  =  1  456 
GO  TO  50  457 

C  458 

C      GECMETRY  459 

C  460 

52  IF  I  I  Al  I) .NE.AG)  .0R.( A( 2) .NE.AE) .OR.(A( 3) .NE.AO)  .OR. ( A ( 4  I .NE . AMA I  )  461 
1  GO  TO  55  462 

KFLAGU2)  '  1  463 

JJ  *    0  464 

CALL  LEFT  (N)  465 

53  CALL  NUMBER  IN, N2, XI, IX)  466 
JJ  *  JJ»1  467 
NP  *  JJ  468 
X(JJ)  «  XI  469 
N  =  N2*l  470 
CALL  NUMBER  (N,N2,X1,IX)  471 
Y(JJ)  =  XI  472 
N  »  N2*l  473 
CALL  NUMBER  (N,N2,X1,IX)  474 
2(JJ)  •  XI  475 
IF  UIN2I .EO. RIGHT)  GO  TO  4  476 
IF  (AIN2I.NE. SLANT)  GO  TO  71  477 
IF  (<A(N2) .EQ. SLANT). AND. (A(N»1I .EO. BLANK))  GO  TO  54  478 

.   N  -  N2»l  479 

GO  TO  53  480 
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5*    READ    (5,76)     A  481 

ICARO    *     ICARD*1  482 

WRITE     (6,77)     ICARO, A  483 

CALL    BLNK    (A)  484 

N    *=    1  485 

GO    TO    53  486 

C  487 

C               INTERVAL    FOR    CALCULATION  488 

C  489 

55  IF  (  I  A<  1)  .NE..AI  )  .OR.  I  A(  2)  .NE.AN)  .OR.(A(  3)  .NE.AT)  .0R.(A(4)  .NE  .AE)  )  490 
1G0  TO  56  491 

KFLAGI21)  *    1  492 

CALL  LEFT  (Nl  493 

CALL  NUMBER  (N,N2,X1,IX)  494 

INT  =  IX  495 

IF  ( A(N2).EQ. RIGHT!  GO  tO  4                                                496 

GO  TO  71  497 

C  498 

C      GROUNO  499 

C  500 

56  IF  (I Al  I)  .NE.AG)  .0R.( A(2)  ,NE.AR».OR.(A(  3I.NE.A0)  .0R.(A(4).NE.AU)  )        501 
1G0  TO  66  502 

KFLAG(25)  =  1  503 

KFLAGI26)  ■=  1  504 

1GRD  =  2  505 

CALL  LEFT  (N)  506 

57  IF  l(AIN).NE.AP)  .OR.  (  A  (N«-l !  .NE.  AE  1  .OR.  (  A ( N*2) .NE . ARI .OR. ( A(N» 31 .NE  507 
■  l.AF))  GO  TO  58  508 

IGPD  '     1  509 

GO  TO  64  510 

58  IF  <  (A(N) .NE.AS) ,0R.( A (N»l 1 .NE . A01 .OR. (  A(N*2) .NE. A01 .OR . ( A(N» 31 .NE  511 
l.AO))  GO  TO  59  512 

ER4  =  30.  513 

SIG4  =  .02  514 

GO  TO  64  515 

59  IF  ((A(N» .NE.AP1  .OR. ( A ( N» 1 I .NE. AOI .OR. ( A (N»2I .NE . AO! .OR. ( A(N»3I .NE  516 
l.ARJI  GO  TO  60  517 

ER4    =    4.  518 

SIG4    =     .001  519 

GO   TO    64  520 

60  IF  ((A(NI.NE.AS1  .0R.(A(N»11  .NE.AE1.0R.( AIN»21.NE.AA!1  GO  TO  61  521 
ER4  »  80.  522 
SIG4  «  4.  523 
GO    TO    64  524 

61  IF     (  ( Al N) .NE.AH)  .0R.(A(N»1)  .NE.AE I.OR.I  A(N*2I .NE.Al I.OR.I AlNOI.NE  525 
1.AGH    GO    TO   62  526 

CALL    EOUAL    IN)  527 

CALL    NUMBER    IN, N2, XI, IX)  528 

HGT    «=    XI  529 

IF    ( AIN2I.EQ.R IGHTI    GO    TO  4                                                                                                                       530 

IF     ( A(N2) .NE. SLANT)    GO    TO  71                                                                                                                     531 

N    «=    N2»l  532 

GO    TO    57  533 

62  IF     ( ( A(N) .NE.AC) .OR.(A(NM) .NE.AOI.OR.I A < N» 2) .NE . AN ) .OR . ( A  IN* 31 .NE  534 
1.A0I)     GO    TO   63  535 

CALL    EQUAL    (N)  536 

CALL    NUMBER    (N,N2,X1,IX)  537 

SIG4    =    XI  538 

IF    (A(N2) .EQ. RIGHT!    GO    TO  4                                                                                                                        539 

IF    (A(N2).NE. SLANT)    GO    TO  71                                                                                                              540 

N    «    N2*l  541 

GO    TO    57  542 

63  IF     ((  A(N-)  .NE. A3)  .OR. IA(N»1). NE.AI). Oft. (A(N*2». NE.AE). OR. «A(N»3I.NE  543 
l.ALU    GO    TO    71  544 

CALL    EQUAL    (N)  545 

CALL    NUMBER     IN,  N2,  XI, IX)  546 

ER4    =    X 1  547 

IF  (A(N2) .EQ. RIGHT)  GO  TO  4                                                548 

IF  (A(N2) .NE. SLANT!  GO  TO  71                                            549 

N  '    N2»l  550 

GO  TO  57  551 

C  552 

C  553 

C  554 

64  00  65  K=N,80  555 
IF  (A(K).EO. RIGHT!  GO  TO  4  556 
N  =  K»l  557 
IF  (AtKI.EQ. SLANT!  GO  TO  57  558 

65  CONTINUE  559 
C  560 
C  561 
C  562 

GO  TO  71  563 

C  564 

C  565 

66  IF  (( A(  ll.NE.ASI .0R.( A(2». NE.AT! .OR. ( A( 3  I .NE. AO! . OR. ( A  I  4  I ,NE . API  1  566 
1G0  TO  67  567 

IFLAG  «  2  568 

RETURN  569 

C  5  70 

67  IF  <( A( lt.NE.ACI .0R.( A(2). NE.AH) .OR.(A( 3 ) .NE . AA ) . OR. ( A ( 4 ) .NE . AN) )  571 
1G0  TO  68  572 

IFLAG  -  3  573 

RETURN  574 


68  IF  l(A(  1). NE.AE). OR. ( A(2). NE.AN)  ,OR.(A(  3I.NE.A0)  )  GOTO  71 


1} 


IFLAG  =  1  578 

RETURN  579 

69  IFLAG  *  5  580 
RETURN  581 

70  IFLAG  *  4  582 
RETURN  583 

71  MSG  *  1  584 
KFLAGI30I  =  ICARO  585 

72  IF  UFLAG.NE.5)  WRITE  (6,781  [jf6, 
IFLAG  '  5  588 
RETURN  589 

C  590 

73  IFLAG  =  6  591 
ICARO  •  ICARO-l  592 
RETURN  59} 

C  594 

C  595 

°   74  FORMAT  (5X.80A1)  „.-_-,,,.                                             597 

75  FORMAT  (////5X,«0ATA  CARDS*//)  598 

76  FORMAT  (80A1)  599 

77  FORMAT  <6X  ,  I  2  ,2X  .  80A1  J  ,..„_,. CTnD  r.on  „ ,  c<  |Nf-««*»i  1  600 

78  FORMAT  (•  i**»t  END  CARO/STOP  CARD  MISSING**"  I  got 
END  79 


RITE 


PURPOSE:       to  generate  a  list  of  branch  currents  from 
the  input  loop  currents. 

METHOD:        The  generation  of  branch  currents  is 
accomplished  in  the  DO  LOOP  ending  at  statement  2.  The 
branch  currents  are  stored  in  CJ  (I)  by  the   latter  part  of 
the  DO  LOOP  ending  at  statement  3.  If  the  branch  currents 
are  reguested  for  output  (IWR  positive) ,    the  DO  LOOP  ending 
at  statement  5  accomplishes  this. 

CALLED  BY:      GANT1 

GFFLD 

CALLS  TO:       NONE 
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SUBROUTINE  RITE  (  I  A  ,  I B  ,  INM,  IWR,  I  I , I  2 , 1 3 , MO. ND.NM ,C J ,CG,  1GRO)  1 

COMPLEX  CJIll  .CG(1  ), CJA, CJB  2 

.   DIMENSION  14(1),  16(11.  Ulll,  12(1).  13(1).  M0UNM.4),  NO(l)             3 

AMAX  »  .0  4 

C  5 

C  6 

DO  3  K=l ,NM  7 

KA  *  IA(K>  8 

KB  «  IB(K)  9 

CJA  ■  I .0,-0)  10 

CJB  ■  ( .0..0)  11 

NDK  ■  NOIK)  12 

C  13 

C  14 

00  2  1 1 -I .NDK  15 

1  '  MD(K, tl)  16 
FI  «  1.  17 
IF  (KB. EQ.I2I  I))  GO  TO  1  18 
IF  (KB. EO, 11( I) )  FI  — 1.  19 
CJA  =  CJA*FI*CJ(I  )  20 
GO  TO  2  .                                            21 

1  IF  (KA.E0.I3(  I))  FI  — I.  22 
CJB  '    CJB*FI»CJ« I)  23 

2  CONTINUE  2* 
C  25 
C  26 

CG(K)  =  CJA  27 

KK  =  K*NM  28 

CG(KK)  =  CJB  29 

ACJ  «=  CABS(CJA)  30 

BCJ  =  CABS(CJB)  31 

IF  (ACJ.  GT. AMAX)  AMAX-ACJ                                                       32 

IF  (BCJ. GT. AMAX)  AMAX-BCJ                                                       33 

3  CONTINUE  34 
C  35 
C  36 

IF  (IWR.GT.O)  GO  TO  4                                                        37 

RETURN  38 

4  IF  (AMAX.LE.O.)  AMAX-1.  39 
WRITE  (6.8)  40 
NMG  ■  NM  41 
IF  (IGRD.GT.O)  NMG  •  NM/2  42 

C  43 

00  5  K=1,NMG  44 

CJA  -  CG(K)  45 

KK  ■  K*NM  46 

CJB  -  CG(KK)  47 

CCJA  -  CABS(CJA)  48 


CCJB  ■  CABS(CJB)  49 

ACJ  =  CCJA/AMAX  50 

BCJ  *  CCJB/AMAX  51 

PA  «=  .0  52 

PB  =  .0  53 

IF  (ACJ.GT.O.)  PA  »  57.29578*ATAN2(AIMAG(CJA) .REAL1CJA) )  54 

IF  (BCJ.GT.O.)  PB  ■  57.29578*ATAN2(AIMAG(CJB),REAL(CJB) )  55 

5  WRITE  (6,7)  K.1A(K). CJA, CCJA. ACJ.PA.IB(K), CJB. CCJB. BCJ.PB  56 

C  58 

WRITE  (6,6)  59 

RETURN  60 

C  61 

C  62 

6  FORMAT  I1H0)  63 

7  FORMAT  (2X,I2,2(2X,[2,2X,E11.5,1X,E11.5,  U.E11.5,  1X.E11.5, 1X.F6.1)  64 
1)  65 

8  FORMAT  (/2(46X, "NORMALIZED", 5X)/«  SEG',2("  NODE" , 4X, "RE AL « ,6X , • I  MA  66 
1GINARY" .SX.'MACNITUDE'taXt'MAGNITUOE ',3X1 'PHASE') )  67 

END  66 
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SART 


PURPOSE:       to  sort  data  for  polar  plot. 

METHOD:        This  subroutine  sorts  the  values  of  the 

points  to  be  plotted  by  the  polar  plot  package  starting  with 

the  greatest  positive  value  of  y  to  the  greatest  negative 

value.   In  the  DO  LOOP  ending  at  statement  1,  the  value  of 

(x  ,y  )  is  interchanged  with  the  value  of  (x  ,y  )  if  y   is 
i   i  3      J      J 

greater  than  y  . 
i 

CALLED  BY:      POLPRT 


CALLS  TO:       NONE 


SUBROUTINE  SART  < OAT AX .DATAY, N)  1 

DIMENSION  OATAXI500I,  0ATAVI500)  2 

C  3 

C      THIS  ROUTINE  SORTS  DATA  IN  DATAY  BY  MAGNITUOE                              4 

C  5 

NN  ■  N-l  6 

C  7 

00  2  1=1, NN  8 

NM  •  1*1  9 

C  10 

00  1  J=NM.N  11 

IF  IDATAyII) .GE.OATAYI Jll  CO  TO  1  12 

STOR  ■  DATAY! I »  13 

DATA  Y(  I)  -  OATAYU)  14 

DATA  Y( Jl  «  STOR  15 

STOR  -  DATAXI I)  16 

DATA  XIII-  DATAX(J)  17 

OATA  X( J)  -  STOR  18 

1  CONTINUE  19 
C  20 

2  CONTINUE  21 
C  22 

RETURN  23 

END  24 
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SGANT 


PURPOSE:       to   calculate  the  mutual  impedance  between 
filamentary  monopoles. 


SOURCE 
DIPOLE 


RECEIVING 
DIPOLE 


Z  =  P11  ♦  P12  +  P21  +  P22 


METHOD:         In  the  induced  emf  formulation,  the  mutual 
impedance  of  coupled  dipoles  is 


Z  =  - 


/ 


I  (t)  E  (t)  dt 

2         1    ' 


where  I  (t)  denotes  the  current  distribution  (normalized  to 
2 

unit  terminal  current)  on  dipole  2,  and  E  (t)  is  the  field 

1 

of  dipole  1  when  it  transmits  with  unit  terminal  current. 

Distance  along  the  axis  of  dipole  2  is  denoted  by  the 

coordinate  t.  E   may  be  expressed  as  the  sum  of  the  fields 

1 

from  each  of  the  monopoles  comprising  dipole  1.  Furthermore, 

the  integral  is  the  sum  of  the  integrations  over  each  of  the 

monopoles  comprising  dipole  2.   Thus,  the  dipole-dipole 

mutual  impedance  may  be  expressed  as  the  sum  of  four 

monopole-monopole  impedances. 

It  may  be  convenient  to  draw  the  above  figure  in  terms 
of  monopoles  with  the  current  distribution  shown  as  dotted 
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lines.   (The  monopole  letters  remain  the  same.) 


P11 


P12  i 


\ 


P21 


B 


P22  * 


The  surface  impedance  is  calculated  just  above  statement 

2.  B01  denotes  J   /J   where  J   and  J   are  the  Bessel 

0     1        0      1 

functions  of  order  zero  and  one  with  complex  argument,  ZABG. 

It  is  assumed  that  all  the  wire  segments  have  the  same 

radius,  conductivity  and  surface  impedance. 


In  the  DO  LOOP  ending  with  statement  3,  SGANT  calculates 
the  segment  lengths  D ( J) .  DMIN  and  DMAX  denote  the  lengths 
of  the  shortest  and  longest  segments.  If  the  wire  radius  or 
the  segment  lengths  are  clearly  beyond  the  range  of 
thin-wire  theory,  N  is  set  to  zero  at  statement  4  followed 
by  RETURN  to  the  main  program  to  abort  the  calculation. 

At  statement  5,  the  program  selects  a  segment  K,  and  a 
few  statements  below  this  it  selects  another  segment  L.  K  is 
a  segment  of  test  diple  I,  and  L  is  a  segment  of  expansion 
mode  J.  The  mutual  impedance  between  segments  K  and  L  is 
obtained  by  calling  subroutine  GGS  or  GGMM.  In  statement  18, 
this  impedance  is  lumped  into  C  (MMM) .  The  mutual  impedance 

Z    between  dipoles  I  and  J  is  the  sum  of  four 
ij 

segment-segment  impedances. 

The  variables  IFLAG  and  JFLAG   are  used  if  a  ground 
plane  is  present  for  the  calculation  of  the  mutual  impedance 
elements.   If  IFLAG  is  equal  to  JFLAG,  the  mutual  impedance 
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terms  will  not  have  the  effects  of  a  ground  plane  since  both 
monopoles  lie  on  the  same  side  of  the  ground  interface.  If 
the  monopoles  are  on  the  opposite  sides  of  the  interface 
(IFLAG  not  egual  to  JFLAG) ,  the  reflection  coefficient 
correction  must  be  applied  to  the  mutual  impedance  elements. 
This  same  technigue  is  applied  in  subroutines  GNFLD  and 
GFFLD. 

In  SGANT,  segment  K  has  endpoints  KA  and  KB,  and  segment 

L  has  endpoints  LA  and  LB.  It  is  convenient  to  think  of  KA 

and  KB  as  points  1  and  2  on  segment  K,  and  LA  and  LB  as 

points  1  and  2  on  L.   The  four  segment-segment  impedances 

can  be  defined  as  P(IS,JS).   The  first  subscript  IS  refers 

to  the  terminal  point  on  segment  K,  and  the  second  subscript 

JS  refers  to  the  terminal  point  on  L.  Thus  IS=1  or  2  if 

dipole  I  has  its  terminal  point  12  (I)  at  KA  (point  1)  or  KB 

(point  2) ,  respectively.  Similarly,  JS=1  or  2  if  mode  J  has 

its  terminal  point  12  (J)  at  LA  or  LB.  The  impedances 

P(IS,JS)  are  defined  with  the  following  reference  directions 

for  current  flow:  from  point  1  toward  point  2  on  each 

segment.  If  dipole  I  has  this  same  reference  direction  on 

segment  K,  FI=1;  otherwise  FI=-1.  Similarly  FJ=1  or  -1  in 

accordance  with  the  reference  direction  for  mode  J  on 

segment  L.  In  statement  18,  P(IS,JS)  is  multiplied  by  FI  and 

FJ  before  its  contribution  is  added  to  Z 

ij 

Subroutine  GGMM  calculates  the  impedances  Q(KK,LL)  which 
are  like  the  P(IS,JS)  but  have  different  conventions  for 
reference  directions  and  subscript  meaning.  The 
transformation  from  the  Q  impedances  to  the  P  impedances  is 
accomplished  in  the  DO  LOOP  ending  with  statement  13. 

If  the  wire  has  finite  conductivity,  the  appropriate 
modification  is  applied  to  the  impedance  matrix  just  above 
statement  15.   The  terms  arising  from  the  dielectric  shell 
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on  an  insulated  segment  are  obtained  from  subroutine  DSHELL 
just  above  statement  16.  Finally,  the  lumped  loads,  ZLD,  are 
added  to  the  diagonal  elements  of  the  impedance  matrix  in 
the  DO  LOOP  ending  at  statement  23. 

K  is  a  segment  of  test  dipole  I,  and  L  is  a  segment  of 
expansion  mode  J.  When  the  segment  numbers  K  and  L  are 
egual,  SGANT  calls  GGMM  to  obtain  the  mutual  impedance 
between  two  filamentary  electric  monopoles.  These  monopoles 
are  parallel  and  have  the  same  length.  Monopole  K  is 
positioned  on  the  axis  of  the  wire  segment,  and  monopole  L 
is  on  the  surface  of  the  same  wire  segment.  Thus,  the 
displacement  is  egual  to  the  wire  radius.  The  two  monopoles 
are  side-by-side  with  no  stagger. 

When  segments  K  and  L  intersect,  SGANT  again  calls  GGMM 
for  the  mutual  impedance  between  the  two  filamentary 
monopoles.  Monopole  K  is  situated  on  the  axis  of  wire 
segment  K,  and  monopole  L  is  on  the  surface  of  wire  segment 
L.  The  axes  of  segments  K  and  L  define  a  plane  P,and 
monopole  K  lies  in  this  plane.  Monopole  L  is  parallel  with 
plane  P  and  is  displaced  from  it  by  a  distance  egual  to  the 
wire  radius. 

CALLED  BY:      MAIN 

CALLS  TO:       CBES 

DSHELL 

GGMM 

GGS 
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SUBROUTINE    SGANT     (  I  A  ,  I  6  ,  I  NM  .  IN*T  ,  I  SC  .  I  1  ,  I  2,  I  3,  JA,  JB  ,MD,  N  ,  NO  ,NH  ,NP,  A  1 

1H,BM,C,CG0,CMM, 0,EP2,EP3.ETA,FHZ,GAM,SGO,X,Y, Z  , ZLD , ZS, ERR, 1 &RO)  2 

CCHPLEX    ERR  3 

COMPLEX    ZG.ZH.ZS.EGO.GO.CGDStSGOS.SGDT.BOl  4 

COMPLEX    Pll,Pl2,P21,P22.011,ai2,Q2l,Q22tEP2,EPiETA,GAM,EP3  5 

COMPLEX    EPSILA.CWEA.BETA.ZARG  6 

COMPLEX    P(2,2),QI 2,2>,Cg6(1 ) , SGD 1 1 ) ,C ( 1 ) ,ZLD(1)  7 

DIMENSION    X (  I  1  .     Ylil,     Z(l)t     0(1),     IA(1>.      16(1),     Hp(INM.4|  8 

DIMENSION     11(1),     12(1).     13(1).     JA(l),     J 6(  1  )  ,    NO(l),     IS(!(1)  9 

OATA    E0,TP,U0/8.854E-12,6.283i8,1.2566E-6/  10 

EP    »    EP3  11 

ICC    »    <N*N»N)/2  12 

13 

00    1    1=1, ICC  1* 

1  C(  I)    -    (  .0, .0)  15 

16 

■     ZS    =    I.0..0)  17 

IF     (CMM.LE.O.)     GO    TO    2  IB 

OMEGA    =     TP*FHZ  19 

EPSILA    -    CMPLX(E0,-CMM*1.E6/0MEGA)  20 

CHEA    •     <  .0,1 . »*OHEGA*EPSILA  21 

BETA    •    OMEGA«SQRT(U0)*CSQRTIEPSILA-EP)  22 

ZARG    ■    BETA»AM  23 

CALL    CBES     (ZARG.801 )  24 

ZS    •    B("TA*BOl/C*EA  25 

2  ZH  •  ZS/( TP*AM»GAMI  26 
OMIN  ■  1.E30  27 
OMAX    •     .0  28 

?9 

DO    3    J'l.NM  30 

K    •    IA( j[  31 

L    '    1B(J)  32 

0(J)    -    SQRTUX(K)-X(L))**2»(Y(K)-YCLi)»*2HZtK)-Z(L>)**2)  33 

IF     (D(  JJ.LT.OMIN)     OMIN«D(J)  34 

IF    (D(  JJ.GT.DMAX)     DMAX-DUJ  35 

EGO   =    CEXP(GAM*DI J))  36 

CGD(J)    »    (EGD»l./EGOI/2.  37 

3  SGDU)  •  <EG0-l./EG0)/2.  38 

39 

IF  <0M1N.LT.2.»AM)  GO  TO  4  40 

IF  (CABS(GAM»AM) .GT.0.06)  GO  TO  4  41 

IF  (CABSI GAM»0MAX).GT.3.)  GO  TO  4  42 

IF  (AM.GT.O.)  SO  TO  5  43 

4  CONTINUE  44 
N  =  0  45 
WRITE  (6,241  AM.DMAX.DHIN  46 
WRITE  (6,25)  47 

48 


5  00  19  K*1,NM  **> 
IFLAG  *  0  50 
IF  (( ISRD.GT.O) .AND.IK.GT.NM/2) )  IFLAG-1  51 
NOK  ■  NO(K)  52 
KA  «  IA(Kl  |3 
KB  «  IB(K)  54 
OK  «  01 K)  55 
CGDS  «  CGD(K)  56 
SGDS  ■  SGO(K)  57 

J   b 

00  19  L-l.NM  59 

JFLAG  =0  69 

IF  (( IGRD.GT.O) .AND.IL.GT.NM/2) »  JFLAG-1  61 

NOL  «  NOIL)  6Z, 

LA  «  IA(L)  J| 

LB  -  IB (LI  ** 

DL  *  OIL)  £5 

SGDT  •  SGO(L)  66. 

NIL  «  0  6J 

o  o 

00  19     I  1  =  1. NOK  69 

1  «  HOIK. it)  7° 
MM  -  tI-l)*N-(l*I-l)/2  l\ 
FI  «  1.  H 
IF  (K6.E0.I2I II)  GO  TO  6  J3 
IF  1KB.  EO.IU  I)  )  FI  — 1.  It 
IS  «  1  H 
GO    TO    7  1% 

6  IF  (KA.E0.I3I  I)  )  FI  — 1-  H 
IS  '  2  78 
'     £  79 

7  DO  19  JJ=1,NDL  f° 
J  *  MO(L.JJ)  gl 
MMM  '  MM*J  82 
IF  (I.GT.J)  GO  TO  19  83 
FJ  '  1.  84 
IF  (LB.EQ.I21J))  GO  TO  8  85 
IF  (LB.EQ.IKJ)  )  FJ— 1.  86 

GO  TO  9  88 

8  IF  (LA.EQ.I3UI  I  FJ— 1.  89 
JS  •  2  90 

9  IF  (NIL.NE.O)  GO  TO  18  21 
NIL  •  I  H 
IF  (K.EQ.L)  GO  TO  14  ,. _  _,  21 
IND  -  ILA~KA)«(LB-KA)»UA-KBJ»ILB-KB)  ** 
NGRO  "  IGRO  g* 
IF  (IFLAG. EQ. JFLAG)  IGRO  — 1  9* 
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IF  (INO.EQ.OI  GO 
SEGMENTS  K  AND  L 


TO  10 
_  SHARE  NO  POINTS 

Call  ggs  (x<ka),v(ka»,z(ka).x(k8),y(kbj,z<k8),xua),y(lai,Z(la).x< 

iLB».V<LBi.zlLBi.AM.DK.CGDS.SGDSt6LtSG0T,lNT,ETA,GAM,Pa,n.Ml,Jl, 

2Pl2.it, pi  >.2I,err,1gr6j 

JGRO    ■    NGRD 

SEGMENTS    K   ANQ   L    SHARE   ONE    POINT    (THEY    INTERSECT) 
10  KG  -    0 
JM    -    KB 
JC    -    KA 
KF    -    1 
INO   -    (KB-LA}*<KB-LBI 


IF    (  IND.NE.O)    GO    TO 
JC   -    KB 
KF    -   -I 
JM  -   KA 


11 


KG  -    3 


LG 

JP 

IF 
JP 
LF 

LG 
12    SGN 


3 

-  LA 
■  -1 
(L8.EQ.JC) 

-  LB 

-  1 
•    0 

KF*LF 


GO    TO    12 


CPSI    >     l(X(JP)-X(JC))*(X(JM)-X(JC 
J-ZIJCJH/ 

'.ALL    i»t>HM     I. Of  UK 

11, Q<2, 11,0(2, 21) 


11  JP|-Z(JC»I*<Z(JM 
CALL    GGMM    ' 


IIMY1 


JP)-VIJC))*(V(JN)-V(JC)*)»«Z 


.C&DsTs&OS.SGDT.CPSI.ETA.GAM, 


00    13    KK-1,2 
KP    «    IABS(KK- 


KGI 


00    13    LL-1,2 
LP    «    IABS(LL-LG) 
P(KP.LP)     -    SGN*0«KK,LL) 
13   CONTINUE 


14 


IGRD-NGRO 
GO  TO 
K 


0  TO  18 

"L   (SELF  I 


GO 
ZG  < 

Oil 


REACTION  OF  SEGMENT 

..01 
CMM.IE.O.)  GO  TO  15 

GAM*OK 

ZH/IS6DS»*2) 
•  ZC*(SGDS*CGDS-GD)/2. 


K) 


97 
98 
99 

100 
101 

\°i 

103 

104 

05 

Ob 

07 

_08 

109 

10 

•l* 

12 

3 
4 
15 
16 
17 
8 
19 
20 
21 

?2 
123 

l24» 

125 

126 

127 

128 

29 

30 

31 

32 

133 

\% 
36 
37 

IS 

40 

1 


012  -  ZG*(GD*CG0S-SG0S)/2. 

15  ISCK  •  ISC(K) 

Pil  -  J. 5, .ol 

P12  -  1 .0..GI 

IF  (ISCK.EQ.O)  GO  TO  16 

IF  (BM.LE.AN)  GO  TO  16 

CALL  OSHELL  J  AM,  BM .OK.CGOS, SG0S»EP2,EP.ETA, GAM.Pl 1 , P12) 

16  Oil    -    PI  1*01 1 


012   -_P12*Q1 


;aLL    GGMM    J.O.DK, .0,0K .AM.CGOS, SGOS.SGOS . 1. , ETA, GAM, PI 1 , P12, P21.P2 

Oil   -    PIK-Qll 
012   -    P12»012 

?  1:11  :  8.2 

P(2,l)    -    012 
PI2.2)    «    Gil 


17 


112   - 

\\\\\ 

'(2,11 
(2,21 
F  iKI 
10   TO 

pit:!!  :  :81f 

p(2,n  «  -oil 


•(2.21    »    Gil 

IF    (KA.NE.LA)    GO    TO    17 

GO   TO    18 


P(. 

P(2,2)    -    -012 

18  C(MMN)    -    CIMMH»*FI*FJ*PIIS, JS) 

19  CONTINUE 


00  23    1-1, N 

MM   «    (I-lJ»N-( 1*1-11/2 

1  J   ■    MM* I 
JJA    «    JA(I) 
Jl    -    JJA 
112   >    121  I  I 
III    •    IKI) 

IF    (II2.EQ.IB(J1I)    Jl-JUNM 
JJB  -    JB( II 

H  (II2.EQ.IBI J2) 
C(1J)  •  C(  IJ1*ZLD 
JJJ   -    JJA 


J2»NM 


I    J2-JJ 
<J1)*ZLD(J2) 


00   22  K-1,2 
NOJ  -   NO? JJJ I 

00   21    JJ-l.NOJ 

J  -  moIjjj.jji 

IF    tJ.EQ.li    GO   TO    21 

IF    U2(  JI.NE.II2I    GO   TO  21 


1*5 
146 
147 
148 
49 

I? 
52 
53 
54 
.55 
56 
57 

•  5> 

59 
60 

li 

64 
65 
66 
67 
68 

$9 
70 
71 
72 
73 
74 
75 
.76 
77 
78 
79 
80 


83 
1 84 

85 

§? 

a! 

90 

92 
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IJ  -  MM»J 

FI  «  1. 

IF  (K.EQ.2)  60  TO  20 

IF  (IKJI.NE.II1)  FI  — 1. 

C(IJ)  -  C(IJ)»FI*ZL0U1> 

GO  TO  21 

20  IF  (I3(J).NE.I3<m  FI  — 1. 
CIIJ)  «  C(  IJ»*FI»ZL0( J2I 

21  CONTINUE 

22  JJJ  -  JJB 

23  CONTINUE 
RETURN 

24  FORMAT    (3Xt'AM   -    • tE10,3i3Xt *OMAX   -    • ,E 10. 3.3X. • DMIN   -    •1E10.3> 

2  5   FORMAT    (•    WARNING    ****••*••**»*•*•»•****•****•***»***»•»*»****«••/ 
It*    THIS    PROBLEM   |XCEE0  LIMIT   OF   THIN  MIRE  CONOITIONt    THE   RESULTS 

END 


89 


SORT 


PURPOSE:       to  define  the  set  of  dipole  modes. 

METHOD:        In  the  DO  LOOP  ending  at  statement  3,  the  set 
of  dipoles  is  defined  by  filling  the  vectors  11(1)  and  13(1) 
(the  endpoints  of  dipole  I)  ;  12  (I)  (  the  terminal  point  of 
dipole  I);  and  the  vectors  JA  (I)  and  JB(I)  (the  monopoles 
comprising  dipole  I)  with  the  node  numbers  and  segment 
numbers,  respectively.  The  DO  LOOP  ending  at  statement  8 
determines  MD(J,K)  (the  list  of  dipoles  sharing  segment  J) 
and  ND  (K)  (the  number  of  dipoles  sharing  segment  J). 

CALLED  BY:      MAIN 
CALLS  TO:       NONE 
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00  7  L=li2 
NOU)  ■  NO(J)U 

K    *     1 

M  «  0 

5 

MJK  '     MD(J.K) 

IF  (MJK.NE.O)  SO  TO 

M  *  1 

HD(J.K)  ■  I 

6 

K  ■  K*l 

IF  (K.GT.4)  GO  TO  7 

IF  (M.EQ.O)  GO  TO  5 

c 

7 

J  «  JB(  I) 

8 

CONTINUE 

c- 

MIN  >  100 
MAX  ■  0 

SUBROUTINE    SORT    ( I  A, 1 B , 1 1 , I  2 , I  3 , J A, JB, MD ,ND, NM, NP, N.MAX , Ml N, I C J , IN  1 

1M)  2 

DIMENSION    JSPI20I  3 

DIMENSION     11(1),  12(1).     1311),     JA(1),     JB(l)                                                                                        4 

OIMENSION     lAlli,  lBdi,    NDUI,     M0(lNM,4i                                                                                          5 

I    '    0  6 

C  7 

DO    3    K=1,NP  8 

NJK    «    0  9 

C  10 

DO  I  J=1,NM  11 

IND  ■  (  IM Jl-K)*( IB(J)-K)  12 

IF  I IND.NE.01  GO  TO  1                                                        13 

NJK  »  NJK*l  1* 

JSP(NJK)  -  J  15 

1  CONTINUE  16 
C  17 

MOD  ■=  NJK-1  18 

IF  (MOD.LE.O)  GO  TO  3                                                    19 

C  .                     20 

DO  2  IMDM.MOD  21 

I  =  1*1  22 

IF  (  l.&T.ICJ J  GO  TO  2                                                    23 

I  PD  «  IMDH  2* 

JAI  »  JSP( IMD)  25 

JA(I)  =  JAI  26 

JBI  «  JSP(IPD)  27 

JBIII  ■  JBI  28 

11(1)  •  I  AC  JAI )  29 

IF  (IA(  JAD.EO.KI  Il(l)-IBCJAI)  30 

I2( I )  ■  K  31 

13(1)  '  lA(JBI)  32 

IF  (IAC JBII.EO.K)  I3(I)-IB(JBI>  33 

2  CONTINUE  34 
C  35 

3  CONTINUE  36 
C  37 

N  -  I  38 

C  39 

00  4  J-l.NM  40 

NOCJ)  -  6  41 

C  42 

00  4  K- 1,4  43 

4  MOIJ.K)  •  0  44 
C  45 

III  -  N  46 

IF  (N.GT.ICJ)  III  -  1CJ  47 

C  48 


00  8  1*1,111  50 

J  -  JAC  I)  ft 

52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 

DO  9  J-l.NM  19 

NOJ  -  NOU)  77 

IF  (NDJ.GT.MAXI  MAX-NOJ  it 


9  IF  (NOJ. LT. MIN)  MIN-NOJ  }{ 

RET' 

ENO 


RETURN  }l 
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SQROT 


PURPOSE:       to  solve  the  set  of  simultaneous  equations  to 
determine  the  currents  on  the  thin  wire  structure. 

METHOD:        This  subroutine  considers  the  matrix  equation 
ZI  =  V  which  represents  a  system  of  simultaneous  linear 
equations.  NEQ  denotes  the  number  of  simultaneous  equations 
and  the  size  of  the  matrix  Z. 

On  entry  to  SQROT,  S  is  the  excitation  column  V.  On 
exit,  the  solution  I  is  stored  in  S.   Z(I,J)  denotes  the 
symmetric  square  matrix.   Also  on  entry,  the  upper-right 
trianqular  position  of  Z(I,J)  is  stored  by  rows  in  C(K)  with 

K  =  (I  -  1)  *  NEQ  -  (I  *  I)  /  2  +  J  . 


If  112  =  1,  SQROT  will  transform  the  symmetric  matrix 
into  the  auxiliary  matrix  (implicit  inverse) ,  store  the 
result  in  C  (K)  and  use  the  auxiliary  matrix  to  solve  the 
simultaneous  equations.  If  112  =  2,  this  indicates  that  C(K) 
already  contains  the  auxiliary  matrix. 

The  transfromation  from  the  symmetric  matrix  to  the 
auxiliary  matrix  is  accomplished  in  the  DO  LOOP  endinq  at 
statement  5.   The  solution  of  the  simultaneous  equations  is 
accomplished  in  the  remainder  of  the  program. 

CALLED  BY:      GFFLD 
CALLS  TO:       NONE 


92 


SUBROUTINE    SQROT     I  C  ,  S  ,  I WR , I 12.NEQI  1 

COMPLEX    C(  1J.SU  ),SS  2 

N    «    NEO  -3 

IF    (  I12.EQ.2I    GO    TO   6  4 

Cll)    '    CS0RT(C(  Ml  5 

C  6 

DO    1    K«2,N  7 

1  CIKI    -    CtKl/Clll  8 

9  9 

C  10 

DO    5    1*2. N  u 

IMO   ■    I-l  12 

IPO  •    1M  13 

ID   »    (I-l)»N-II»I-I)/2  1* 

II    '    10*1  15 

C  16 

DO    2    L-l.IMO  [7 

LI    •    <L-ll*N-(L»L-L>/2»l  18 

2  CUM    •    C(  II  l-CUl  I  *C  I  L  I  I  19 

C<  1 1 1    »    CSQRTICUI  »»  21 

IF    I  IPO.GT.NI    GO    TO    5  22 

C  23 

DO  4  J* IPO.N  24 

U  *  1D»J  25 

C  26 

DO  3  M-1.IM0  27 

MO  «  IM-1 )*N-IM*M-MI/2  28 

MI  »  MD»I  29 

MJ  «  MD»J  30 

3  C1IJJ    •    C(  IJI-C«MJ»»C(MN  31 


C 


C 


c 


4    CIUI    -    CUJI/CIIII 


K    =    N-I 


KL    =    KD*L 


DO    12    1*1. N 


PH  «    .0 
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II 


C  3* 

5  CONTINUE  35 
C  36 

6  Sill  «  SU)/CU)  37 
C  38 

DO  8  1*2, N  39 

I  MO  ■  I-i  40 

C  41 

00  7  1*1. IMO  42 
LI  •  IL-Il*N-<L«L-LI/2M  43 

7  SU1  ■  Sll)-CILI)*S(LI  44 
C  45 

1  I  -  €1-1 »*N-I 1*1  — 1I/2»I  46 

8  $11)  -  SUI/CUII  47 
C  48 


NN  *  (  IN»ll*N)/2  «? 

SINI  •  SINI/C1NNI  50 

NMO  *  N-l  |1 

DO  10  1*1, NMO  53 


5*. 


KPO  =  K*\  55 

KO  »  IK-l l*N-(K*K-KJ^2  56 

00  9  L=KPO,N  58 


59 


9  SIKI  *  SIK»-CIKLI*SILI  60 

C  61 

KK  *  KO»K  51 

10  SIKI  «  SIK1/CIKK)  63 


64 
IF  I IWR.LE.OI  GO  TO  13  £5 


CNOR  «  .0  fc6. 

C  ° ' 

00  11  1*1, N  £g 

SA  *  CABStSI III  ,                   *' 

11  IF  ISA.GT.CNOR)  CNOR-SA                                                      70 

C  U 


IF  ICNOR.LE.O.I  CN0R«1.  |2 
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SS  -  St  1 1  H 

SA  »  CABSISSI  T,6, 

SNOR  -  SA/CNOR  .'  J7. 
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IF  ISA. GT. 0.1  PH  •  57.29578*ATAN2IAIMAGISS),REAL(SSII  79 

12  WRITE  16,141  l.SNOft.SA.PH.SS  g° 

WRITE  (6,15)                               ,  *{ 

13  RETURN  2* 

84 

14  FORMAT  UX,1I5,1F10.3,IF15.7,1F10.0,2F15.6>  85 

15  FORMAT  I1H0I  g6, 
END  B' 
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•SQROT 
CALLING    SEQUENCE    OF    THE    SUBROUTINES 
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SYMBOL  DICTIONARY 


A 

ABAP 
ABAT 
ABIP 

ABIT 

ACSP 
ACST 

AFFP 
AFFT 
All 

BM 

C 

CG 
CGD 
CJ 

CMM 
D 

ECSP 
ECST 

EFF 

EP 

EPP 

EPPS 

EPTS 

EP2 
EP3 
EP4 
ERR 
ER2 


characters  of  the  input  data  cards 

backscattering  phi  plane  angle  for  plotting 

backscattering  theta  plane  angle  for  plotting 

bistatic  scattering  phi  plane  angle  for 
plotting 

bistatic  scattering  theta  plane  angle  for 
plotting 

absorption  cross  section  for  phi  polarization 

absorption  cross  section  for  theta 
polarization 

far-zone  phi  plane  angle  for  plotting 

far-zone  theta  plane  angle  for  plotting 

radius  of  the  thin  wire  of  the  structure 

outer  radius  of  the  dielectric  shell  of  the 
insulation  of  the  wire 

elements  of  the  open-circuit  impedance  matrix 

branch  currents  for  the  structure 

cosh  /d  for  a  given  segment 

loop  currents  for  the  structure 

conductivity  of  the  wire 

length  of  a  given  segment 

extinction  cross  section  for  phi  polarization 

extinction  cross  section  for  theta 
polarization 

radiation  efficiency 

loop  currents  induced  by  a  phi  polarized  wave 

phi-polarized  far-zone  feild  of  the  dipole 
mode 

scattered  electric  field  in  the  phi  direction 
due  to  a  phi  polarized  wave 

scattered  electric  field  in  the  theta 
direction  due  to  a  phi  polarized  wave 

complex  permittivity  of  insulation 

complex  permittivity  of  ambient  medium 

complex  permittivity  of  ground 

EP4/EP3 

relative  dielectric  constant  of  insulation 
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ER3       relative  dielectric  constant  of  the  ambient 
medium 

ER4       relative  dielectric  constant  of  the  ground 

ET        loop  current  induced  by  a  theta  polarized 
wave 

ETA       intrinic  impedance  of  ambient  medium 

ETPS      scattered  electric  field  in  the  phi  direction 
due  to  a  theta  polarized  wave 

ETT       theta  polarized  far-zone  field  of  the  dipole 
mode 

ETTS      scattered  electric  field  on  the  theta 
direction  due  to  theta  polarized  wave 

EX  near-zone  electric  field  in  x  direction 

EY  near-zone  electric  field  in  the  y  direction 

EZ  near-zone  electric  field  in  z  direction 

EO  8.854E-12 

FHZ  frequency  in  hertz 

FMC  frequency  in  megahertz 

GAM       intrinic  progration  constant  of  the  ambient 
medium 

GG        time-average  power  input 

GPP       power  gain  associated  with  the  phi  polarized 
component 

GTT       power  gain  associated  with  the  theta 
polarized  component 

HGT  height  of  the  structure  above  ground  plane 

IA  first  node  of  a  given  segment 

IB  second  node  of  a  given  segment 

IBISC  indicator  for  bistatic  scatter  calculations 

ICARD  indicator  for  the  data  cards 

ICJ       dimension  corresponding  to  the  number  of 
simultaneous  linear  equations 

IFLAG     indicator  for  program  termination 

IGAIN     indicator  for  antenna  gain  calculations 

IGED      indicator  for  presence  of  the  ground  plane 

INC       indicator  for  the  type  of  far-zone 
calculations 

INEAR     indicator  for  near-zone  calculations 

INM       dimension  corresponding  to  the  number  of 
m on  poles 

INT       number  of  integration  steps 
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ISC  indicator  for  the  insualtion 

ISCAT  indicator  for  backscatter  calculations 

IWR  indicator  for  current  distribution  output 

11  endpoint  node  of  a  given  dipole 
112  indicator  for  auxiliary  matrix 

12  terminal  node  number  of  a  given  dipole 

13  endpoint  node  number  of  a  given  dipole 
JA  first  segment  number  of  a  given  dipoile 
JB  second  segment  number  of  a  given  dipole 
KFLAG  print  indicator 

KGEN  list  of  generator/feed  locations 

LOAD  indicator  for  structure  load 

LZD  list  of  impedance/load  locations 

MAX  maximum  of  the  number  of  segments  connected 
to  any  one  given  node 

MD  list  of  dipoles  sharing  a  given  segment 

MIN  minimum  of  the  number  of  segments  that 
connected  to  any  one  given  node 

MSG  indicator  for  error  printout 

N  number  of  simultaneous  linear  equations 

ND  total  number  of  dipoles  sharing  a  given 
segment 

NGEN  indicator  for  antenna  calculations 

NM  number  of  segments 

NPL  indicator  for  polar  plot 

OMEGA  angular  frequency 

PH  phi  angle  for  far-zone  calculations 

SCSP  scattering  cross  section  for  phi  polarization 

SCST  scattering  cross  section  for  theta 
polarization 

SGD  sinh  /d  of  a  given  segment 

SIG2  conductivity  of  insulation 

SIG3  conductivity  of  the  ambient  medium 

SIG4  conductivity  of  ground 

SPPM  echo  area  phi  incident-phi  scattered  wave 

SPTM  echo  area  phi  incident-theta  scattered  wave 

STPM  echo  area  theta  incident-phi  scattered  wave 

STTM  echo  area  theta  incident-theta  scattered  wave 
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TD2  loss  tangent  of  the  insulation 

TD3  loss  tangent  of  ambient  medium 

TH  theta  angle  for  far-zone  calculations 

TP  2n  (6.28318) 

UO  1.2566E-6 

VG  antenna  complex  driving  voltages 

VOLT  list  of  VG's 

X  x-coordinate  of  each  node 

XNP  list  of  XP's 

XP  x-coordinate  for  near-zone  calculations 

Y  y-coordinate  of  each  node 

YNP  list  of  YP's 

YP  y-coordinate  for  near-zone  calculations 

Y11  complex  power  input 

Z  z-coordinate  of  each  node 

ZLD  complex  load  at  a  given  node 

ZLLD  list  of  ZLD's 

ZNP  list  of  ZP's 

ZP  z-coordinate  for  near-zone  calculations 

ZS  surface  impedance  of  the  wire 

Z11  antenna  input  impedance 
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APPENDIX  B 
USER'S  MANUAL 
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USER'S  MANUAL 

The  Antennas-Scatterers  Analysis  Program  (ASAP)  for  thin 
wire  structures  in  a  homogenous  conducting  medium  performs  a 
frequency  domain  analysis  of  antennas  and  scatters.  The 
program  is  applicable  in  the  presence  of  either  a  perfect  or 
a  finite  ground.   This  appendix  will  describe  and  explain 
the  data  cards  necessary  to  execute  the  compute  program. 
Although  the  program  was  written  for  the  IBM  360  computer 
system,  it  can  be  executed  on  another  system  with  minor 
modifications. 

The  program  utilizes  piecewise  sinusoidal  expansion  for 
the  current  distribution  with  Kirchhoff  Current  Law  enforced 
everywhere  on  the  structure.  If  the  structure  contains  end 
points,  the  currents  at  these  points  are  assumed  to  vanish. 

I.  Program  Limits 

The  thin  wire  assumptions  are  questionable  and  the 
accuracy  and  convergence  deteriorate  if  the  radius  of  wire 
utilized  for  the  structure  exceeds  0.01   of  a  wavelength,  if 
the  longest  segment  is  greater  than  one-fourth  of  a 
wavelength,  if  the  length  ratio  of  the  longest  and  shortest 
segments  exceeds  100,  or  if  the  total  wire  length  is  less 
than  30  times  the  wire  diameter.   If  a  wire  is  bent  sharply 
to  form  a  small  acute  angle  (less  than  30  degrees)  ,  the  thin 
wire  model  is  questionable.  It  is  assumed  that  the  wire 
conductivity  greatly  exceeds  the  conductivity  of  the  ambient 
medium,  for  insulated  wires,  the  dielectric  layer  is  assumed 
to  be  electrically  thin. 

II.  Minimum  Data 

The  minimum  data  necessary  to  execute  the  program  is: 
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a.  description  of  structure 

b.  radius  of  wire  used  for  the  structure 

The  program  will  default  to  the  other  parameters  necessary. 
The  default  parameters  are: 

a.  wire  for  the  structure  is  copper 

b.  frequency  of  operation  is  300  mhz 

c.  homogneous  medium  is  free  space 

A  more  detailed  explanation  of  the  defaults  will  be 
discussed  when  the  data  card  for  the  parameter  is  described. 

III.  Outputs 

In  antenna  problems,  the  output  includes  structure 
currents,  impedance(s)  of  f eed (s) ,  gain,  polar  radiation 
plots,  and  near  field  calculations.  In  bistatic  scattering 
problems,  the  output  includes  structure  currents,  complex 
elements  of  the  polarization  scattering  matrix,  polar 
reradiation  pattern  plots,  and  echo  areas  produced  by  a 
plane  wave.   For  backscattering  problems  the  output  includes 
absorption,  scattering  and  extinction  cross  sections  in 
addition  to  the  outputs  of  bistatic  scattering.  Most  of  the 
outputs  are  suppressed  and  must  be  requested.  Since  the 
program  can  produce  a  large  volume  of  output,  care  should  be 
exercised  until  the  user  is  familiar  with  the  outputs. 

IV.  Data  Cards 

The  Analysis  Program  utilizes  free  format  for  the  data 
cards,  that  is,  the  program  utilizes  character  recognition 
to  determine  which  parameters  are  being  read.  Data  placement 
(location)  on  the  input  card  is  not  critial.  Blank 
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characters,  on  all  input  cards  but  the  COMMENT  data  card, 
are  ignored  and  may  be  used  at  the  discretion  of  the  user. 
Since  character  recognition  is  used,  only  the  first  four 
characters  of  the  key  words  roust  be  present  and  correct. 

The  format  for  the   COMMENT  CARD  utilizes  standard 
FORTRAN  format  (i.e.  ■Cl  in  column  1  followed  by  at  least 
four  blanks) .   The  COMMENT  CARD  is  the  only  type  of  input 
card  that  position  in  the  data  block  is  critial.   This 
(these)  card(s)  must  be  placed  at  the  beginning  of  a  data 
block.  A  data  block  is  a  series  of  related  data  cards. 
Several  data  blocks  may  be  used  to  define  an  analysis 
problem.  This  will  become  clear  when  the  termination  cards 
(END,  STOP,  or  CHANGE)  are  discussed.   There  is  no  limit  to 
the  number  of  comment  cards  that  may  be  used.  As  a  check  for 
the  user,  all  input  data  cards  will  appear  on  the  output  as 
they  appear  in  the  input  deck. 

The  format  of  other  data  can  be  of  one  of  two  forms: 

a.  type  of  card  (option  1/option  2/ ) 

b.  parameter  (value)  . 

The  type  of  format  to  use  will  be  apparent  as  the  individual 
data  cards  are  discussed. 

The  numerical  values  for  the  parameters  may  be  stated  in 
any  one  of  the  following  forms.  The  program  will  translate 
the  number  to  the  proper  form  for  the  specified  parameter, 
either  fixed  or  floating  point.   All  of  the  following 
examples  have  the  same  value. 

0.0001  or  .0001  or  100. U  or  100U  or  .1M  or  0.1M  or  .0000001K 

-6  -3  3 

U  =  10  M  =  10  K  =  10 
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1"  Jill!  This  card  is  used  to  define  the  parameters 
associated  with  the  wire  utilized  by  the  thin  wire 
structure.  Two  options  are  available  and  are  defined  as: 

RADIUS=value  of  the  radius  of  the  wire  in  meters 

COND0CTIVITY=value  in  megamhos  per  meters  . 

The  wire  data  card  must  appear  in  the  first  data  block  to 
define  wire  radius.   The  default  value  of  the  conductivity 
is  50  megamhos/meter  (copper) . 


4»IRE(   RADIUS=.001/   CONDU=28.5) 
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2.  INSOLATION  This  card  is  utilized  to  define  the  parameters 
associated  with  the  insulation  of  the  wire  used  for  the 
structure  to  be  analyzed.  If  this  card  is  omitted,  the 
program  assumes  that  the  structure  is  uninsulated.  Four 
options  are  available  and  are  defined  as: 

RADIUS=value  of  outer  radius  in  meters 

CONDUCTIVITY=value  in  micromhos  per  meter 

DIELECTEIC=value  of  relative  dielectric  constant 

LOSS  TANGENT=value  . 

The  conductivity  and  either  the  relative  dielectric  constant 
or  the  loss  tangent  (but  not  all  three)  options  may  be 
stated. 


'INSULATION (  RADIUS=.015/  C0ND=7. /DIEL=5) 
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3.  EXTERIOR  MEDIUM  This  card  is  utilized  to  describe  the 
homogeneous  medium  surrounding  the  structure.  If  the  medium 
is  free  space,  this  card  may  be  omitted.   Three  options  are 
available  and  are  defined  as: 

DIELECTRIC=value  of  relative  dielectric  constant 
CONDUCTIVITY=value  in  micromhos  per  meter 
LOSS  TANG2NT=value  . 

As  with  INSULATION   card  state  either  conductivity  or  loss 
tangent. 


:XTE(LOSS=.45) 
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4.  DESCRIPTION  This  card  is  utilized  to  describe  the  shape 
of  the  wire  structure  to  the  program.  The  user  must  divide 
the  wire  structure  into  segments  of  the  appropriate  length 
and  number  each  node  starting  at  one.  A  node  is  a  point 
where  a  segment  begins  or  ends.  A  maximium  of  four  segments 
can  meet  at  any  given  node.   An  isolated  wire  must  contain 
at  least  two  segments  and  three  nodes.   Thus  the  DESCRIPTION 
CARD  must  show  at  least  3  consecutive  nodes  for  all  portions 
of  the  wire  structure.   The  structure  is  described  by 
stating  the  node  numbers  that  each  segment  connects.  The 
description  of  a  sguare  loop  might  appear  as; 

DESCRIPTION (1-2/2-3/3-4/4-1)  . 

The  description  of  a  dipole  and  reflector  might  appear  as: 

DESCRIPTION (1-2/2-3/3-4/4-5/6-7/7-8/8-9/9-10)  . 

If  the  description  will  not  fit  on  one  data  card  continue  on 
the  next  card  as  if  the  previous  card  were  longer.  The 
dipole  example  might  appear  as: 

DESCRIPTION ( 1-2/2-3/3-4/4-5/ 

6-7/7-8/8-9/9-10)  . 

Note  that  the  last  character  on  the  card  to  be  continued  is 
a  slant  (/) .  As  many  cards  as  necessary  may  be  used.  The 
maximum  number  of  nodes  permitted  is  fifty.  If  ground  plane 
is  present,  the  maximum  number  is  twenty-five.  If  a  ground 
plane  is  present  and  the  structure  touches  the  ground  plane, 
the  lowest  node  numbers  MUST  be  used  for  the  touching  nodes. 
That  is,  if  the  structure  touches  the  ground  plane  at  two 
points,  node  numbers  1  and  2  MUST  be  assigned  to  these 
nodes. 


DESCR (1-2/2-3/3-4/4-1) 
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5.  GEOMETRY  This  card  is  used  to  state  the  physical  location 
in  rectangular  coordinates  of  each  node  of  the  DESCRIPTION 
CARD  .  The  rectangular  grid  is  in  units  of  meters.  If  node  1 
is  located  at  x1,y1,z1  and  node  2  at  x2,y2,z2  and  node  3  at 
x3,y3,z3,etc. ,  the   GEOMETRY  CARD   might  appear  as: 

GEOMETRY  (x 1 , y 1 , z1/x2, y2, z2/x3, y 3,  z3/ ) 

As  with  the   DESCRIPTION  CARD,  continuation  cards  are 
permitted. 


'GEOM^IrO,.  1/-.  1,0.1/-.  1,0-.1/.1,0,-.1) 
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6.  FEJ?D  For  antenna  analysis  the  feed  point  (s)  and 
voltage (s)  must  be  stated.  In  the  foremention  dipole  and 
reflector  example  if  the  feeds  were  at  node  2  with  a  voltage 
source  of  .5  at  an  angle  of  -90  degrees  and  at  node  4  with  a 
voltage  source  of  .5  at  an  angle  of  +90  degrees  the   FEED 
CARD   might  appear  as: 

FEED (2,. 5, -9 0/4,. 5, +  9  0) 

The  order  of  the  information  for  each  voltage  source  is  node 
number,  magnitude,  and  phase  angle.  This  order  is  repeated 
until  all  sources  are  stated.  If  the  source  information  will 
not  fit  on  one  card,  use  another  card  similiar  to  the 
initial  one;  that  is,  repeat  the  word  "FEED".   If  only  one 
voltage  source  is  applied  to  the  structure,  only  the  node 
number  must  be  stated.  In  the  dipole  example,  if  the  drive 
is  at  node  3,  the   FEED  CARD   might  appear  as: 

FE£D(3) 

A  default  source  of  one  volt  at  zero  degree  phase  is 
assumed.  Voltage  sources  should  only  be  stated  for  nodes 
with  only  two  segments. 


FEED  (2,. 5, -9 0/4,. 5, +90) 
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7.  LOAD  This  card  is  used  to  describe  the  loads  to  be  placed 
at  various  locations  on  the  structure.  The  format  for  this 
card  is  similiar  to  that  of  the   FEED  CARD,  that  is,  the 
word  "LOAD"  is  used  in  the  place  of  "FEED".   The  order  of 
the  information  on  the  card  is  the  same.   Since  this  card  is 
frequency  dependent,  it  must  be  changed  if  the  frequency  of 
operation  is  changed.  No  default  parameters  are  available. 
The  structure  is  assumed  unloaded  unless  this  card  is  used. 
Once  the  structure  is  loaded,  it  will  remain  loaded  for  the 
remainder  of  the  data  block  series.  To  unload  the  structure 
the  following  card  may  be  used: 

L0AD<-1) 


'LOAD  (1,  120,-45/3,1 20,  *H  5) 
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8.  OUTPUT  This  card  is  used  to  request  output  data.  Most  of 
the  output  is  in  tabular  form.   More  than  one   OUTPUT  CARD 
is  permitted  per  data  block,  but  not  for  the  same  type  of 
output.   If  only  the  antenna  input  impedance,  antenna 
efficiency,  or  time-average  power  input  is  of  interest,  no 
OUTPUT  CARD  is  necessary.  These  parameters  are  automatically 
printed  if  a   PEED  CARD   or  GENERATOR  CARD   is  utilized.  One 
or  more  of  the  following  options  may  be  used  to  request  the 
various  outputs  available. 

FAR  FIELD=phi  initial,  phi  final,  theta  initial,  theta 
final 

This  option  gives  the  components  of  the  electric  field 
intensity  in  the  far  field  as  phi  and  theta  varies  between 
limits  specified  in  one  degree  divisions. 

BACKSCATERING=phi  initial,  phi  final,  theta  initial, 
theta  final 

This  option  gives  the  absorption,  scattering,  and  extinction 
cross  sections,  and  the  complex  elements  of  the  polarization 
scattering  matrix  for  an  incident  plane  wave  illuminating 
the  structure  from  the  spherical  direction  of  phi,  theta  as 
both  vary  between  limits  specified  in  one  degree  divisions. 

BISTATIC=phi  inital,  phi  final,  theta  initial,  theta 
final 

This  option  gives  echo  area  and  the  complex  elements  of  the 
polarization  scattering  matrix  for  an  incident  plane  wave 
illuminating  the  structure  from  the  spherical  direction  phi, 
theta  final  of  the  backseat tering  output  option,  reradiated 
in  the  phi,  theta  direction  as  both  vary  between  limits 
specified  in  one  degree  divisions.  A  bistatic  output  request 
must  be  accompanied  with  a  backscattering  request  in  the 
same  data  block. 

STEP=value  in  degrees 
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This  option  will  cause  any  of  the  above  output  options  to  be 
stepped  at  a  different  interval  size.  That  is,  if  one  of  the 
above  options  is  to  be  stepped  at  ten  degrees  intervals,  use 
this  option.  This  option  overrides  the  one  degree  stepping. 

CURRENT 

This  option  gives  the  currents  on  the  structure  which  are 
produced  by  the  feed/generator  voltages  and/or  the  incident 
plane  wave  of  the  backscattering  reguest. 

NEAR=x1,y1,z1 

or 

NEAR=(x1,y1,z1/x2,y2,z2/x3,y3,z3/etc ) 

This  option  gives  the  value  of  electric  field  components  in 
the  near  field  for  the  antenna  at  the  point  or  points 
specified. 


/OUTPUT  (FARF=45,50,25,50) 
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9-  PLOT  This  card  will  produce  normalized  polar  plots  in  the 
specifed  plane  for  the  stated  option.   The  plane  is 

specified  by  stating  either  "PHI= "  or  "THETA= " .  The 

PLOT  CARD  overrides  the  limits  of  the   OUTPUT  CARD  for  the 
same  option.  If  only  a  normalized  pattern  is  of  interest, 
only  a   PLOT  CARD  is  necessary.  If  a  table  of  values  and  a 
normalized  pattern  is  desired,  both  a   PLOT  CARD  and  OUTPUT 
CARD   must  be  used.   Only  one   PLOT  CARD   is  permitted  per 
data  block.   The  following  pattern  plots  are  available: 

FAR  FIELD/plane 

This  option  will  plot  the  far  field  intensity  for  each 
component  of  the  electric  field. 

BACKSCATTERING/plane 

This  option  will  plot  the  normalized  magnitude  of  each  of 
the  elements  of  the  polarization  scattering  matrix. 

BISTATIC/plane 

This  option  will  plot  the  normalized  magnitude  of  each  of 
the  elements  of  the  polarization  scattering  matrix  produced 
by  the  incident  plane  wave  stated  by  final  limits  of  the 
backscattering  option  of  the  output  reguest. 


PLOT  (FARF/THET=90) 
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10.  GROUND  This  card  is  used  to  describe  the  ground 
parameters  if  a  ground  plane  is  present.  If  no  ground  plane 
is  present,  the  structure  is  assumed  to  be  in  free  space  or 
the  homogeneous  medium  of  the   EXTERIOR  MEDIUM  data  card. 
Seven  options  are  available  and  are  defined  as: 

PERFECT 

This  option  will  analyze  the  structure  over  a  perfect  ground 
plane. 

GOOD 

This  option  will  analyze  the  structure  over  a  good  ground 
plane  where  the  conductivity  of  the  ground  is  .02  mhos/meter 
and  the  relative  dielectric  constant  is  30. 

POOR 

This  option  will  analyze  the  structure  over  a  poor  ground 
plane  where  the  conductivity  of  the  ground  is  .001 
mhos/meter  and  the  relative  dielectric  constant  is  4. 

SEA 

This  option  will  analyze  the  structure  over  salt  water  where 
the  conductivty  of  the  water  is  4.  mhos/meter  and  the 
relative  dielectric  constant  is  80. 

HEIGHT=value  in  meters 

This  option  will  analyze  the  structure  with  origin  of  the 
GEOMETRY   card  this  height  above  the  ground  plane.  The 
lowest  point  of  the  structure  must  not  lie  below  the  ground 
plane.  It  may  lie  on  the  ground  plane. 

CONDUCTIVITY=  value  in  mhos/meter 

This  option  is  used  to  state  the  value  of  conductivity  of 
the  ground  plane  if  the  default  values  mentioned  above  are 
not  utilized. 
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DIELECTBIC=  value 

This  option  is  used  to  state  the  relative  dielectric 
constant  of  the  ground  plane  if  the  default  values  mentioned 
above  are  not  utilized. 


'GROUND  (HEIG=10/COND=.002/DIEL=10) 
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11.  INTERVAL  FOR  CALCULATION  This  card  is  used  to  state  the 
number  of  intervals  to  be  used  for  calculating  the  elements 
of  the  impedance  matrix  with  Simpsons-rule  integration.  A 
large  value  for  the  number  improves  the  accuracy  at  the 
expense  of  greater  execution  time.  For  most  problems  a 
suitable  combination  of  speed  and  accuracy  is  obtained  with 
a  value  of  four,  the  default  value.  If  the  rigiorous 
closed-form  impedance  expressions  in  terms  of  the 
exponential  integrals  is  desired,  set  this  value  to  zero. 

INTEBVAL=value 


'lNTE(6) 
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12.  GENERATOR  This  card  is  similiar  to  the  FEED  CARD  in  use, 
except  that  the  segment  numbers  are  stated  instead  of  the 
node  numbers.  This  is  useful  if  three  or  four  segments  meet 
at  a  node.  The  positive  terminal  of  the  generator  is 
connected  to  the  specified  segment  such  that  current  is 
forced  in  the  the  positive  direction.  The  positive  direction 
of  current  flow  is  from  the  first  stated  node  number  of  that 
segment  toward  the  second  stated  as  ordered  on  the 
DESCRIPTION  CARD. 


/GENE(2r.5,-9  0/Ur  .5, +9  0) 
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13.  IMPEDANCE  This  card  is  similiar  to  the  LOAD  CARD  in  use, 
except  that  the  segment  numbers  are  stated  instead  of  the 
node  numbers.  As  with  the  GENERATOR  CARD,  this  is  used  if 
three  or  four  segments  are  connected  to  a  node.   The 
impedance  will  be  connected  to  the  positive  terminal  of  the 
specified  segment.  The  format  of  this  card  is  the  same  as 
the  LOAD  CARD. 


IMPE (1,120,-45/3, 120, +4  5) 


14.  FREQUENCY  This  card  is  used  to  specify  the  frequency  in 
megahertz  if  it  is  to  be  other  than  the  default  value  of  300 
MHz. 


FREQ(12.5) 
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15.  CHANGE  This  card  at  the  end  of  the  data  block  signals 
the  program  that  the  following  data  cards  are  changes  to  the 
previously  read  data,  for  the  next  run.  If  a  "CHANGE  CARD" 
is  used,  the  outputs  must  be  requested  again  in  the  next 
data  block.   The  "CHANGE  CARD"  cannot  be  used  to  change 
"DESCRIPTION  CARD"  or  "GEOMETRY  CARD"  data  when  operating 
with  a  "GROUND  CARD".  Use  an  "END  CARD"  to  make  changes  when 
a  "GROUND  CARD"  is  used. 


16.  END  This  card  signals  the  program  that  this  is  the  end 
of  a  data  block  series  and  to  reinitialize  data  for  the  next 
problem.   An  "END  CARD"  cannot  be  used  with  a  "CHANGE  CARD". 


17.  STOP  This  card  signals  the  program  that  all  of  the  data 
cards  have  been  read  and  to  terminate  itself  when  execution 
is  completed.   This  card  must  be  used  as  the  last  card  in 
place  of  the  "END  CARD"  of  the  last  data  block  series.   A 
"STOP  CARD"  cannot  be  used  with  an  "  END  CARD"  in  the  same 
data  block. 
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C      AN  EXAMPLE  PROBLEM 

C 

C      V  ANTEMNA 

C 

WIRE(RADIUS=1M) 

GEOMiOt-. 18, +.18/0, -.09, *. 09/0, 0,0/0, 0.09,. 09/0,. 18,. 18) 

DESCU-2/2-3/3-4/4-5) 

FEEDO) 

OUTPUT(FARF=4  5,50,65,80/STEP=5) 

CHANGE 

OUTPUT(BIST=45,45,45,^5/BACK=0,0,10,12) 

OUTPUT(CURRENT) 

CHANGE 

C 

C  CHANGE    STRUCTURE    SHAPE    TO    DIPOLE 

C 

GEOM(0,-. 25, 0/0, -.125, 0/0, 0,0/0, .125, 0/0, .25,0) 

PLOT(FARF/PHI=90> 

GROUND(HEIGHT=.25/G000) 

STOP 


THE  ABOVE  DATA  DECK  WILL  PRODUCE  THE  OUTPUT  ON  THE 
FOLLOWING  PAGES. 
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